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INTRODUCTION 


The  tc^ic  of  the  experimental  research  performed  under  Air  Force  contract  was 
magnetic  field  line  reconnection  in  a  large  laboratory  plasma.  The  motivation  was 
to  test  theoretical  models  and  ideas  a^^lied  to  reconnection  in  space  pleianas. 

Particular  attention  has  been  focused  on  the  magnetic  field  topology  which  is  thought 
to  change  by  tearing  instabilities  resulting  in  impulsive  conversion  of  magnetic 
energy  to  kinetic  particle  energies. 

The  first  report,  Reference  1,  describes  measurements  and  interpretations  of 
triggered  tearing  modes  in  a  neutral  sheet.  In  these  experiments  multiple  magnetic 
X-  and  O-points  were  imposed  on  the  neutral  sheet  and  growth/decay  of  these  structures 
were  followed  in  space  and  time.  No  major  disruptions  were  observed. 

Reference  2  sunmarizes  the  nonlinear  interactions  between  fields  and  particles 
in  a  neutral  sheet.  Fully  resolved  particle  distribution  functions  in  three-dimensional 
velocity  space  have  been  measured.  Likewise,  fluctuations  in  the  magnetic  field 
have  been  resolved  by  three-dimensional  correlation  measurements.  As  shown  in  detail 
in  Reference  3,  these  fluctuations  consist  of  oblique  whistler  waves  which,  as  demon¬ 
strated  in  Reference  4  are  driven  unstable  by  anisotropic  tails  in  the  electron 
distribution  function.  These  tails  are  generated  by  particle  acceleration  in  the 
reconnection  electric  field  along  the  separator,  i.e.  by  conversion  of  magnetic  field 
energy  via  Faraday's  law. 

Reference  5  describes  experiments  on  current  sheet  disruptions  such  as  those 
which  may  arise  in  magnetic  substorms  or  solar  flares.  Most  striking  is  the  observation 
of  the  formation  of  double  layers  in  the  current  path.  At  these  double  layers,  the 
meignetic  field  energy  is  converted  to  particle  beams  which,  upon  injection  into  the 
surrounding  plasma,  generate  microinstabilities  and  heat  electrons  euid  ions. 

Reference  6  summarizes  our  understanding  of  reconnection  in  a  laboratory  plasma 
and  relates  the  finding  to  observations  in  space.  It  is  pointed  out  that  MHD  theories 


cannot  account  for  the  plasnta  dynamics  near  null  points  since  the  ions  are  essentially 
unmagnetized  and  the  magnetic  f ields/currents  are  mainly  controlled  ty  the  electrons. 
For  example,  magnetic  field  transport  occurs  by  electron  vrfiistler  waves  rather  than 
ionic  Alfven  waves. 

Due  to  space  charge  coupling  of  electrcais  euid  ions,  an  electron  current  system 
can  be  drastically  modified  by  the  ion  dynamics  such  as  demonstrated  in  Reference  7. 

It  is  observed  that  strong  current  flow  in  narrow  flux  tubes  leads  to  space  charge 
electric  fields  which  punp  the  ions/plasma  out  of  the  flux  tube  leading  to  a  current 
disruption.  Such  phenomena  may  account  for  the  frequently  observed  turbulence  and 
fine  scale  structure  of  currents  in  neutral  sheets. 

StMIARY 

In  sunmary,  extensive  data  on  magnetic  field  line  reconnection  have  been  obtained 
in  a  laboratory  plasma  which  have  helped  to  understand  reconnection  processes  close  to 
magnetic  null  regions.  These  regions  are  difficult  to  study  in  the  magnutotail  due 
to  the  variability  of  the  solar  wind,  and  they  have  not  been  studied  theoreticcLLly  due 
to  the  difficulty  to  use  kinetic  theory  coupled  with  electromagnetism.  We  think  the 
laboratory  effort  contributed  to  the  understanding  of  reconnection  in  a  parameter 
regime  not  accessible  by  other  means. 
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Triggered  Tearing  Modes  in  a 
Magnetic  Neutral  Sheet 

R.  L.  Stenzel  and  W.  Gekelman 


PPG  756 


November  1983 


Department  of  Physics 
University  of  California 
Los  Angeles,  CA  90024 


Abstract 

An  experimental  investigation  of  tearing  mode  instabilities 
in  a  magnetic  neutral  sheet  has  been  performed.  A  long  flat 
current  sheet  is  generated  in  a  large  laboratory  plasma.  It  is 
perturbed  by  pulsing  a  current  through  a  fine  wire  along  the 
separator  so  as  to  trigger  the  formation  of  magnetic  0-and 
X-points.  Two  magnetic  islands  are  excited.  Their  growth, 
propagation  and  merging  within  the  neutral  sheet  is  observed. 
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Introductlon 

Tearing  mode  instabilities  are  considered  to  be  responsible 
for  the  dli,..uptlon  of  current  sheets  In  solar  flares  and 
■agne tospheric  substoras  (Blrn  and  Schindler,  1981;  Hones,  1980). 
Although  tearing  nodes  have  been  investigated  theoretically  In 
great  detail  (Furth  et  al.,  1963;  Coronlti,  1977;  Park  et  al., 
1983)  their  ezperlaental  Investigation  has  been  aostly  restricted 
to  Tokaaaks  (Von  Goeler  et  al.,  1974;  Robinson  and  McGuire, 

1979).  These  toroidal  devices  have  magnetic  field  topologies 
substantially  different  froa  those  In  space  and  In-sltu 
aeasureaents  of  the  magnetic  flelda  are  not  possible  due  to  high 
teaperatures  and  densities.  In  the  UCLA  reconnection  experiments 
a  long  flat  reverse  field  geometry  has  been  established  and 
diagnosed  In  detail  (Stenzel  et  al.,  1983a;  Gekelman  et  al., 
1982).  The  current  sheet  Is  expected  to  be  unstable  to  tearing 
and  indeed,  the  spontaneous  formation  of  magnetic  Islands  has 
been  observed  In  a  certain  parameter  range.  In  order  to  study 
tearing  modes  In  a  controlled  fashion  It  Is  desirable  to  perturb 
a  marginally  unstable  neutral  sheet  with  an  external  trigger 
pulse  and  then  to  observe  Its  space-time  evolution.  Such 
experiments  will  be  described  here. 

Experiment 

A  linear  discharge  plasma  column  shown  In  Fig.  1  Is 

generated  with  aim  diameter  cathode.  Typical  plasma  parameters 

12  -3 

are  a  density  n  »  10  cm  ,  temperature  kT  s:  10  eV,  kT.  s  leV, 

-4 

Argon  gas  pressure  p  -  2  x  10  Torr,  and  an  axial  confining  dc 
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magnetlc  field  lOG.  The  plasma  is  uniform,  collisionless 

and  highly  reproducible  In  pulsed  mode.  Typical  pulse  widths  are 

t  =5  msec  and  repetition  times  t  -  1  sec. 

P  r 

In  order  to  generate  a  current  sheet  a  pulsed  transverse 

magnetic  field  (x,y,t}  Is  applied.  It  Is  generated  by  pulsing 

axial  currents  (I  =  lOkA,  -  lOOpsec)  through  two  plates 

parallel  to  the  plasma  column.  The  magnetic  field  topology 

contains  In  vacuum  an  X-type  neutral  point  at  the  center. 

However,  In  the  presence  of  plasma  electron  currents  are  Induced 

which  are  so  large  (1^  ^  lOOOA)  Chat  Che  magnetic  field  topology 

Is  strongly  modified.  While  the  current  prefers  to  flow  In  the 

null  regions  of  the  transverse  magnetic  field  It  can  rearrange 

the  neutral  region  In  a  self-consistent  manner.  As  predicted  by 

MHO  theories  (Dungey,  1958)  we  observe  during  the  rise  of  the 

external  magnetic  field  the  formation  of  a  classical  neutral 

sheet  shown  In  Fig.  2a.  The  corresponding  current  sheet,  >  Vx 

has  a  thickness  ==  5  cm  which  lies  In  between  the 

collisionless  skin  depth  0.8  cm)  and  Che  ion  Larmor 

radius  (r  .  “  30  cm), 
cl 

Current  sheets  can  break  up  into  current  filaments  through 
Che  process  of  tearing  mode  Instabilities.  Fig.  2b  shows  the 
observation  of  a  large  magnetic  Island  which  has  grown  In  the 
original  neutral  sheet  at  larger  current  densities.  In  order  to 
study  in  a  controlled  manner  Che  growth  of  such  modes  the  current 
sheet  Is  kept  near  threshold  for  tearing  and  a  trigger  magnetic 
Island  Is  excited.  As  shown  In  Fig.  3  this  is  accomplished  by 
inserting  a  fine  wire  (0.5  mm  diam.)  parallel  to  the  neutral 
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sheet  along  the  full  length  of  the  device  and  applying  a  current 
pulse  to  It  (900  A,  20ysec).  The  wire  Is  Insulated  from  the 
plasma  and  causes  only  a  small  local  density  perturbation. 

Time  and  space  resolved  measurements  of  the  magnetic  field 
topology  are  obtained  with  a  movable  magnetic  probe  pair  (B^, 

B^).  For  repeated  pulses  the  time  history  of  Sj_(t)  Is  digitally 
recorded  at  typically  300  different  positions  In  the  x-z  plane 
transverse  to  the  axial  currents  J  .  The  large  data  flow 

y 

requires  the  extensive  use  of  modern  digital  data  processing 
equipment  Including  fast  anolog~to-dlgltal  converters  (100  MHz, 

32  k  words),  on-line  minicomputer  (LSI  11/23),  an  array  processor 
(MAP  200),  and  a  data  link  to  a  large  off-slte  computer  (Cray  1) 
for  final  data  analysis  and  display.  Not  only  magnetic  probe 
signals  are  analyzed  In  this  manner  but  also  signals  from 
particle  detector,  Langmuir  probes  and  wave  detectors  (Stenzel  et 
al.,  1983b).  Fig.  4  shows  a  schematic  summary  of  the  diagnostics 
and  data  processing  setup. 

Triggered  Tearing  Modes 

At  a  time  t  -  30psec  after  the  generation  of  the  current 

sheet  the  trigger  pulse  Is  applied  to  the  axial  wire.  Its 

waveform  Is  that  of  a  half-sine  wave  of  duration  T  -  tt /u  =20iisec 

and  peak  current  I  ==  900  A  and  peak  voltage  V_  ■  L  dl/dt  <== 

max  max 

1000  V. 

A  sequence  of  vector  magnetic  field  maps  at  various  times  t 
Is  presented  In  Fig.  5  a-1 .  Just  after  turn-on  of  the  current 
pulse  (t  *  30.2usec,  Fig.  5a)  the  perturbation  magnetic  field 
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becomes  visible.  The  x-z  plane  of  measurement  Is  chosen 
Inadver tedly  too  high  with  respect  to  the  device  axis  or  trigger 
wire.  At  t  s  36usec  (Fig.  5b}  the  perturbation  magnetic  field 
stands  clearly  out  from  the  neutral  sheet.  The  point  of 
measurement  Is  at  the  center  of  each  vector  whose  length  scales 
linearly  with  the  local  field  strength.  The  data  points  are 
separated  by  Ax  -  Az  -  2  cm.  Due  to  the  rapid  variation  of 
near  the  wire  the  measured  field  appears  somewhat  random  but 
close  Inspection  shows  that  the  field  lines  circle  around  the 
wire  located  at  x  ^  33  cm,  z  s:  6  cm.  The  current  In  the  wire  flows 
normal  to  the  x-z  plane  In  the  same  direction  as  the  plasma 
current,  l.e.  toward  the  cathode. 

Fig.  5c  and  d  show  the  growth  of  the  trigger  magnetic  Island 
up  to  its  maximum  at  t  ^  40  (isec.  The  addition  of  the  large 
magnetic  Island  to  the  neutral  sheet  has  generated  two  magnetic 
X-polnts  at  s  17  cm,  x^  »  50  cm  and  2  ~  ^  growing 

magnetic  Island  tends  to  decrease  the  plasma  current  since,  by 
Lenz's  law  the  Induced  currents  are  opposite  to  the  applied 
current.  The  axial  plasma  current  Is  greatly  Impeded  by  the 
strong  transverse  magnetic  field  of  the  Island  -  10... 100  G). 

An  ^  X  ^  plasma  drift  leads  to  the  expulsion  of  plasma  radially 
away  from  the  wire.  At  typical  ion  velocities  v^  =  c^  =  5  x  10^ 
cm/sec  the  plasma  can  move  in  At  lOpsec  over  distances  Ar  5 
cm  which  leads  to  significant  plasma  losses  in  the  Island.  The 
X-points  shift  outward  with  the  expanding  plasma.  However,  since 
the  separatrices  at  the  X-pointa  intersect  at  nearly  90°  little 
current  flows  in  these  regions. 
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Besldes  Che  transverse  magnetic  field  there  is  also  an  axial 

magnetic  field  component,  =  10  G,  applied  to  the  plasma.  The 

12  -3 

peak  plasma  pressure  nkT  =  10  cm  x  lOeV  is  larger  chan  the 

magnetic  field  pressure,  nkT =  4  This  implies  a 

significant  modification  of  Che  axial  field.  Observations  show 

Chat  the  axial  field  does  not  completely  vanish  but  can  locally 

decrease  to  50Z  of  its  external  value.  Nevertheless,  the  field 

lines  around  Che  wire  are  helical  with  varying  pitch  like  in  a 

flux  rope.  The  current  may  actually  arrange  Itself  in  a 

force-free  configuration  (V  x  ^  II  Three-dimensional  field 

maps  [B  .  ,,  B  /  ,  and  B  /  .]  would  be  required  to 

'■  X  (x,y,x)*  y  (x,y,z)  a  (x,y,z)'  ^ 

conclusively  establish  Che  field  topology.  The  capability  for 
performing  such  measurements  is  presently  being  build  up.  It 
requires  a  fully  three-dimensional  probe  motion,  three  vector 
component  measurements  and  a  full-size  computer  (VAX  11/750)  for 
an  order  of  magnitude  larger  data  arrays. 

Fig.  Se  and  f  show  the  field  topology  during  the  decrease  of 
Che  current  applied  to  Che  wire.  At  this  time  the  inductive 
electric  field  has  changed  sign  and  Induces  plasma  currents  in 
the  same  direction  as  Che  wire  currents.  These  currents  flow 
near  the  neutral  points  whose  separatrix  now  intersects  at  angles 
different  from  90°  (V  x  B  0). 

Fig.  5g  and  h  are  the  most  interesting  cases  showing  the 
excitation  of  two  magnetic  Islands.  These  islands  are  induced  in 
the  current  sheet  in  response  to  the  time-varying  wire  current. 

It  is  Important  to  note  that  the  induced  plasma  current  flows  in 
two  current  channels  located  in  the  original  neutral  sheet  rather 
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than  In  a  coaxial  return  current  distributed  symae tr ically  around 
the  wire.  This  feature  indicates  the  Inherent  instability  for 
the  current  sheet  to  filament  or  to  tear. 

While  in  Fig.  5h  the  0-and  X-points  are  clearly  separated 

one  observes  in  Fig.  51  and  j  the  merging  of  both  null  points. 

The  current  in  the  axial  wire  has  reversed  sign  and  is  growing  in 

magnitude.  The  trigger  field  lines  are  reversed  from  those  of 

the  first  half  cycle.  The  resultant  field  topology  forms  neutral 

points  above  and  below  the  wire.  The  normal  component  B  in  the 

z 

original  neutral  sheet  is  enhanced  so  that  the  null  region 
vanishes . 

Finally,  Fig.  Sk  and  1  shows  the  change  in  the  field 
topology  during  the  decrease  of  the  trigger  current.  The  null 
region  above  the  trigger  wire  has  separated  into  two  X-polnts 
(Xj^  s  15  cm,  X2  -  30  cm,  z  s  20  cm).  Over  a  large  region 
around  the  wire  the  field  lines  are  linked  with  it  indicating  a 
substantial  plasma  current  parallel  to  the  wire  current.  It  is 
worth  noting  that  at  t  =  SOiisec  the  applied  current  through  the 
plate  electrodes  has  reached  its  maximum  yet  the  plasma  current 
has  reversed  the  applied  component.  Except  near  the 
boundaries  the  field  inside  the  plasma  is  established  self- 
consistently  by  V  x  ^  ■  ^o— * 
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Flg.  1  Schenatlc  picture  of  Che  experlmenCel  arrangement.  (a) 
Crosa-sectlonal  view  showing  parallel  plate  electrodes 
with  pulsed  currents  Ig  and  magnetic  field  lines 
without  plasma.  (b)  Side  view  of  the  device  with  main 
electrodes,  currents  (In,  Ig)  electric  fields  (I  --I  - 
V(^p),  and  magnetic  fields  (S  ~  Bj_  tyo)*  The  coordinate 
system  common  In  magnetospherlc  physics  has  been  adopted 
where  y  Is  along  the  neutral  line  (device  axis),  x  Is 
along  the  horlxonCal  neutral  sheet,  and  z  Is  normal  Co 
the  sheet. 

Fig.  2  Transverse  magnetic  field  topologies  during  field  line 

reconnection.  (a)  Magnetic  neutral  sheet  (b)  Magnetic 
Island  (0-polnt  with  adjacent  X-polnCs)  due  to  tearing 
Instability . 

Fig.  3  Schematic  experimental  arrangement  for  triggering  tearing 
modes.  k  fine  wire  Is  stretched  along  the  device  axis 
and  a  current  pulse  Is  passed  through  It.  It  excites 
magnetic  0-  and  X-polnts  In  the  current  sheet. 

Fig.  4  Block  diagram  of  Che  plasma  diagnostics  performed  with 
various  probes  and  a  digital  data  acquisition  system. 

Fast  mass  data  handling  with  computers  allows  one  to 
perform  time  and  space  resolved  vector  measurements, 
statistical  analysis  and  particle  distribution  function 
measurements . 

Fig.  5  Time  evolution  of  the  transverse  magnetic  field  topology 
.  In  response  to  Che  trigger  magnetic  Island. 

The  excitation  of  two  0-polnts  and  X**polnts  due  to 
current  sheet  tearing  Is  clearly  visible  In  Fig.  5h. 
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Stimulated  by  important  dynamic  processes  in  space  and  fusion  plasmas,  the  reconnection  of  magnetic  fields  in  plasmas 
is  investigated.  A  well-controlled  large  laboratory  plasma  is  established,  and  direct  observations  of  magnetic  fields  and  particle 
distributions  are  performed.  High  speed  computers  are  interfaced  with  the  experiment  so  as  to  permit  the  analysis  of 
multidimensional  functions.  These  include  the  average  magnetic  field  topology  <A(r,  /)>  and  the  tensor  correlation  function 
of  magnetic  fluctuations,  as  well  as  the  electron  phase  space  distribution  function  /(r,  r,  ().  The  observations  show 
that  the  characteristic  field  topology  of  a  neutral  sheet  is  subject  to  a  high  level  of  magnetic  turbulence  consisting  of  low 
frequency  whistler  waves.  These  are  excited  by  anisotropies  in  the  electron  distribution  functions  inside  the  current  sheet,  in 
particular  by  energetic  electrons  which  run  away  in  the  electric  field  along  the  separator.  The  observed  microinstabilities  and 
kinetic  processes  cannot  be  described  by  classical  fluid  theories  or  MHD  simulations  but  are  important  in  determining  the 
effective  resistivity  and  reconnection  rate. 


1.  Introduction 

Solar  flares  [1]  and  magnetic  substorms  [2] 
present  examples  from  space  plasma  physics  where 
a  strong  nonlinear  interaction  of  magnetic  fields 
with  plasmas  takes  place.  These  dynamic  processes 
are  both  characterized  by  a  rapid  energy  transfer 
from  stored  magnetic  field  energy  into  plasma 
kinetic  energy  [3].  The  events  appear  to  violate  the 
classical  MHD  picture  which  predicts  that  mag¬ 
netic  fields  are  essentially  “frozen  in”  [4]  highly 
conducting  fluids  and  cannot  be  dissipated  ext^t 
in  small  regions  near  magnetic  null  points  where 
finite  resistivity  permits  diffusion,  hence  recon¬ 
nection  of  field  lines  and  dissipation  [^.  The 
discrepancy  between  the  observed  impulsive  energy 
exchange  and  the  classical  steady-state  theories  has 
long  been  recognized,  and  models  for  various 
plasma  instabilities  near  null  points  have  been 
proposed  [6].  These  involve  either  macroscopic 
changes  in  the  magnetic  field  topology  by  tearing 
instabilities  [7]  or  microscopic  instabilities  which 
generate  small  scale  fluctuations  and  effectively 
enhance  the  plasma  resistivity  [8].  Either  case  leads 
to  nonlinear  phenomena  since  field  and  plasma 


properties  are  highly  coupled.  In  spite  of  many 
years  of  theoretical  work,  numerical  simulation 
effort  and  observations  from  spacecraft,  the  com¬ 
plex  reconnection  process  is  still  an  open  funda¬ 
mental  problem  in  plasma  physics. 

Its  significance  is  not  restricted  to  space  plasmas. 
In  various  fusion  devices  such  as  tokamaks,  reverse 
field  pinches,  spheromaks,  and  laser  fusion,  mag¬ 
netic  null  point  geometries  arise  [9],  In  most 
confinement  devices  there  is  evidence  that  tearing 
instabilities  lead  to  a  loss  of  confinement  of  par-: 
tides  and  energy  even  though  the  details  are  not 
known. 

A  common  problem  in  fusion  and  space  plasmas 
is  the  difficulty  of  making  direct  detailed  obser¬ 
vations  of  the  reconnection  processes  near  mag¬ 
netic  null  points.  The  basic  physical  processes  are 
best  observed  in  a  well-controlled  laboratory 
plasma  designed  for  optimizing  diagnostic  access. 
A  large  diameter  discharge  plasma  has  been  used 
at  UCLA  for  reconnection  studies  for  several  years 
[10-15]. 

Although  one  cannot  model  in  a  laboratory 
device  global  configurations  of  cosmic  plasmas,  the 
characteristic  magnetic  null  point  geometry  can  be 
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reproduced  [16-18],  Here  we  will  show  the  for¬ 
mation  of  a  classical  magnetic  neutral  sheet  [5],  It 
presents  the  interface  between  two  oppositely  di¬ 
rected  magnetic  field  regions,  and  it  is  maintained 
by  a  plasma  current  sheet.  Its  stability  has  been 
investigated  and  is  found  to  contain  a  high  level  of 
magnetic  turbulence.  It  consists  of  a  broad  spec¬ 
trum  of  unstable  whistler  waves.  The  source  of 
these  instabilities  lies  in  the  non-equilibrium  elec¬ 
tron  distribution  function  in  the  current  sheet. 
Anisotropies  in  the  form  of  runaway  electron  tails 
are  observed  from  direct  measurements  of 
/ (r,  r,  /).  The  significance  of  these  findings  to  the 
reconnection  physics  is  twofold:  (i)  Non- 
Maxwellian  distribution  functions  modify  the  clas¬ 
sical  transport  parameters;  for  example,  the 
effective  resistivity  is  determined  by  the  runaway 
electrons  rather  than  the  bulk  electron  “tem¬ 
perature”.  (ii)  Velocity  space  instabilities  create  a 
high  level  of  turbulence  which  interacts  with  reso¬ 
nant  particles.  These  wave-particle  interactions  are 
outside  the  models  described  by  MHD  theory,  yet 
may  explain  “anomalous”  resistivities,  enhanced 
reconnection,  and  energy  transfer. 

Most  of  the  physical  quantities  measured  in  this 
experiment  involve  higher  dimensions  and  multiple 
variables.  In  order  to  be  able  to  investigate  a  wide 
parameter  space  the  experiment  has  been  inter¬ 
faced  with  high  speed  computers.  This  approach 
opens  a  new  dimension  in  experimental  work. 
Without  digital  data  processing  it  would  be  essen¬ 
tially  impossible  to  investigate  the  seven-variable 
distribution  function  or  the  Fourier  transformed 
cross-correlation  tensor  of  magnetic  fluctuations. 
We  will  display  some  of  these  functions  in  stereo¬ 
scopic  form;  the  details  for  viewing  such  figures  are 
described  in  the  appendix. 

2.  Experimental  arrangement 

The  present  reconnection  experiment  (see  fig.  1) 
employs  a  linear  discharge  plasma  column  pro¬ 
duced  with  a  large  diameter  (1  meter)  cathode  [19]. 
Typical  densities  are  10'^ cm”’,  temperature 


kT^»  10 eV,  axial  magnetic  fields  By»  10 G.  The 
plasma  is  uniform,  quiescent,  nearly  collisionless, 
and  highly  reproducible  in  pulses  of  duration 
/p  »  5  ms  repeated  every  r, «  2  s.  After  establishing 
the  plasma  a  time  varying  transverse  magnetic  field 
»  lOG)  is  applied.  Its  topology  in  vacuum 
contains  an  X-type  neutral  point  on  the  axis  of  the 
device  (see  fig.  la).  In  the  presence  of  the  plasma 
electron  currents  (/p  a:  1000  A)  are  induced  which 
flow  preferentially  in  regions  of  BiiaO  but  are 
sufficiently  lafge  to  modify  the  applied  field  topol¬ 
ogy.  During  the  period  of  rising  magnetic  fields 
(0<r  <  100/is)  the  self-consistent  field  distribu¬ 
tion  is  that  of  a  classical  neutral  sheet  [S].  The 
plasma  current,  driven  by  inductive  and  space 
charge  electric  fields,  E  =  —A  —  Vil>p,  is  mainly 
carried  by  electrons  which  drift  from  the  cathode 
toward  the  end  anode.  The  current  closes  via  the 
grounded  chamber  wall  to  the  cathode.  No  net 
currents  flow  to  the  insulated  plates  or  to  the  radial 
chamber  walls  which  are  separated  from  the 
plasma  by  vacuum. 

The  plasma  processes  are  analyzed  with  various 
probes  which  provide  time  and  space  resolved 
data.  A  digital  data  acquisition  system  consisting 
of  8  bit,  100  MHz  analog-to-digital  (A/D)  con¬ 
verters,  minicomputer  (LSI  11/23),  array  processor 
and  direct  data  link  to  large  size  computers  (CDC 
7600,  Cray  1)  is  employed.  Probes  include  magnetic 
loops,  Langmuir  probes,  high-frequency  rf  probes, 
and  a  novel  directional  velocity  analyzer  [20].  In  a 
typical  magnetic  field  measurement  the  signals  of 
orthogonal  loops  are  recorded  vs.  time  at  300 
spatial  locations.  At  each  position  the  experiment  is 
repeated  up  to  80  times  so  as  to  obtain  statistical 
averages,  i.e.,  mean,  standard  deviation,  and  cor¬ 
relation  functions.  Extensive  space-time  maps  may 
involve  >  100  million  numbers  which  can  only  be 
handled  by  high  speed  computers. 

3.  Magnetic  fieM  topologies 

A  function  of  multiple  variables  such  as  the 
vector  field  B(x,  y,  z,  i)  cannot  be  visualized  at 
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Fig.  I.  Schematic  pictuK  of  the  experimental  arnmgement.  (a)  Cross-sectional  view  showing  parallel  plate  electrodes  with  pulsed 
currents  /,  and  magnetic  field  lines  without  plasma,  (b)  Side  view  of  the  device  with  main  electrodes,  currents  (fp,  /J  electric  fields 
(E  =  -  A  —  F^p),  and  magnetic  fields  (B  =  B^  +  B^).  The  coordinate  system  common  in  magnetospheric  physics  has  been  adopted 
where  y  is  along  the  neutral  line  (device  axis),  x  is  along  the  horizontal  neutral  sheet,  and  2  is  normal  to  the  sheet. 


once  but  only  in  a  reduced  parameter  space.  We 
start  to  show  the  two-dimensional  field 
in  the  transverse  plane  (x,  z)  at  a 
given  time  t,  ensemble  averaged  at  each  point  over 
N  =  25  repeated  experiments.  Fig.  2a  depicts  the 
field  topology  of  a  neutral  sheet,  the  classical 
configuration  treated  in  most  steady-state  theories 
of  reconnection  [21],  We  observe  this  field  pattern 
during  the  rise  of  the  applied  currents.  The  induced 
plasma  current  tends  to  slow  down  the  penetration 
of  the  external  field  and  its  reconnection  in  the  null 
region.  The  current  sheet  7  =  P  x  has  a  half 
width  of  dz  %  3  cm  which  is  on  the  order  of  the 
collisionless  skin  depth  c/cupe « 0.7  cm.  Current 
sheets  investigated  in  space  are  often  broader  than 
an  ion  Larmor  radius  so  that,  in  comparison,  we 
are  resolving  the  6ne  structure  of  current  sheets. 
The  field  topology  is  macroscopically  stable  for 
many  Alfven  transit  times  across  the  plasma 
«  7  )is  in  He,  20  ns  in  Ar).  However,  in  the  case 
of  very  large  current  densities  (/  %  nev,)  disruptive 
instabilities  are  observed  [22].  These  will  not  be 
considered  in  this  paper. 

In  response  to  the  induced  currents  and  J  x  B 
forces,  the  plasma  is  vertically  compressed  and 
horizontally  expelled.  We  have  measured  the  fluid 
motion  with  orthogonal  differential  ion  velocity 
analyzers  [12].  Fig.  2b  depicts  the  characteristic  ion 
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Fig.  2.  Measured  magnetic  field  and  flow  field  during  recon¬ 
nection.  (a)  Transverse  field  B^(x,  z)  at  a  fixed  time  t,  position 
y  =  137  cm.  showing  the  classical  neutral  sheet  topology,  (b) 
Transverse  ion  flow  velocity,  Vj^(x,  z),  normalized  to  the  local 
sound  or  Alfven  speed,  Cj  = 

Vertical  compresrion  and  horizontal  jetting  are  in  qu^tative 
agreement  with  reconnection  models. 

flow  pattern  which  is  observed  to  develop  after  an 
initial  phase  of  turbulence.  The  measured  flow  field 
clearly  shows  the  fluid  jetting  at  Alfvenic  speeds, 
which  is  predicted  by  theory  and  computer  simu¬ 
lations  [23]  and  thought  to  explain  the  generation 
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of  high  energy  particles  in  solar  flares  and  the 
magnetotail. 

We  now  turn  to  the  time  dependence  of  the 
magnetic  field  topology.  Two  time  scales  can  be 
distinguished:  a  slow  one  corresponding  to  the  rise 
time  of  the  external  fields  (t  ^  100 /is),  and  a  fast 
time  scale  associated  with  fluctuations  and  in¬ 
stabilities.  The  slowly  varying  processes  are  well 
reproducible  over  repeated  experiments.  Thus,  en¬ 
semble  averages  such  as  shown  in  fig.  2  are  gener¬ 
ated  in  typically  1024  successive,  small  time  in¬ 
crements.  The  time  evolution  of  the  field  topologies 
is  then  viewed  in  sequence  in  the  form  of  a 
computer  processed  movie.  Here,  we  can  display 
only  a  few  frames  depicting  the  major  features. 

Fig.  3  shows  magnetic  field  lines  in  the  transverse 
x-z  plane  at  three  different  time  steps,  t  =  57,  58, 
and  59  fts.  Instead  of  the  vector  field  of  fig.  2a,  we 
have  integrated  the  B  field  and  display  contours  of 
constant  vector  potential,  =  const.,  where 
B  =  F  X  A.  The  motion  of  magnetic  field  lines 
through  the  plasma  is  best  followed  by  the  two 
inner  field  lines  dravm  heavier.  With  increasing 
time  the  opposing  field  lines  move  vertically  ( ±  2) 
towards  the  center,  merge  at  the  null  point,  and 
reconnect  to  form  a  new  field  line  pair  which  moves 
horizontally  ( ±  x)  outward.  This  process  of  field 
line  reconnection  occurs  continuously  during  the 
external  current  rise.  Magnetic  flux  generated  at 
the  two  current  carrying  plates  is  transferred  across 
the  separatrix  into  regions  of  common  flux  where 
field  lines  encircle  both  current  plates.  During  the 
current  decrease  the  process  reverses  but  the  re¬ 
sultant  vertical  current  sheet  usually  tears  to  form 
a  large  magnetic  island  [10].  The  motion  of  field 
lines  through  plasmas  is  highly  descriptive  but 
becomes  physically  significant  only  through  the 
coupling  of  plasma  and  field  lines.  However,  near 
the  null  region,  diffusion  due  to  finite  resistivity 
dominates.  The  field-line  “velocity”  is  then  a  mea¬ 
sure  of  the  axial  inductive  electric  field  which  is 
often  taken  as  another  measure  for  the  recon¬ 
nection  rate  [21].  The  effective  plasma  resistivity  is 
an  important  parameter  since  it  determines  the  size 
of  the  diffusion  region,  the  reconnection  rate,  and 
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Fig.  3.  Snapshots  of  magnetic  field  lines  (X,  =  const.)  during 
reconnection  in  a  plasma.  The  two  inner  field  lines  merge  and 
reconnect  to  form  a  pair  of  horizontally  outward  moving  field 
lines.  The  entire  time  evolution  has  been  displayed  in  a  motion 
picture. 

the  plasma  heating.  Direct  measurements  of  the 
resistivity  [13]  show  a  level  well  above  the  classical 
one,  and  an  inhomogeneous  pattern  in  the  neutral 
sheet.  These  features  may  be  understood  after 
having  discussed  wave  and  particle  measurements. 

The  present  discussion  was  so  far  limited  to  the 
two-dimensional  field  topology.  Because  of  the 
axial  component  the  total  magnetic  field  is 
three-dimensional,  B  =  B^,  B^.  Fig.  4  is  a 

display  of  the  three-dimensional  vector  field  in  the 
transverse  plane  It  is  presented  as  a  picture 
pair  for  stereoscopic  viewing,  the  techniques  for 
which  are  explained  in  the  appendix.  The  field  lines 
are  straight  only  in  the  current  sheet,  and  exhibit 
increasing  shear  both  above  and  below  it.  The 
neutral  sheet  is  a  null  region  only  with  respect  to 
the  transverse  field  (B,  =  fl,  =  0)  but  not  the  total 
field  (|B|  =  Byo^  10 G).  Thus,  the  electrons  in  the 
current  sheet  always  remain  magnetized.  In  a 
different  investigation  we  have  varied  the  axial  field 
strength  from  B^  =  0  to  By^^  |Bx|  but  here  we 
limit  ourselves  to  the  case  « |Si|.  Both  in 


32 


I 


R.L.  Stenzel  and  If.  Gekelman  I  Magnetic  field-line  reconnection  in  plasmas 


137 


Fig.  4.  Stereoscopic  display  of  the  three-dimensional  vector  magnetic  field  B  =  in  the  transverse x -  z plane.  For  reference, 

bottom  plane  gives  a  projection  of  the  transverse  field  lines  /4j,(jii;z)  =  const.  (y=  137  cm).  See  appendix  for  viewing  stereo  pictures. 


tokamaks  and  in  the  magnetosphere  the  fields  are, 
in  general,  three-dimensional,  while  most  the¬ 
oretical  models  assume  By  =  0,  for  simplicity. 
When  the  axial  field  component  is  absent  the 
magnetic  field  gradients  are  largest.  Fig.  5  presents 
a  stereoscopic  view  of  the  field  strength 
along  the  neutral  sheet  which  forms  a  deep  long 
valley.  In  such  inhomogeneous  fields  the  particle 
orbits  are,  in  general,  nonadiabatic  and  the  distri¬ 
butions  are  anisotropic  [24]. 

We  now  turn  to  the  rapidly  varying  magnetic 
field  fluctuations  {dBjB  «  10%).  These  result  in 
rippled  and  distorted  field  lines.  Although  with  a 
single  probe  we  cannot  map  the  instantaneous  field 


lines,  we  can  study  the  nature  of  the  fluctuations 
statistically,  e.g.,  by  two-point  cross-correlations. 
The  first  step  in  this  investigation  is  to  determine 
the  fluctuation  spectrum  in  terms  of  characteristic 
modes.  Once  the  waves  arc  identified  the  instability 
mechanism  and  the  saturation  processes  have  to  be 
identified,  which  requires  the  knowledge  of  particle 
distributions.  The  final  step  is  to  assess  the  effect 
of  the  instability  on  the  macroscopic  plasma  pa¬ 
rameters,  i.c.,  the  transport  properties. 

Cross  correlation  studies  are  performed  with 
two  identical  pairs  of  magnetic  probes,  one  mov¬ 
able  ( #  1)  and  one  stationary  ( #  2).  At  a  given 
position  an  ensemble  of  fV  =  80  repeated  mea- 


Fig.  S.  Stcfcoscopic  view  of  the  transverse  magnetic  field  strength  z)|  showing  steep  gradients  adjacent  to  neutral  sheet 
(0<  iB  <  15G,  -25  <  jc  <25  cm,  -12<z  <  12  cm,  t  =  IS  pi).  See  appendix  for  viewing  stereo  pictures. 
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surements  of  the  vector  components  B„  B,  is 
recorded  vs.  time.  Due  to  the  slow  time  variation 
the  fluctuation  processes  are  weakly  nonstationary 
so  that  statistical  averages  must  be  obtained  by 
ensemble  averaging  rather  than  time  averaging 
[2S].  From  the  stored  data  we  calculate  on-line  the 
ensemble  mean 

the  root  mean  square  value  of  the  fluctuations, 
/I  " 

.  (2) 

and  the  cross  correlation  function 

(3) 

W  j.i 


We  note  that  the  correlation  of  two  vectors  results 
in  a  tensor  of,  at  least,  rank  R  =  2.  Here,  its 


components  are  given  by 


{BM\ 

{b,,b^)) 


(4) 


In  order  to  identify  the  mode  spectrum  we  Fourier 
transform  the  time  dependence  of  the  fluctuations, 
select  a  narrow  spectral  component,  and  analyze  its 
spatial  properties.  The  resultant  cross  spectral 
function  is  a  function  of  probe  separation 
Ar  =  (Ax,  Ay,  Ag),  delay  time  t  between  the  two 
probes,  and  a  weaker  function  of  the  time  t  during 
the  reconnection  process,  <J?iB2>„(dr,  t,  t). 

In  fig.  6  we  display  an  example  of  a  diagonal 
tensor  component  (.B^iB^y  in  the  axial  x-y  plane 
at  a  given  frequency  (co/2n  =  1  MHz),  time 
(I  =  21.4  fis),  and  delay  time  (t  =  1.3  ns).  The  up¬ 
per  part  of  the  figure  shows  a  topological  view  of 
the  cross  spectral  function  while  the  lower  part 
presents  a  corresponding  contour  plot.  The  dis¬ 
plays  on  the  right-and  left-hand  side  are  rotated  by 
ss  6"  for  purpose  of  stereoscopic  viewing  (see  ap¬ 
pendix  A). 

The  measurement  indicates  that  even  at  a 
selected  frequency  a  broad  spectrum  of  spatial 
modes  is  excited.  This  behavior  is  also  found  at 


Fig.  6.  Stereoscopic  view  of  the  cross-spectral  tensor  components  of  magnetic  field  fluctuations.  The  data  are  taken 

at  t  =  2l.4^s  during  reconnection  with  r  =>  1.3  /is  time  delay  between  the  two  probe  signals,  each  filtered  at  a>l2n  ==  I  MHz.  A 
two-dimensional  spatial  mode  spectrum  is  observed.  See  appendix  for  viewing  stereo  pictures. 
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Other  frequencies  in  the  range  investigated  persion  relation  for  whistlers  [26] 

(100  kHz  </  <  1  MHz).  The  characteristic  electro¬ 
magnetic  modes  in  this  frequency  band  are  the  /kc 
whistler  waves  [26].  These  modes  are  not  only 
dispersive  (to/k  ^  const.)  but  also  anisotropic,  i.e., 
at  a  given  to  the  wave  number  k  depends  on  the  where  tojln  =  1  MHz,  electron  plasma  frequency 

direction  of  wave  propagation  with  respect  to  the  mp/2B  =  20  GHz,  and  cyclotron  frequency 

magnetic  field  [k  =  k(0)].  From  our  cross-spectral  tojln  =  28  MHz  are  typical  experimental  parame- 

measurements  we  can  obtain  the  wavenumber  ters.  One  can  see  that  the  unstable  modes  fall  close 

spectrum  by  a  spatial  Fourier  transformation  from  so  the  theoretical  refractive  index  {kcjto)  surface 

r-space  into  A-space.  Since  the  time  and  space  for  whistlers.  An  exact  fit  cannot  be  expected  due 

dependence  is  recorded  the  complex  fast  Fourier  to  magnetic  field  and  density  inhomogeneities, 

transform  yields  both  magnitude  and  direction  of  Besides  the  dispersion  properties  we  have  also 

k.  Using  the  data  from  two  orthogonal  planes  investigated  the  polarization  of  the  unstable 

(x-y,  x-z)  we  have  constructed  the  three-  modes.  Theory  predicts  that  the  magnetic  field  of 

dimensional  wave  vector  spectrum  from  the  two  whistler  waves  has  a  right  hand  circular  polar- 

projections  in  the  k^-k^  and  k^-k^  planes.  ization  [26].  By  applying  a  digital  filter  to  the 

An  example  of  such  data  analysis  is  shown  in  fig.  measured  k  vector  spectrum  we  select  one  mode  at 

7.  It  is  a  stereoscopic  display  of  the  dominant  a  time,  inverse  Fourier  transform  it,  and  inspect  its 

modes  of  in  the  full  three-dimensional  k  space-time  dependence  (dr,  t).  This  procedure  is 

space.  The  relative  intensity  of  the  modes  is  pro-  applied  to  the  two  tensor  components 

portional  to  the  dot  size.  A  more  detailed  ampli-  and  which  have  two  orthogonal  com- 

tude  profile  is  provided  by  contour  plots  in  the  two  ponents  of  vector  B]  correlated  with  the  same 

orthogonal  planes.  In  order  to  compare  the  mea-  reference  signal  A  display  of  <B,Brt>  versus  r 

surements  with  theory  we  also  display  the  wave  at  a  given  Jr  (or  vice  versa)  provides  a  hodogram 

vector  surfaces  ±A,«)  calculated  from  the  dis-  which  reveals  the  polarization. 


Three  Dimensional  k  Spectra  Three  Dimensional  k  Spectra 


Fig.  7.  Three-dimenMonal  wavevector  space  (*„  *,)  with  measured  modes  of  the  cross-spectral  fuDCtiou  and  the  theoretical 

dispersion  surfaces  ±ik,„  for  whistler  waves  {to /In  =  1  MHz,  tojln  =  28  MHz,  <oJ2n  =  20  GHz).  Vertical  axis  k,  is  along  average 
origin  k  =  0  is  indicated  by  square  dot.  The  mode  intensity  is  proportional  to  the  dot  size,  the  mode  location  k  is  constructed  from 
measurements  of  <S,|S,2>  in  two  orthogonal  spatial  planes,  Fourier  transformed  into  k  planes  and  presented  as  contour  plots  adjacent 
to  the  k  space.  The  observations  indicate  that  the  magnetic  turbulence  consisu  of  oblique  whistler  modes.  See  appendix  for  viewing 
stereo  pictures. 


co(a)  —  tOc  cos  9)  ’ 
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Fig.  8  is  a  stereoscopic  view  of  the  hodogram 
vs.  T  for  a  mode  propagating  nearly  along 
the  average  magnetic  field.  The  helical  trajectory 
indicates  that  the  mode  is  essentially  right-hand 
circularly  polarized.  However,  for  long  delay  times 
(t  »  3  /is)  phase  mixing  due  to  the  finite  filter 
bandwidth  (dm/co  %  10~')  tends  to  randomize  the 
polarizations.  For  modes  propagating  at  an  angle 
to  the  magnetic  field  we  first  transform  the  wave 
A-field  into  a  plane  normal  to  A  (A  *  A  =  0)  and 
then  inspect  the  polarization.  Most  modes  which 
fall  close  to  the  dispersion  surface  are  right-hand 
polarized. 

Thus,  from  both  the  dispersion  and  polarization 
properties  the  eigenmodes  which  make  up  the 
magnetic  field  fluctuations  are  found  to  be  whistler 
waves. 

We  will  next  investigate  the  cause  for  the  in¬ 
stability  of  whistler  waves.  According  to  theo¬ 
retical  models,  anisotropies  in  the  electron  velocity 
distribution  can  lead  to  whistler  instabilities  [27]. 


4.  Dfstribution  fanctioas 

The  distribution  of  particles  in  velocity  space  v, 
real  space  r,  and  time  I,  f(v„r,  t),  provides  the 
fundamental  description  of  plasmas  on  a  mic¬ 


roscopic  level.  All  macroscopic  plasma  properties 
can  be  derived  from  it  by  forming  moments,  i.e., 
averages  over  the  velocity  distribution.  Kinetic 
theories  and  numerical  particle  simulations  all 
depend  on  the  knowledge  of  the  distribution  func¬ 
tion.  In  spite  of  the  theoretical  significance  of  this 
function  there  arc,  surprisingly,  few  experimental 
observations  of,  e.g.,  electron  distributions  in  labo¬ 
ratory  plasmas.  In  order  to  measure /(r)  a  velocity 
analyzer  with  directional  sensitivity  is  required. 
Such  detectors  are  available  for  low  density,  high 
energy  particles  in  space  [28]  but  have  only  recently 
been  developed  for  denser  laboratory  plasmas 
[20,29]. 

Our  electron  detector  consists  of  a  small  (3  mm 
radius)  retarding  potential  velocity  analyzer  whose 
directional  sensitivity  (dQ%4nxlO~’  sterad) 
arises  from  filtering  electrons  through  a  micro- 
channel  plate.  It  can  be  rotated  through  the  two 
spherical  angles  (0<0<7t,  0<i^< 2n)  so  as  to 
collect  electrons  from  any  direction  in  velocity 
space.  The  three-dimensional  distribution  function 
f{Vy9,4>)  is  assembled  from  individual  flux  mea¬ 
surements  at  typically  200  different  directions  9,  ifi. 
Furthermore,  at  each  angular  position  the  experi¬ 
ment  is  repeated  N  times  so  as  to  form  the  ensem¬ 
ble  average  </>  and  standard  deviation  anal¬ 
ogous  to  the  case  of  the  magnetic  field 


Fig.  8.  Polarization  properties  of  a  single  (m.k)  mode  of  the  magnetic  turbulence  spectrun.  With  increasing  delay  time  r  the 
correlation  vector  <.81,812)  rotates  initially  in  a  right-hand  circular  sense  as  expected  for  whistler  waves  {co/2n=  I  MHz).  See 
appendix  for  viewing  stereo  pictures. 
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Fig.  12.  Comparison  between  the  average  distribution  function  and  its  fluctuations  in  the  neutral  sheet. 

Upper  panel:  Logarithmic  display  of  the  average  value  <^v||,  Vi)>  and  two  one-dimensional  cuts 
along  v||  =  v,.  axis  and  Vi  =  VjXxis.  Note  tail  of  runaway  electrons.  [/f0.0)ai  10  '^cm  ‘  s’;  reference 
plane  at /=  10  cm  <’  s’.] 

Middle  level:  Contours  of  constant  average  log</(v||,  Vi)>  (3  dB  contour  spacing). 

Bonom  level:  Contours  of  root-mean-square  fluctuations  log/nn^  >1|.  fi)  =  const.  (3  d  B  per  contour, 
/mu  mix/ <fimx>®  20%).  Note  large  fluctuations  in  regions  of  runaway  electron  tails.  [0<|  vj,  |  v  |  < 
6xl0*cm/s.] 
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measurements  [see  eqs.  (I),  (2)].  Besides  the  veloc¬ 
ity  space  dependence  we  also  measure  the  real 
space  dependence  by  moving  the  analyzer  to  vari¬ 
ous  positions  in  the  neutral  sheet.  The  mea¬ 
surements  are  also  time  resolved  to  within 
dt%3/is,  the  time  required  to  sweep  out  the 
current-voltage  characteristics.  The  data  pro¬ 
cessing  for  a  function  of  up  to  seven  variables  is 
made  possible  only  through  the  use  of  high  speed 
computers  with  large  memory. 

Fig.  9  presents  measurements  of  the  electron 
distribution  function  in  velocity  space.  In  two 
dimensions  the  distribution  / (u„  «„)  can  either  be 
displayed  topographically  (fig.  9a)  or  as  a  contour 
plot  (fig.  9b).  In  three-dimensional  velocity  space  a 
display  analogous  to  fig.  9a  cannot  be  visualized 
but  we  can  extend  fig.  9b  and  plot  surfaces  of 
constant  /(r)  as  done  in  fig.  9c.  The  maximum 
value I0~'^cm‘‘s’  generates  a  surface  near 
the  origin.  Decreasing  values  of  f(v)  produce 
nested  surfaces  of  expanding  radius  in  velocity 
space,  three  of  which  are  shown  in  fig.  9c  where  the 
outer  surfaces  are  cut  open  for  purpose  of  display. 
Using  time  as  the  variable  /we  have  displayed  the 
continuous  variation  0  </  <  in  a  computer 
generated  movie. 

Anisotropic  distribution  functions  exhibit  non- 
spherical  surfaces  of/(e)  =  const.  Fig.  10  shows  a 
stereoscopic  view  of  a  surface  near  the  tail  of  the 
distribution.  We  choose  to  display  v\f  (»)  =  const, 
rather  than  /(r)  =  const,  to  stress  the  importance 
of  a  seemingly  small  tail  (see  fig.  9)  on  the  moments 
of/(i»).  The  quantity  displayed  is  the  kernel  for  the 
zeroth  moment,  i.e.,  the  density 


2k  ft  00 


”e  =  j/(»)dr  =  0)dn  sin  0  dO  dtf) . 


0  0  0 


(6) 


Higher  moments  such  as  those  leading  to  the 
current  density  y  =  jr/(e)de,  mean  energy 
£  =  j\mvY(v)dv,  and  heat  flow  ^  = 
J \nw^vf  (r)  dr  involve  progressively  higher  powers 
of  V  which  puts  a  strong  weight  on  the  fast  tail 
electrons.  For  example,  integration  over  the  distri- 


Fig.  9.  Measured  electron  distribution  function  in  the  neutral 
sheet,  (a)  Linear  display  of/(B|,Oi);  10-'‘cm-‘s’.  (b) 

Contour  plot/(i)(,  i'^)  =  const.;  =  10%  per  contour,  (c) 

Distribution  function  in  three-dimensional  velocity  space  dis¬ 
played  as  surfaces  of/(i)„  Vy,  v,)  =  const.  Outer  surfaces  are  cut 
for  purpose  of  display.  Continuous  variation  of  0  </  <f^  has 
been  displayed  in  a  motion  picture. 
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Fig.  10.  Stereoscopic  display  of  a  surface  uV(u)  =  const  (50  cm"*  s)  in  three-dimensional  velocity  space.  The  distortion  on  the  left- 
hand  side  reflects  an  anisotropy  caused  by  runaway  electrons  in  the  electric  field  Ey\Byo  of  the  current  sheet.  t)j-axis  along  maximum 
protrusion  to  the  left  (180°  rotation  with  respect  to  fig.  9c).  See  appendix  for  viewing  stereo  pictures. 

button  function  of  fig.  9  reveals  that  the  asym-  nection.  The  tail  of  high  energy  (50-100  eV)  elec- 

metric  tail  electrons  have  a  relative  density  trons  is  thought  to  arise  from  the  free  acceleration 

njn^ »  10%,  energy  EJE » 22%,  and  are  re-  of  essentially  collisionless  tail  electrons  by  the 

sponsible  for  the  current  and  heat  flow  in  the  electric  field  Ey  along  By^.  Outside  the  neutral  sheet 

plasma  which  is  at  a  small  angle  (A6  »  15°)  with  the  transverse  components  B^  #  0)  prevent  the 

respect  to  the  local  magnetic  field.  free  acceleration.  This  is  clearly  visible  in  fig.  11b 

The  distribution  of  figs.  9  and  10  are  observed  in  which  displays  distribution  functions,  log/ (Vy,  pj, 

the  middle  of  the  neutral  sheet  during  recon-  at  two  different  positions  (x,z)  pointed  out  by 


Fig.  II.  Spatial  variation  of  the  distribution  function  in  the  neutral  sheet,  (a)  Magnetic  field  topology  with  locations  where/(v)  was 
measured,  (b)  Electron  distribution / (t>,,  v,)  on  a  logarithmic  scale  with  reference  plane  three  orders  of  magnitude  below  the  maximum. 
Note  the  presence  of  an  energetic  electron  tail  in  the  middle  of  the  neutral  sheet  and  absence  of  runaway  electrons  at  the  sheet  edge. 
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arrows  in  the  simultaneously  measured  magnetic 
field  map  BJ^x,z)  of  fig.  11a.  It  is  interesting  to 
note  that  the  bulk  electron  temperature  is  higher  at 
the  side  of  the  neutral  sheet  than  in  the  middle 
where  the  energy  is  transferred  to  runaway  elec¬ 
trons. 

The  presence  of  runaway  electron  tails  in  a 
current  carrying  plasma  can  give  rise  to  various 
kinetic  instabilities.  For  example,  the  emission  of 
plasma  waves  and  lower  hybrid  waves  in  low 
density  tokamak  discharges  is  well  known  [30]. 
Auroral  kilometric  radiation  is  associated  with 
energetic  electron  precipitation  [31].  Oblique  whis¬ 
tler  waves  can  be  excited  by  parallel  beams  [32]. 
Here,  we  are  particularly  interested  in  low  fre¬ 
quency  electromagnetic  instabilities  (0.1  < 
/  <  1  MHz)  in  order  to  establish  the  relation  with 
the  observed  unstable  whistler  modes.  Whistler 
instabilities  arise  from  wave-particle  interactions 
by  cyclotron  resonance  [(cu  —  nco^lky  «  Oj]  and 
Cerenkov  resonance  [aalky  =  Uj]  [33].  We  have, 
therefore,  investigated  low  frequency  fluctuations 
in  the  distribution  function  in  the  velocity  regime 
satisfying  wave-particle  resonance. 

Fig.  1 2  (color  plate)  displays  in  the  bottom  plane 
[(i’ll  t’l)  =  (Pyy  *’r)]  the  rms  fluctuations  of  the  elec¬ 
tron  distribution,  given  by 

f  1  IV  -iMi 

=  1^1  [/(«’) -</(«’)>]^|  (7) 

The  average  distribution,  </(p)>  = 

is  shown  for  comparison  as  a 
contour  plot  in  the  middle  plane  and  topo¬ 
graphically  in  the  upper  plane.  The  distributions 
are  displayed  on  a  logarithmic  scale  so  as  to  obtain 
a  wide  dynamic  range  in  /(r)  (1000:1,  3dB  per 
contour,  /«,„«,/</ (0)>  «  0.2).  The  observations 
show  that  near  the  runaway  electron  tail  of  <  f(p)> 
the  fluctuation  level  is  enhanced  as  expected  from 
kinetic  instability  theories.  The  characteristic  ve¬ 
locity  of  the  tail  particles  is  i;|«3.5x  10*  cm/s 
close  to  the  fundamental  (n  =  1)  cyclotron  reso¬ 
nance  {coJky»4x  10* cm/s).  Neither  the  second 
resonance  («  =  2)  nor  the  Cerenkov  resonance 


(/j=0.  <o//t|  »  1.5  X  lO^cm/s)  involves  a  sig¬ 
nificant  number  of  resonant  particles.  Thus,  it 
appears  that  the  tail  electrons  excite  whistlers  by 
cyclotron  resonance.  This  plausible  model  awaits  a 
final  proof,  i.e.,  the  cross  correlation  between 
particle  and  wave  fluctuations, 

5.  Discussion  and  conclusion 

In  the  present  experiment  the  classical  recon¬ 
nection  problem  has  been  investigated.  Because  of 
a  well-controlled  laboratory  plasma,  refined 
plasma  diagnostics  and  the  power  of  digital  data 
processing  we  could  extract  information  previously 
not  available.  We  found  that,  to  first  order,  the 
magnetic  field  topology  and  fluid  flow  are  estab¬ 
lished  in  accord  with  theory.  However,  mea¬ 
surements  of  the  resistivity  and  fluid  acceleration 
showed  orders  of  magnitude  discrepancies  with 
classical  fluid  parameters  [12, 13].  A  study  of  the 
microscopic  processes  such  as  fluctuations  and 
distributions  has  been  undertaken.  The  electron 
distribution  in  the  current  sheet  is  found  to  be 
anisotropic  due  to  a  tail  of  runaway  electrons.  This 
is  a  very  important  result  since  it  shows  that  MHD 
theories  break  down  in  the  neutral  sheet.  For 
example,  the  resistivity  now  is  determined  by  the 
collisionality  of  the  runaway  electrons  rather  than 
the  bulk  electron  temperature.  The  resistivity  may 
become  spatially  very  inhomogeneous  due  to 
changes  in  the  distribution  function. 

Runaway  electrons  are  a  source  of  various  mi¬ 
croinstabilities  one  of  which  has  been  analyzed  in 
great  detail.  The  magnetic  fluctuations  have  been 
identified  from  their  dispersion  and  polarization  to 
consist  of  a  broad  spectrum  of  whistler  modes.  The 
growth  of  the  wave  spectrum  modifies,  in  turn,  the 
electron  distribution.  The  tail  electrons  lose  some 
of  their  energy  and  momentum  to  the  waves  which 
results  in  an  effective  resistivity  for  the  plasma 
current.  The  background  electrons  acquire  energy 
through  wave  absorption  which  saturates  the  in¬ 
stabilities  [34].  Ions  are  coupled  to  electrons  by 
space  charge  electric  fields. 
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Thus,  the  energy  transfer  during  reconnection  in 
collisionless  plasmas  involves  various  nonlinear 
wave-particle  processes  which  cannot  be  analyzed 
in  the  framework  of  fluid  theory.  Refined  experi¬ 
ments  are  necessary  to  clarify  these  “anomalous” 
transport  processes. 
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Appendix  A 

In  order  to  view  figs.  4-8  and  10  stereo- 
scopically,  a  simple  viewing  device  is  helpful.  It 
consists  of  two  identical  lenses  of  focal  length 
/»  15  cm  and  diameter  d  »  2.5  cm,  mounted  ap¬ 
prox.  7  cm  apart  and  held  about  10-15  cm  sym¬ 
metrically  above  the  figure.  A  divider  between  the 
two  figures  and  lenses  may  be  helpful  in  separating 
the  two  images. 

A  commerically  made  inexpensive  viewer  has 
been  purchased  by  the  authors  from  the  Taylor 
Merchant  Corporation,  212  W.  35th  Street,  New 
York,  N.Y.  10001  (Stereopticon  Model  707). 
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Magnetic  Field  Line  Reconnection  Experiments 
6.  Magnetic  Turbulence 

W.  Gekelman  and  R.  L.  Stenzel 

Department  of  Physics,  Unioersity  of  Cal^omia,  Los  Angeles 

Extensive  statistical  analysis  of  the  vector  components  of  fluctuating  magnetic  fields  have  been  per¬ 
formed  in  a  time  dependent  neutral  magnetic  sheet.  Cross  spectra]  analysis  indicates  a  variety  of  wave 
numbers  present  for  each  frequency  investigated.  Comparison  of  Fourier  components  of  the  cross  spec¬ 
tral  function  with  dispersion  surfaces  in  k  space  demonstrates  the  waves  are  large  amplitude  whistlers; 
this  is  verified  by  pojarization  analysis,  which  shows  the  random  waves'  magnetic  fields  to  be  right-hand 
circular.  Ion  acoustic  and  Langmuir  turbulence  are  also  observed  along  with  bursts  of  microwave 
radiation.  Measurements  of  the  electron  distribution  function  /(v,  r.  ()  and  its  fluctuations  in  velocity 
space  relate  wave  and  particle  activity. 


Introduction 

Fluctuations  are  a  pervasive  property  of  any  active  medium. 
Noise  that  may  be  considered  to  be  a  superposition  of  many 
waves  has  been  long  observed  in  both  laboratory  and  space 
plasmas.  When  free  energy  is  available,  as  in  a  non- 
Maxwellian  particle  distribution,  wave  amplitudes  can  grow 
and  turbulence  ensues. 

Magnetic  fluctuations  have  been  observed  by  satellites  in 
the  plasma  sheet  and  its  boundary  layer  [McFerron,  1979] 
and  in  the  magnetopause  [Sonnerup,  1976;  Russell  and  Elphic, 
1978;  Gurnett  et  al.,  1979a;  Anderson  et  al.,  1982].  “Lion 
roars”  or  intense  bursts  of  electromagnetic  radiation  have 
been  explored  in  the  magnetosheath  [Tsurutani  et  al.,  1982; 
Sudan  and  Denavit,  1973]  and  have  been  identified  as  whistler 
waves  [Smith  and  Tsurutani,  1976].  In  addition,  whistlers  have 
been  observed  in  the  bow  shock  [Coroniti  et  ai,  1982].  Satel¬ 
lites  have  detected  fluctuating  magnetic  fields  in  the  mag¬ 
netosphere  [Russell  et  al.,  1972].  The  Voyager  1  spacecraft 
data  established  the  presence  of  whistler  waves  in  the  Jovian 
magnetosphere  [Scarf  cmd  Gurnett,  1977].  These  waves  may  be 
caused  by  a  powerful  perturbation  such  as  a  lightning  stroke 
[Gurnett  et  al.,  19796]  or  by  an  anisotropy  in  the  electron 
distribution  function  [Senimtin  et  al.,  1983].  The  underlying 
physics  of  whistler  waves  has  been  treated  in  the  classical 
works  of  Storey  [1953]  and  Helliwell  [1%5]. 

The  electron  distribution  function  in  a  current  carrying 
plasma  may  often  develop  into  non-Maxwellian;  such  a 
(riasma  is  subject  to  a  variety  of  instabilities  [Kadomtsev, 
1965].  The  ion  acoustic  instability  and  ion  sound  turbulence 
has  been  extensively  investigated  in  the  laboratory  [Gekelman 
and  Stenzel,  1978,  and  references  therein]  and  observed  in  the 
solar  wind  [Gurnett  and  Frank,  1978]. 

In  this  paper  we  wish  to  report  observations  of  these  insta¬ 
bilities  as  well  as  Langmuir  turbulence  in  a  magnetic  field  line 
reconnection  experiment  containing  a  neutral  sheet.  Although 
these  waves  are  not  unique  to  the  present  magnetic  field  topol¬ 
ogy,  they  are  a  natural  consequence  of  the  currents  along  the 
neutral  sheet  and  the  anisotropic  distribution  functions 
therdn.  These,  in  turn,  are  in  part  due  to  the  sheared  magnetic 
fields.  Once  these  waves  become  unstable  their  feedback  onto 
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the  reconnection  event  is  quite  important.  Ion  sound  turbu¬ 
lence  can  cause  turbulent  electron  heating  [Oum,  1978]  and 
dramatically  raise  the  resistivity  above  the  classical  Spitzer 
value  in  the  reconnection  region  [Wild  et  al.,  1981].  In  these 
experiments  the  Hall  term  (V  x  B)  in  Ohm’s  law  was  mea¬ 
sured  and  found  to  be  negligible  with  respect  to  other  terms. 
An  estimate  of  the  inertial  conductivity  places  it  orders  of 
magnitude  larger  than  that  due  to  Coulomb  oscillations.  The 
resistivity  governs  the  reconnection  rate  [Fosyliunos,  1975; 
Sonnerup,  1979];  in  fact,  without  finite  resistivity,  the  field 
lines,  in  accord  with  Maxwell’s  equations,  would  pile  up  at  the 
separatrix  and  never  reconnect  at  all.  To  simulate  reconnec¬ 
tion  on  computers,  ad  hoc  resistivities  are  introduced  [Sato 
and  Hasegawa,  1982]  and  are  required  for  theoretical  tearing 
mode  calculations  [Furth  et  ai,  1963]. 

The  fluctuating  magnetic  fields  associated  with  the  turbu¬ 
lent  whistlers  affect  the  electron  orbits.  Although  the  ions  near 
the  neutral  sheet  are  essentially  unmagnetized,  they  are  cou¬ 
pled  to  the  electrons  through  space  charge  fields  such  that  the 
ion  collision  rate  may  become  anomalously  high  [Stenzel  and 
Gekelman,  1981a]. 

In  a  controlled  laboratory  experiment  it  is  possible  to  study 
the  structure  of  the  magnetic  and  electrostatic  turbulence  in 
great  detail.  The  variety  of  things  happening  simultaneously  in 
a  turbulent  plasma  implies  that  an  overwhelming  amount  of 
data  must  be  taken  to  sort  them  out.  The  introduction  of 
relatively  inexpensive,  fast  computers  has  led  to  a  quantum 
jump  in  what  is  possible  experimentally.  Fully  three- 
dimensional  measurements  of  the  electron  velocity  distri¬ 
bution  have  been  performed  [Stenzel  et  ai,  1983a]  and  will  be 
discussed  here  and  in  a  forthcoming  paper  in  this  journal.  We 
point  out  that  spacecr^'A  measuring  wave  activity  detect  the 
waves  at  a  Doppler  shifted  frequency  owing  to  a  combination 
of  their  motion  and  possible  fluid  drifts  relative  to  the  earth. 
Single  probes  take  data  at  one  spatial  location  and  cannot 
determine  the  wave  number  spectra  from  correlation  analyses 
such  as  those  reported  in  this  work. 

The  present  reconnection  experiment  has  already  shown  de¬ 
tailed  diagnostics  of  the  space  time  evolution  of  the  magnetic 
fields  [Stenzel  and  Gekelman,  19816],  the  plasma  parameters 
[Gekelman  and  Stenzel,  1981],  the  fluid  dynamics  [Gekelman  et 
al.,  1982a],  the  energy  flow  [Stenzel  et  al.,  1982],  and  current 
disruptions  with  double  layer  formations  [Stenzel  et  al., 
19836].  These  will  be  referred  to  throughout  the  text  as  JGR 
parts  1-5.  An  overall  review  of  these  experiments  is  given  in  a 
recent  issue  of  Physica  Scripta  [Gekelman  et  al.,  19826]. 
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Fig.  I.  Schematic  view  of  the  experimental  device,  (a)  Cross-sectional  view  with  transverse  ma^etic  field  line  pattern 
in  vacuum.  (6)  Side  view  with  characteristic  fields  and  currents.  Note  that  the  axis  of  the  device  is  in  the  >  direction  while 
the  transverse  coordinates  are  x,  z,  as  customary  in  magnetospheric  physics. 


Experimental  Arrangement 

The  reconnection  experiments  are  performed  in  a  cylindrical 
vacuum  chamber  (l.S-m  diameter,  2-m  length)  in  which  a  low 
pressure  {p  10"*  torr,  A,  Hj,  He)  discharge  is  produced  with 
a  1-m  diameter  oxide-coated  cathode  (Figure  1).  The  emission 
properties  of  such  special  cathodes  [StenzW  and  Daley,  1980] 
provide  a  high-quality  laboratory  plasma  of  the  following 
properties:  Uniformity  over  large  soUe  lengths  (An/n  e:  10% 
over  80  cm  across  Bo,  180  cm  along  Bg),  nearly  collisionless, 
high  beta  parameters  (density  n, «  10'^  cm'^,  temperature 
kT,  Si  10  eV,  electron-ion  mean  free  path  si  200  cm,  axial 
magnetic  field  =:  12  G,  ^  =  nkTlB^/lpo)  si  2),  high  repeti¬ 
tion  rate  (t,,,  =  2  8),  and  good  reproducibility  {Sn/n  S  5%)  of 
the  pulsed  discharge  (F^,,  s  40  V,  /^j,  a;  1500  A,  t^j,  a;  5  ms). 
The  discharge  pulse  length  is  more  than  sufficient  to  reach 
steady-state  plasma  conditions,  at  which  time  the  reconnec¬ 
tion  experiment  is  initiated.  This  involves  establishing  a  trans¬ 
verse  magnetic  field  B^^  =  (B,,  B,)  with  a  null  point.  Note  that 
we  have  adopted  x,  y,  i  coordinates  commonly  used  for  the 
magnetotail  where  the  y  direction  is  along  the  separator,  x  is 
along  the  neutral  sheet,  and  z  is  normal  to  it.  A  large  sinu¬ 
soidal  current  pulse  (/,  2;  20  kA,  rise  time  i,  s:  100  /is)  is 
drawn  axially  through  two  parallel  aluminum  plates  (75  cm 
wide,  33  cm  spacing,  180  cm  length)  located  on  top  and 
bottom  of  the  plasma  column.  The  plates  are  insulated  from 
the  plasma,  which  is  confined  to  the  space  between  them. 
Axial  and  average  transverse  magnetic  fields  are  of  com¬ 
parable  strength. 

The  transverse  magnetic  field  geometry  in  vacuum  is  shown 
schematically  in  the  cross  sectional  view  of  the  device  (Figure 
Id).  Without  plasma  there  is  an  X  type  neutral  point  on  the 
axis  of  the  device.  However,  with  plasma,  axial  plasma  cur¬ 
rents  are  induced  (/,  1000  A)  that  modify  the  magnetic  field 

topology.  The  self-consistent  current  distribution  and  field  to¬ 
pology  can  result  in  the  formation  of  a  classical  neutral  sheet 
[Dungey,  1958],  The  plasma  current  driven  by  inductive  and 
space  charge  electric  fields  (E  =  -A  -  Td^,)  is  mainly  carried 
by  the  electrons.  During  the  initial  phase  of  external  current 
rise  [dljdt  >0,  0  <  t  100  /<s)  the  electrons  drift  from  the 
cathode  toward  the  end  anode.  The  current  closes  via  the 
grounded  chamber  walls  to  the  cathode.  No  net  current  flows 
to  the  insulated  plates  or  to  the  radial  chamber  walls  which 
are  well  separated  from  the  plasma  column. 

The  cathode  emits  as  many  electrons  as  are  collected  at  the 
end  plate  so  as  to  maintain  approximate  space  charge  neu¬ 
trality  prior  to  the  reconnection  event  However,  as  the  ap¬ 
plied  inductive  electric  field  (A  £  0.5  V/cm)  can  exceed  the 


classical  runaway  field  [Dreicer,  1959]  (£,  =:  i  0.1 

V/cm,  where  =  IkTJm,)  and  since  the  space  charge  limit¬ 
ed  emission  [Langmuir,  1929]  does  not  provide  for  current 
densities  j  nev,  the  plasma  sets  up  an  axial  space  charge 
electric  field  —  opposite  to  the  induced  field.  This  shield¬ 
ing  effect  limits  the  electron  drift  to  <  v„  at  least  on  the 
average  over  the  entire  plasma  cross  section.  However,  by 
decreasing  the  plasma  density,  the  space  charge  contribution 
may  be  substantially  lessened. 

The  evaluation  of  statistical  quantities  such  as  rms  field 
fluctuations  and  correlation  tensors  is  made  possible  through 
a  data  acquisition  system  and  processor.  The  wave  forms  are 
digitized  by  fast  (Transiac  100  MHz,  LeCroy  20  MHz)  analog 
to  digital  (A/D)  converters  that  are  controlled  by  an  LSI  1 1/23 
computer.  The  computer  also  drives  stepping  motors  that 
position  probes  within  the  device.  The  computer  is  too  slow  to 
perform  the  necessary  arithmetic  such  as  fast  Fourier  trans¬ 
forms,  digital  filtering,  and  correlations  during  the  1.5  s  be¬ 
tween  shots.  An  array  processor  (CSPl  Map  200),  therefore, 
has  been  interfaced  to  the  computer.  Upon  performing  all  the 
calculations,  the  MAP  places  the  results  in  the  core  of  the  LSI 
II;  it  is  subsequently  placed  on  disc.  Upon  completion  of  a 
data  run  the  final  data  set  is  shipped  via  a  direct  link  to  large 
(Cray,  CDC  7600)  offsite  computers. 


Previous  Experimental  Results 

Before  addressing  ourselves  to  fluctuations  it  is  helpful  to 
review  the  average  magnetic  field  topology  (part  I)  as  well  as 
other  macroscopic  properties.  Plate  1  is  a  stereo  slide  pair  that 
should  be  examined  with  a  proper  viewer.  (Note:  Plates  1-3, 
Plates  5  and  6,  and  Figure  15  are  best  viewed  in  three  dimen¬ 
sions  because  of  the  complexity  of  the  data.  This  may  be 
accomplished  with  a  viewer  available  through  the  Taylor  Mer¬ 
chant  Corporation,  212  W.  35  Street,  New  York,  New  York 
10001  (Stereopticon  707).  One  may  be  constructed  out  of  two 
14-cm  focal  length,  2.S-cm  diameter  lenses  and  cardboard.  The 
lens  centers  should  be  placed  7  cm  apart  and  13.5  cm  above 
the  page,  with  a  cardboard  divider  midway  between  them. 
Color  reprints  of  this  article  and  viewers  are  available  for 
S4.50  U.S.,  in  care  of  the  authors  at  Department  of  Physics, 
University  of  California,  Los  Angeles,  405  Hilgard  Avenue, 
Los  Angeles,  California  90024.) 

The  vector  potential  has  an  axial  component  A,  =  J  B,  dx 
—  f  B,dz,  d/dy  s:  0,  which  is  calculated  from  the  transverse 
components  of  the  magnetic  field  and  displayed  on  the  lower 
x-z  plane.  Small  (1  cm  in  diameter)  orthogonal  magnetic  loops 
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Fig.  2.  (a)  Schematic  view  of  recording  the  magnetic  field  versus 
time  in  form  of  a  hodogram.  (h)  Single  shot  hodogram  taken  near  the 
neutral  sheet  identifies  the  fluctuations  as  loops,  i.e.,  a  rotating  field 
component  associated  with  hydromagnetic  waves,  (c)  Superposition  of 
three  hodograms  taken  at  same  position  as  (b)  shows  dominance  of 
magnetic  fluctuations,  (d)  Superposition  of  four  hodograms  when 
probe  is  removed  from  neutral  sheet  region.  He  3  x  10'*  torr,  = 
20G,y  =  0. 


detect  the  changing  B  field,  which  is  recorded  at  each  of  300 
spatial  locations  (x,  z)  as  a  function  of  time  (1024  time  steps. 
At  =  100  ns)  and  averaged  over  an  ensemble  of  25  shots.  One 
sees  that  a  classical  neutral  sheet  topology,  associated  with  an 
axial  current  sheet,  d,  =  (V  x  Bj^)/tto>  been  formed. 

Superimposed  upon  the  transverse  magnetic  field  is  a  con¬ 
stant  axial  bias  field  (0^  =  12  G)  such  that  the  total  magnetic 
field  is  sheared.  This  is  illustrated  as  a  three-dimensional 
vector  field  B(x,  y,  z)  =  B^  -I-  Bj^  in  Plate  1.  For  clarity,  only  a 
fraction  of  the  vectors  are  plotted.  In  the  center  of  the  neutm! 
sheet  the  vectors  are  vertical  since  the  transverse  component 
of  B  vanishes,  and  at  the  boundaries  (z  =  ±12  cm)  the  vectors 
are  inclined  at  45°  since 

From  rapidly  swept  planar  Langmuir  probe  traces  analyzed 
with  our  digital  data  acquisition  system,  the  ensemble 
averaged  density  and  electron  temperatures  are  mapped  as  a 
function  of  space  and  time  (part  2).  The  electron  temperature 
increases  strongly  but  nonuniformly  during  the  reconnection 
process,  ranging  from  T,  5  eV  in  the  center  of  the  neutral 
sheet  to  T,  2;  15  eV  in  argon  (30  eV  in  helium)  near  the  hori¬ 
zontal  edges  of  the  neutral  sheet.  A  density  minimum  develops 
in  the  center  of  the  neutral  sheet  (n,  <:  lO’^/cm’)  and  maxima 
(n,  a:  2  X  lO'^/cm’)  near  horizontal  edges  of  the  neutral  sheet. 
By  combining  density  and  temperature  to  form  the  plasma 
pressure  we  found  that  the  forces  on  the  fluid  are  not  in  bal¬ 
ance  (J  X  B  —  Vp  #  0).  This  leads  to  an  acceleration  and  jet¬ 
ting  of  fluid  from  the  neutral  sheet.  The  resulting  flow  fields 
were  measured  with  differential  ion  velocity  analyzers  (part  3) 
and  found  to  be  highly  turbulent  in  space  and  time  during  the 
first  30  /iS  of  the  experiment  when  the  induced  electric  field  is 
large.  When  the  turbulence  level  has  decreased  (the  neutral 
sheet  topology  is  present  throughout  this  time)  the  classical 
J  X  B  flow  establishes  itself  with  inflow  velocities  c,  i  0.2 
and  outflow  velocities  d,  ^  0.5  i\,  where  is  the  local  Alfven 
speed  based  on  the  total  magnetic  field.  From  the  measured 
fluid  flow,  pressure,  fields,  and  currents  we  have  quantitatively 
obtained  the  phenomenological  scattering  term  Af  =  mnv*v  in 
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Fig  3.  Histograms,  which  compare  the  number  of  occurrences  of  a  fluctuation  as  a  function  of  its  amplitude,  are 
displayed  for  the  magnetic  field  and  white  noise  using  an  ensemble  of  2000  events.  Both  variables  may  be  fit  to  the 
superimposed  Gaussian  distributions  with  99%  confidence. 
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the  fluid  equation 


CORRELATION  MEASUREMENTS 


mn  —  =  (J  X  B  -  Vp)  -  Af  (1) 

at 

The  components  of  the  collision  tensor  v*  were  evaluated 
and  found  to  be  large  within  the  neutral  sheet  during  times  of 
turbulent  flow  and  subsequently  vanish. 

Magnetic  Field  Fluctuations 

Although  the  initial  plasma  parameters  and  applied  cur¬ 
rents  are  highly  reproducible  from  shot  to  shot,  we  observe 
magnetic  field  fluctuations  inside  the  plasma.  Figure  2  shows  a 
magnetic  field  hodogram,  i.e.,  the  trajectory  of  the  vector  tip 
versus  time  during  a  180-/i$  time  interval.  Just  as  in  a 
satellite  measurement,  the  hodogram  is  taken  at  a  single  spa¬ 
tial  position  that,  in  Figures  2b  and  2c,  is  within  the  neutral 
sheet  where  the  noise  is  largest.  The  fluctuations  manifest 
themselves,  in  part,  as  small  loops  in  the  vector  plot  (Figure 
2b).  These  can  be  interpreted  as  a  rapidly  rotating  magnetic 
field  component  associated  with  a  wave.  The  random  nature 
of  these  waves  is  evident  in  Figure  2c,  which  has  hodograms 
from  three  repeated  events  superimposed.  At  this  location  the 
fluctuations  are  comparable  to  the  average  value  of  If  the 
probe  is  removed  from  the  neutral  sheet  (Figure  2d)  the  irre- 
producible  segment  of  the  hodogram  collapses.  From  Figures 
2b  and  2c  one  can  estimate  the  fluctuating  magnetic  field  to  be 
of  order  1-2  G  and,  therefore,  =  10-20%. 

The  spatial  properties  of  the  fluctuations  within  one  shot 
cannot  be  explored  with  a  single  probe.  Correlation  functions 
that  require  taking  a  large  number  of  shots  at  many  positions 
necessitate  the  use  of  a  digital  data  acquisition  system.  For 
example,  we  have  recorded  2S  shots  of  B(r)  at  each  of  300 
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Fig.  4.  The  autocorelation  function  for  digitally  filtered  (Q  =  6) 
white  noise  as  a  function  of  S,  the  number  of  shots  in  an  ensemble. 
As  A  ->  00,  the  373-kHz  correlation  signal  should  decay  to  1/e  in  six 
periods.  An  80-mefflber  ensemble  chosen  experimentally  always 
guaranteed  convergence. 
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Fig.  S.  Schematic  of  probe  arrangement  for  correlation  measure¬ 
ments.  A  fixed  reference  probe  is  located  in  the  x,  z  plane  at  y  ==  0  cm. 
A  second  probe  is  mounted  to  move  in  the  x-z  plane  (x  >  x„,^, 
y  =  0,  - 12  cm  s  z  s  12  cm)  or  in  the  x,  y  plane  (x  >  x,|„j,  -26 
cm  <  y  <  +26  cm,  z  =  0).  For  reference  a  linear  vector  plot  of  is 
placed  in  the  x.  z,  y  =  0  cm  plane.  He  3  x  10"*  torr,  =  20  G, 
I  =  27  /IS.  The  current  carrying  sheets  are  located  at  z  =  ±16  cm.  The 
cathode  is  located  at  y  =  —  1 37  cm. 

positions  and  calculated  on-line  both  the  ensemble  mean 
values  of  the  field  magnitude  (.B^y,  direction  <d)  and  the  root 
mean  square  values  B^^  and  The  rms  fluctuation  or 
autocorrelation  that  contains  all  frequencies  within  the  probe 
and  amplifier  bandwidth  (1  kHzg/^2  MHz)  maximizes 
within  the  neutral  sheet,  which  is  also  the  location  of  large 
current  density.  These  are  of  order  =:  1.  The  spatial 

extent  of  the  magnetic  waves  is  frequency  dependent;  for  I 
MHz  the  disturbance  fills  the  entire  transverse  plane. 

Nonstationary  Turbulence 

Since  the  magnetic  turbulence  affects  the  particle  motion 
and  ultimately  the  transport  properties,  detailed  studies  of  it 
have  been  undertaken.  The  analysis  involves  measurements  of 
cross  correlation  functions.  If  the  time  averaged  mean  value  of 
a  fluctuating  magnetic  field  <B(r)>  is  independent  of  the  choice 
of  the  averaging  interval,  the  random  variable  B(I)  is  station¬ 
ary.  In  this  case  the  correlation  function  may  be  calculated 
from  a  sufficiently  long  temporal  record  of  each  held  compo¬ 
nent.  If  a  single  damped  mode  is  a  part  of  the  random  vari¬ 
able,  there  exists  a  technique  developed  by  Matthaeus  and 
Goldstein  [1982]  that  subtracts  this  mode  and  treats  the  re¬ 
mainder  as  stationary  noise,  A  process  is  called  weakly  non¬ 
stationary  when  the  mean  values  of  B(x,  y,  z)  change  slowly  in 
comparison  with  the  fluctuation  frequencies  as  the  neutral 
sheet  is  formed.  For  nonstationary  processes  the  statistical 
properties  such  as  correlation  functions  must  be  calculated 
from  ensembles  of  random  events  that  we  establish  by  repeat¬ 
ing  the  experiment  under  identical  initial  conditions. 

Each  component  of  the  fluctuating  magnetic  field  is  a 
random  variable  whose  probability  distribution  should  ap¬ 
proach  a  Gaussian  with  increasing  number  of  events  [Loeve, 
I960].  This  basic  property  has  been  verified  by  generating  a 
histogram  of  taken  by  a  probe  in  the  center  of  the  neutral 
sheet,  shown  in  Figure  3<i,  for  an  ensemble  of  2(X)0  events. 
Superimposed  on  the  data  is  a  Gaussian  distribution  with  the 
same  mean  and  standard  deviation.  A  chi  square  test  [Taylor, 
1982]  indicates  the  distribution  of  to  be  normal  with  a 
99%  probability.  Shown  in  the  adjacent  Figure  3b  is  the  result 
using  white  noise  as  the  input  signal,  which  is  essentially  iden¬ 
tical  to  that  of 

In  evaluating  the  correlation  functions,  one  must  determine 
how  many  events  constitute  an  adequate  ensemble.  The  result 
determines  the  data  acquisition  time  and  required  computer 
memory. 
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Fig.  6.  Contours  of  a  component  of  CSF  tensor 

(f  =  I  MHz)  at  four  delay  times  r  in  fis  with  respect  to  the  measure¬ 
ment  t  =  33.2  pts  into  the  reconnection  event.  The  peaks  oscillate  at 
the  wave  frequency  but  do  not  propagate,  which  is  characteristic  of  a 
standing  wave.  The  average  wavelength  is  of  order  10  cm.  The  lo¬ 
cations  of  the  peaks  referred  to  in  the  text  A  =  y4(x,  y)  =  A(—22.2, 
4.3),  B  =  (0.0,  3.3),  C  =  (2.5,  -  5.5). 


This  was  approached  heuristically  by  a  convergence  test. 
Output  from  a  white  noise  generator  (20  kHz  </ <  I  MHz) 
was  fed  into  the  20-MHz  analog  to  digital  converters  and 
digitally  filtered  [Liu,  197S]  in  the  same  way  as  the  magnetic 
fluctuations.  The  quality  factor  Q  =  ml  duo  of  a  filter  deter¬ 
mines  the  number  of  periods  it  rings  during  an  e  folding  decay 
time  in  response  to  a  pulsed  input.  The  normalized  auto¬ 
correlation  function  is  given  by 


^  Z  C'»(fo)Cit(fo  + 


where  C((o  -I-  t)  is  a  signal  with  no  mean  values,  measured  at 


delay  time  r  with  respect  to  the  start  of  measurement  and  N  is 
the  number  of  records  of  the  ensemble.  For  delay  times 
T  >  Q(2n/m}  and  the  ensemble  numbers  large  enough  for 
phase  mixing  to  occur  the  correlation  decays,  R^x)-*  0.  This  is 
illustrated  in  Figure  4,  where  the  autocorrelation  is  displayed 
for  three  ensemble  sizes.  The  input  signal  is  white  noise  digi¬ 
tally  filtered  at  373  kHz  with  Q  =  6.  RfO)  is  one  as  expected 
but  only  for  N  ^  SO  R(t)  decays  properly.  It  was  therefore 
decided  to  use  an  80  shot  ensemble  for  the  cross  spectral 
analysis  since  this  represented  an  upper  limit  on  the  computer 
memory  for  storage  of  four  components  and  seemed  to  con¬ 
verge  well. 

Cross  Spectral  Function  Analysis 

Assuming  the  noise  to  be  composed  of  a  large  number  of 
randomly  generated  waves,  a  correlation  analysis  is  performed 
in  order  to  identify  modes  from  their  dispersion  relation  cu(k) 
and  polarization  properties. 

Two  identical  probe  sets  are  used  to  measure  the  vector 
components  of  B^fr,  t)  at  two  different  positions  r.  During 
each  reconnection  event  the  signals  are  digitized  by  four 
20-MHz  analog  to  digital  converters,  then  fast  Fourier  trans¬ 
formed  (FFT)  and  digitally  filtered  by  the  on-line  array  pro¬ 
cessor  (MAP200).  An  inverse  FFT  is  done  to  generate  the  time 
history  of  the  filtered  signals  in  the  frequency  band  of  interest. 
Eighty  consecutive  shots  are  stored  in  the  MAP  memory,  and 
the  following  operations  are  performed  on  the  signals  filtered 
at  frequency  m  with  a  filter  Q  =  10 

-  <Bi.Jt)>  (3) 

where  i  =  x,  z  denotes  the  vector  component,  8,"“  is  the  total 
field  measured  during  the  nth  shot,  and  <Bj>  =(\fN) 

8,,'°*  is  the  average  field  at  a  time  t.  As  discussed  in  the 
previous  section,  an  ensemble  number  N  =  80  is  sufficiently 
large  for  a  convergence  of  the  mean  values. 

The  cross-spectral  function  for  nonstationary  fluctuations 
[BenJat  and  Piersol,  1971]  calculated  on-line  is  given  by 
I  so 

t,r)  =  —  Y.  Bi,.,jT,  t)Bj2.,,Jr  +  &.T,  t  +  r) 
w  1 

(4) 

Here  I  and  2  refer  to  the  two  probe  sets.  The  correlation  is 
large  if  in  moving  from  the  fixed  probe  2  over  distance  Ar  a 
signal  arrives  after  time  t  in  (or  180°  out  of)  phase  at  probe  I. 
If  there  is  a  decay  of  the  signal  strength,  a  rotation  of  the  field 
vector,  or  a  loss  of  coherence  with  Ar  and  t,  the  correlation 
will  die  out. 

The  cross-spectral  function  (CSF)  may  be  normalized  in  the 
form 


Pij(Ar,  t,  T,  a>)  = 


_ <B.iBj2>(Ar,  t,  T,  to) _ 

<B„,„"(r,  t)y'\Bj2,Ar  +  Ar,  t  -)-  t)>‘« 


(5) 


so  that  it  lies  between  p  =  ±  1.  Since  the  two  terms  in  the 
denominator  are  stored  separately,  it  is  possible  to  evaluate 
both  the  normalized  and  the  unnormalized  correlation  func¬ 
tions. 

The  normalized  CSF  magnifies  the  low  level  waves  since  the 
denominator  of  (5)  is  proportional  to  B,^  Bz™.  which  can  be 
misleading  in  estimating  the  spatial  properties. 

We  point  out  that  the  filtered  signal  from  a  pickup  loop  is 
proportional  to  oiB,  where  m  is  the  center  frequency  of  the 
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digital  filter.  The  signals,  therefore,  were  not  integrated  since 
the  common  factor  a  drops  out  in  (5)  and  is  an  overall  multi¬ 
plicative  constant  for  the  unnormalized  CSF  (equation  (4)). 

As  observed  with  a  single  pickup  loop,  the  fluctuations  in 
B,  and  B,  are  of  comparable  amplitude.  The  individual  probes 
are  connected  to  isolation  transformers  and  diflerential  ampli¬ 
fiers  that  exhibited  slightly  diflerent  responses  at  higher  fre¬ 
quencies  such  as  1  MHz.  The  magnitudes  and  delay  time  t  of 
the  CSF  were  therefore  corrected  by  a  calibration  magnetic 
field  set  up  by  a  test  wave  to  be  described  below. 

With  two  vector  components  for  each  probe  the  CSF  be¬ 
comes  a  tensor  of  rank  2,  given  by 

p  =  (6) 

The  off-diagonal  elements  are  large  if  over  the  displacement 
(Ar,  t)  an  x  component  of  the  field  rotates  into  the  z  direction 
or  vice  versa.  The  trace  describes  the  correlation  of  the  field 
strengths  and  converges  for  Ar  =  r  =  0  to  Tr  <B|B2>  =  B^^. 
The  total  number  of  on-line  calculations  to  get  p„  at  300 
spatial  positions  in  a  plane,  64  delay  times  r,  eight  starting 
times  t  of  intervals  AT  =:  3.2  /is  is  staggering  (>2  x  10*). 
Nevertheless,  we  have  measured  p„  at  several  frequencies  and 
in  two  orthogonal  planes. 

The  experimental  layout  for  this  measurement  is  shown  in 
Figure  S.  A  stationary  reference  probe  is  positioned  at  z  -  0  in 
a  transverse  plane  located  137  cm  from  the  cathode.  Shown  to 
scale  is  a  measured  vector  plot  of  the  transverse  magnetic 
field.  A  second  probe  is  free  to  move  everywhere  in  either  the 
same  transverse  x-z  plane  or  in  the  orthogonal  x-y  plane  con¬ 
taining  the  neutral  sheet.  Each  probe  shaft  has  two  loops 
oriented  to  measure  B,  and  B,. 

Plate  2  shows  the  cross-spectral  intensity  <B,,B,2>  at /  =  I 
MHz  in  the  horizontal  x-y  plane.  The  measurements  were 
insensitive  to  the  absolute  position  of  the  fixed  reference  probe 
2  within  the  neutral  sheet.  For  all  the  data  described  here  the 
fixed  probe  is  located  at  (x,  y,  z)  =  ( —  20  cm,  0,  0)  to  maximize 
the  area  scanned.  Note  that  the  axis  of  the  device  goes 
through  the  center  of  the  plane  of  Figure  5.  The  cathode  in 
this  coordinate  system  is  at  y=— 137  cm.  The  three- 
dimensional  display  in  the  upper  part  and  its  two-dimensional 
projection  as  a  contour  map  are  taken  at  t  =  36.4  ps  after  the 
start  of  the  reconnection  pulse  with  signal  B.j  delayed  by 
I  =  1.3  /IS  with  respect  to  B,,.  The  observations,  Plate  2,  show 
a  wavelike  pattern  in  the  x-y  plane.  Variations  in  the  delay 


cj  Spectrum 


0  0.5  1 

cj/2n  fMHz) 


Fig.  7.  Frequency  spectra  during  a  single  reconnection  event 
obtained  from  an  on-line  fast  Fourier  transform.  The  magnetic  noise 
extends  to  the  limiting  frequency  of  the  detector  s  2  MHz  and  is  a 
factor  of  10  above  instrument^  noise  at  1  MHz.  The  lower  hybrid 
hequency  is  200  kHz,/,,  3  1  KHz, /„  =  28  MHz 


j  c)  <BzlBx2>  Spectra 


Fig.  8.  Wave  number  spectra  obtained  from  a  two-dimensional 
spatial  Fourier  transform  of  the  cross  spectral  function  as  shown  in 
Plate  2.  In  general,  the  data  indicate  a  spread  in  k  space  even  at  a 
single  frequency  component  of  the  magnetic  fluctuation  spectrum  and 
are  multimode  in  both  the  x-y  and  x-z  planes.  Two  tensor  compo¬ 
nents  and  <B,,B,,>  part  (C)  have  a  dominant  mode  in  at 

least  one  plane.  This  feature  makes  possible  an  unambiguous  polar¬ 
ization  analysis  of  the  mode.  The  spectra  are  normalized  such  that  0 
S  CSF(k)  s  I  with  10%  increments. 

time  indicate  no  preferred  direction  of  wave  propagation.  The 
other  components  of  the  tensor  have  similar  spatial  properties. 
The  fluctuations  in  B,  and  B,  are  of  equal  amplitude  when 
observed  with  a  single  pickup  loop. 

In  the  simplest  case  of  a  single  wave  number  and  frequency, 
the  expected  pattern  would  be  a  succession  of  wave  fronts 
with  a  maximum  amplitude  at  the  reference  probe  location 
[Gekelman  and  Stenzel,  1978].  If  the  wave  propagation  were 
spatially  isotropic,  the  pattern  would  appear  like  the  ripples  in 
a  pond  after  a  pebble  was  dropped  into  it.  Since  the  present 
case  is  more  complicated,  we  begin  to  make  further  inferences 
about  the  waves  by  analyzing  the  CSF  contours  as  a  function 
of  delay  time  t. 
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TEST  WAVE  k- SPECTRA 
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Fig.  9.  Wave  number  spectra  of  a  test  wave  in  the  plane  contain¬ 
ing  the  axial  magnetic  held,  =  10  G,  and  neutral  sheet  (z  =  0), 
generated  from  a  two-dimensional  spatial  Fourier  transform  of  the 
CSF  displayed  in  Plate  4<i.  A  dominant  broadband  mode  lc„  =  0.3S 
cm“‘,  <c,  =  —0.29  cm''  propagates  at  fl  =  tan“‘  =  -40‘.  The 
spectrum  in  the  transverse,  x-z  plane  is  similar  but  centered  at  the 
origin,  =  =  0,  indicating  no  preferred  direction  of  propagation. 


Figure  6  i$  a  contour  map  of  the  CSF  tensor  component 
1  versus  r.  The  start  time  is  33.2  ps  into  the 
reconnection  event.  To  aid  the  eye,  three  peaks  labeled  A,  B, 
and  C  have  been  singled  out  and  are  seen  to  diminish  in 
amplitude  in  the  100  ns  between  Figures  6a  and  6b.  Roughly  a 
quarter  period  after  the  sequence  start  (Figure  6c),  peak  A  has 
entirely  vanished;  B  and  C  are  nearly  gone.  At  approximately 
a  half  cycle  from  the  outset  (Figure  6d)  the  peaks  have  re¬ 
versed  sign  but  none  of  them  have  moved  spatially.  The  CSF 
pattern  has  the  characteristics  of  a  standing  wave,  the  inter¬ 
ference  of  oppositely  propagating  modes  at  this  frequency. 
There  is  no  preferred  direction;  on  the  average,  waves  travel 
toward  and  away  from  the  fixed  probe  with  equal  probability. 
Measurement  of  the  separation  between  adjacent  minima  and 
maxima  permits  an  estimate  of  the  average  wavelength, 
</l>  =:  10  cm.  Frequency  and  wavelength  are  typical  for  whis¬ 
tler  waves  in  the  present  parameter  regime. 

Since  the  digital  filter  in  this  experiment  has  a  value  Q  =  10, 
a  significant  amplitude  decay  of  the  CSF  is  not  seen  during 
the  3-ps  delay  time  with  respect  to  each  start  time  t.  The 
average  amplitude  of  the  CSF  does,  however,  diminish  as  the 
reconnection  event  progresses.  All  components  are  a  factor  of 
2  smaller  at  f  =  50  /<$  since  the  induced  electric  field  and 
neutral  sheet  current  also  decay  in  time. 

Data  taken  in  the  x-z  plane  have  a  similar  structure  to 
those  of  the  x-y  plane.  In  principle,  a  three-dimensional  corre¬ 
lation  measurement  in  a  volume  of  space  would  be  required  to 
determine  the  true  wave  number.  Because  of  the  nearly  iso¬ 
tropic  wave  propagation  two  orthogonal  planes  are  adequate. 

At  lower  frequencies  the  wavelengths  approach  the  system 
boundary  size.  The  set  of  data  at  o)/2n  =  500  kHz  had  a  simi¬ 
lar  structure  to  that  in  Figure  6;  however,  the  wavelengths 
were  longer  and  for  cu/2it  =  200  kHz,  which  is  just  above  the 
average  lower  hybrid  frequency,  the  CSF  is  largest  in  the  neu¬ 
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tral  sheet  region  and  appears  to  be  a  standing  wave.  To  avoid 
any  complications  due  to  bounded  modes,  most  of  the  analy¬ 
sis  was  focused  on  the  highest  frequency  (1  MHz)  mode. 

Wave  Number  Spectra 

The  analysis  of  the  spatial  structure  of  the  CSF  does  not 
lead  to  a  positive  identification  of  the  modes.  Our  study  pro¬ 
ceeds  with  a  two-dimensional  Fourier  transformation  of  the 
spatial  distribution  so  as  to  obtain  a  wave  number  spectrum 
(1^  t)  at  a  given  frequency  to.  This  information  is 
complementary  to  the  frequency  spectrum,  a  typical  example 
of  which  is  shown  for  6,  in  Figure  7. 

At  each  sample  time  (  the  CSF  is  calculated  for  64  delay 
times  T  (At  =  50  ns).  There  is  sufficient  temporal  resolution  so 
as  to  calculate  terms  of  the  form  =  cos 

(k  •  r  —  tor)  -F  i  sin  (k  •  r  —  tor),  from  data  points  one  quarter 
period  out  of  phase.  Using  a  complex  kernel  in  the  two- 
dimensional  FFT  enables  the  wave  number  spectrum  to  be 
calculated  in  a  plane  without  symmetry  assumptions.  The 
two-dimensional  wave  number  spectra  of  the  tensor  compo¬ 
nents  and  are  displayed  in  Figure  8.  The 

data  in  Figures  8u  and  8c  refer  to  the  x-y  plane,  which  con¬ 
tains  the  axial  guide  field  {B^  =  12  G)  and  runs  along  the 
center  of  the  neutral  sheet  (z  =  0).  This  k  spectrum  is  acquired 
at  r  =  33.2  /is  into  the  reconnection  event  when  the  turbulence 
is  fully  developed.  The  center  axis  of  the  grid  corresponds  to 
4,  =  k,  =  0.  Since  it  is  not  possible  to  measure  infinite  wave¬ 
lengths,  the  values  on  these  axes  correspond  to  averages  taken 
from  those  above  and  below.  We  point  out  that  the  intervals 
along  the  abscissa  and  ordinate  are  not  equal  because  the  data 
are  taken  on  a  rectangular  grid  (3.14  cm  £  \k,\  ^  24  cm;  6.85 
cm  ^  |a,|  ^  50  cm  in  this  display).  The  probes  can  resolve 
structures  of  order  1  cm.  The  wave  number  spectrum  obtained 
from  the  x-y  plane  exhibits  a  dominant  mode,  as  shown  in 
Figure  8c  for 

The  magnetic  turbulence  exhibits  a  broad  range  of  a  values 
even  at  a  single  frequency  <a.  In  the  frequency  regime  between 
the  lower  hybrid  and  the  electron  cyclotron  frequency  (col„  < 
<<>«  cUj,  where  zx  0.1  MHz,  to,/2)t  =;  28  MHz)  oblique 

whistler  modes  have  the  properties  of  exhibiting  a  broad  range 
of  wave  numbers  depending  upon  the  angle  of  wave  propaga¬ 
tion  with  respect  to  the  magnetic  field;  they  are  right-hand 
circularly  polarized  electromagnetic  modes  {Helliwell,  1965]. 

The  data  are  still  too  complex  for  comparison  with  theory, 
so  that  a  further  mode  selection  is  required.  The  cross  spectral 
functions  are  filtered  in  k  space.  The  dominant  modes  of  the 
CSF  components  are  selected  and  individually  examined  after 
an  inverse  Fourier  transformation  into  coordinate  space.  Plate 
3a  displays  in  real  space  a  single  mode  =  48  cm,  2,-5  cm) 
of  <B,iB,2>  a  multi  .node  spectrum  not  unlike  those  of 
Figure  8.  The  function  is  sinusoidal  with  phase  fronts  nearly 
parallel  to  B^.  As  will  be  discussed  below,  the  phase  fronts 
may  be  easily  tracked  as  a  function  of  t  and  the  propagation 
studied.  If  a  second  mode  is  included  (Plate  36)  the  inter¬ 
ference  becomes  obvious,  although  the  contour  plots  projected 
below  are  still  symmetric.  For  four  modes  (Plate  3c)  the  pat¬ 
tern  begins  to  approach  that  of  Plate  2,  and  tracking  the 
motion  of  peaks  as  well  as  comparison  with  other  tensor  com¬ 
ponents  of  the  cross-spectral  function  becomes  as  arduous  as 
what  was  previously  described.  By  transforming  away  all 
modes  save  one,  spatial  randomness  is  lost  although  the  prop¬ 
agation  direction,  wavelength,  and  frequency  are  known  with 
certainty.  We  are  now  in  a  position  to  analyze  the  wave  polar¬ 
ization  by  comparing  amplitudes  and  phase  of  two  orthogonal 
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Plate  1 .  A  stereo  pair  that  illustrates  the  sheared  magnetic  field  B(x,  y.  z)  =  B^(x.  z)  +  The  projection  of  the 
three-dimensional  field  lines  A,  is  drawn  in  the  x-z  plane.  The  transverse  magnetic  neutral  sheet  within  which  B^  s  0  is 
clearly  visible.  He  3  x  10  *torr.  B,o  =  120. 


Plate  2.  Stereo  pair  illustrating  the  wave  pattern  of  the  cross  spectral  func'ion  at  I  MHz  in  the  x-y  plane  (B  = 
z  -  0).  The  complex  pattern  is  due  to  the  spatial  interference  of  many  oblique  modes.  The  x  and  y  axes  have  been  shilled 
by  +  24  and  -h  10  cm  in  this  display. 
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Plate  3.  Spatial  variations  of  one  or  a  few  wave  vector  components  of  the  cross-spectral  function.  These  are  selected 
from  the  peaks  of  Figure  6  and  inverse  Fourier  transformed,  (u)  Mode  A  with  A.  =  48  cm,  ^,  =  S  cm.  (h)  Mode  A  and 
mode  B  =  10  cm,  =  12.5  cm),  (c)  Modes  A,  B,  C  (2,  =  8.3  cm,  2,  =  10  cm),  and  D  (2,  =  7.14  cm,  k^—  5  cm).  The 
decomposition  allows  one  to  investigate  the  polarization  of  single  modes  by  comparing  orthogonal  tensor  components. 
The  X  and  y  axes  have  been  shifted  by  +24  and  + 10  cm,  respectively,  in  this  display. 
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Plate  4a 


Test  Wave 


^  N  I 


Plate  46 

Plate  4.  Two  dimensional  interferometer  traces  of  the  B,  components  of  a  magnetic  wave  at/=  1  MHz.  The  wave  is 
excited  at  x  =  -20  cm,  y  =  0,  z  =  0,  and  t  =  21.4  la  into  the  reconnection  event  by  a  magnetic  loop  antenna  with  B 
(0, 0,  B.)  (o)  propagating  in  the  x-y  plane;  =  lOG,  (6)  propagating  in  the  x-z  plane. 
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Plate  56 


Plate  5.  Color  stereo  pairs  of  the  three-dimensional  wave  number  spectra  at  /  =  1  MHz.  The  vertical  green  contour 
map  is  of  the  wave  spectra  in  the  k,-k^  plane,  the  lower  contour  map  is  tr,.nsverse  k^-k^  plane.  The  yellow  markers  with 
size  proportional  to  amplitude  are  projections  of  the  peaks  into  3  space.  The  magenta  surfaces  are  the  cold  plasma 
dispersion  curves  for  oblique  whistler  waves  at  B,  =  10  G,  n,  =  4.9  x  10‘Vcm*.  The  vertical  white  line  is  drawn  through 
k,  =  k,  =  0  and  a  blue  diamond  at  the  origin  k  =  0.  (a)  Spectrum  of  <B„,B,2)  illustrating  prolific  mode  structure.  All  fall 
within  the  limiting  dispersion  surface.  (6)  Spectrum  for  <B,,B,2>.  There  are  only  three  peaks  present.  It  is  possible  to  single 
out  one  mode  (also  see  Figure  Sr)  for  polarization  analysis  in  this  case  (i  =  33.2  /«). 
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Plate  6<i 


Plate  6b 


Plate  6.  Electron  distribution  functions,  (a)  At  the  edge  of  the  neutral  sheet  (.>.  =  20  cm,  r  =  0)  where  no  tails  exist,  and 
^h)  in  the  center  of  the  neutral  sheet  (x  =  0,  z  =  0)  where  runaway  electrons  are  observed.  Upper  panels  show  ensemble 
average  values  of/((;|,),/(ei),  the  surface  is  on  a  logarithmic  scale  (I0~  ”  </<  10“  "*  cm’*’  s’,  0  <  |r|  <  7  x  10" 

cm/s).  Bottom  plane  shows  rms  fluctuation  vj/f{0)  =  constant  (3  dB/contour). 
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tensor  components  of  the  CSF.  This  analysis  is  not  possible  in 
a  multimode  pattern  such  as  is  shown  in  Figure  So. 

Prior  to  the  polarization  analysis  we  wish  to  review  some 
basic  properties  of  monochromatic  whistlers  excited  by  an¬ 
tennas  in  the  neutral  sheet. 

Test  Wave  Experiments 

The  CSF  of  the  random  1-MHz  waves  exhibit  good  corre¬ 
lations  through  the  entire  x-y  plane.  Besides  weak  damping, 
this  fact  suggests  that  the  excitation  sources  are  extended  so  as 
to  excite  plane  waves.  For  comparison,  we  now  launch  whis¬ 
tler  waves  at  the  same  frequency  under  identical  conditions 
from  an  antenna  placed  into  the  neutral  sheet.  Such  wave 
experiments  have  been  performed  earlier  in  a  quiescent,  uni¬ 
formly  magnetized  plasma  with  similar  parameters  (tu,  » 

<  £0  <  o>c)  except  without  a  current  sheet  IStenzel,  1975]. 
In  these  experiments  a  magnetic  loop  launched  whistlers  with 
diverging  phase  fronts,  which  causes  an  amplitude  decay  due 
to  the  divergence  of  the  wave  energy.  The  presence  of  a  den¬ 
sity  gradient  causes  considerable  wave  refraction  and  ampli¬ 
tude  variations.  In  the  reconnection  geometry  the  magnetic 
field  is  sheared  (Plate  1)  and  the  density  within  the  neutral 
sheet  is  highly  nonuniform  (part  2).  In  order  to  facilitate  com¬ 
parison  with  the  random  waves,  the  test  signal  is  launched 
near  the  center  of  the  neutral  sheet  during  the  reconnection 
event. 

The  exciter  is  a  6-cm  diameter,  five  turn  loop,  launching  a 
wave  with  B,  «  3  G  at  1  MHz.  The  applied  signal  is  attenu¬ 
ated  and  is  used  as  the  reference  signal  in  the  correlation  setup 
that  is  identical  with  an  interferometer  circuit.  Owing  to  the 
phase  coherence  the  test  wave  signal  is  easily  distinguished 
(tom  the  1-MHz  background  fluctuations.  Surfaces  of  con¬ 
stant  interferometer  output  [ocBfr)  cos  (k  •  r)]  taken  approxi¬ 
mately  30  fis  into  the  reconnection  event  are  shown  in  the  x-z 
plane  in  Plate  46.  Since  the  movable  probe  could  not  be  posi¬ 
tioned  beyond  the  exciter  loop,  the  data  are  symmetrized  in 
the  X  direction  about  the  z  axis  going  through  the  antenna 
location.  The  radiation  pattern  in  the  x-y  plane  along  the 
dominant  axial  magnetic  field  component  B^  is  shown  in 
Plate  4<i.  It  indicates  how  the  nonplanar  phase  fronts  spread 
out  along  and  across  B^.  The  wave  motion  in  time  is  best 
seen  in  a  computer-generated  movie  that  displays  the  mea¬ 
sured  T  dependence  of  Plate  4.  A  dominant  oscillation  with 
rapid  spatial  amplitude  decay  is  observed  in  the  transverse 
plane  (Plate  46).  From  these  results  one  can  estimate  the  typi¬ 
cal  wavelength  to  be  <2>  »  10  cm. 

Wave  number  spectra  generated  by  Fourier  transforming 
the  X,  y  dependence  of  B,  are  shown  in  Figure  9.  Unlike  the 
corresponding  fluctuation  spectra  there  is  one  peak  elongated 
along  Byo  in  the  k^-k,  plane.  The  contours  are  nearly  circular 
in  the  transverse  k,-k,  plane.  The  dominant  mode  in  the  k,-k, 
plane  moves  at  an  angle  of  B=  —32°  with  respect  to  the 
magnetic  field  B^.  In  Plate  46  the  exciter  antenna  at  x  =  —20 
cm  is  always  to  the  left  of  the  receiver  antenna.  The  dominant 
modes  in  the  k„-k,  plane  are  therefore  peaked  at  k,  >  0,  de¬ 
scribing  the  outward  wave  propagation.  The  waves  of  lar^t 
amplitude  have  axial  wavelengths  in  the  range  10  ^  >1,  £  20 
cm,  in  reasonable  agreement  with  the  dispersion  relation  for 
whistler  waves  given  by 

kcV  .  o),^ 

—  I  —  1  —  -  ■  t - - E _ 

(oj  w((0  —  0),  cos  0)  0j((t),  cos  B  —  oi) 

where  cu,,  and  o  are  the  electron  plasma,  cyclotron,  and 
wave  frequencies,  respectively.  The  same  wavelength  regime 


Fig.  10.  Phasor  diagram  for  the  test  wave  generated  by  plotting 
two  components  of  the  CSF  tensor  and  <B,iB,2>  as  a 

liinction  of  delay  time  r,  0  <  r  <  3.0  /ts,  rotated  into  the  frame  in 
which  k  1  (fc,  £  0).  The  signal  used  to  launch  the  wave  is  the 
common  reference  signal  B,,  in  the  interferometer.  The  obliquely 
propagating  mode  is  right-band  circularly  polarized  (c  =  1,  =  90°), 

as  seen  by  the  projection  onto  the  grid  below.  The  arrows  indicate  the 
magnitude  and  direction  of  B(r)  at  several  of  the  64  measured  delay 
times.  The  dots  on  the  trajectory  indicate  that  the  measurement  of 
polarization  was  made  on  the  second  cycle  and  are  projected  below. 


was  found  from  the  CSF  analysis  of  the  fluctuations  that  es¬ 
tablish  the  connection  between  test  whistler  waves  and  the 
unstable  modes. 

The  test  wave  exciter  does  not  generate  plane  waves.  Owing 
to  geometric  spread  of  the  wave  energy  from  a  point  source, 
the  wave  amplitude  decreases  rapidly  with  distance  from  the 
exciter.  In  contrast,  the  unstable  waves  exhibit  a  small  ampli¬ 
tude  decay.  One  can  therefore  conclude  that  the  random 
waves  are  generated  by  extended  sources;  for  example,  long 
current  filaments  within  the  neutral  sheet. 

Wave  Polarization 

We  consider  the  theoretical  dispersion  and  polarization 
properties  of  whistlers  in  the  quasi-longitudinal  approxi¬ 
mation  [Ratcliffe,  1959].  It  is  valid  when  the  angle  of  wave 
propagation  satisfies 

sin^  0  2tD,^ 

- o  <  (8) 

cos  0  COCOj 

For  our  parameters  the  right-hand  side  exceeds  10‘  so  that  the 
approximation  is  valid  for  essentially  all  angles  (0  <  89.9°).  In 
the  quasi-longitudinal  approximation,  the  polarization  of  the 
wave  magnetic  field  in  tae  plane  perpendicular  to  k  is  right- 
hand  circular  irrespective  of  the  angle  of  wave  propagation 
IHeIttwell,  1965;  Fonmsano  and  Kennel,  1969].  The  polariza¬ 
tion  of  the  electric  fields  is  different  and  will  not  be  considered 
here. 

First,  we  analyze  the  polarization  of  the  test  wave.  The  peak 
mode  at  k  •  (0.35  cm "  —  0.286  cm"  *)  is  diosen  from  the 

k,-k,  spectrum  of  Figure  9.  This  mode  fairly  well  satisfies  the 
dispersion  relation  (equation  (7))  for  n,  =  1.3  x  10'^  cm"’ 
±20%,  Byo  •  12  G,  the  experimental  conditions. 

For  analyzing  the  polarization  we  consider  two  orthogonal 
magnetic  field  components,  B,  and  B,.  The  third  component 
B,  does  not  have  to  be  measured  since  it  is  determined  by 
k  •  B  -  0.  Both  components  are  digitally  filtered  at  k„  2x/18 
cm"'  and  k,  =  -2x/22  cm" ',  and  the  signals  are  subsequent- 
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Fig.  II.  Stereo  pairs  of  phasors  generated  from  random  waves.  The  display  is  as  in  Figure  10;  however,  is  the 
random  wave  signal  to  the  common  reference  probe,  (o)  Mode  with  s  >=  0.79  cm~'.  It  is  nearly  circularly  polarized, 
right  handed,  and  satisfies  the  whistler  dispersion  relation.  There  is  insuf^ient  information  to  transform  oblique  modes 
properly,  such  as  is  shown  in  (h),  S  ss  26°  to  the  frame  with  k  J.  B'.  The  phasors  exhibit  bizarre  motions  in  this  case. 


ly  inverse  Fourier  transformed  so  ds  to  obtain 
B(a),  k)  B(o),  k) 

Amplitude  and  phase  of  each  component  is  retained  as  a  func¬ 
tion  of  delay  time  t  so  that  we  can  plot  the  phasor  Be¬ 
cause  we  consideT  the  magnetic  field  polarization  in  the  plane 
perpendicular  to  k  the  field  components  in  the  laboratory 
coordinate  system  (x,  y,  z)  must  be  projected  into  a  rotated 
frame  (x',  z')  with  y  along  k.  This  transformation  yields 

B.'  -  BJeoa  0  B.'  =  B,  B,'  =  0  (9) 

where  the  primes  refer  to  the  rotated  coordinate  system. 

Figure  10  is  a  phasor  diagram  By,,'  for  the  test  wave.  The 
vector  components  are  noimalized  in  amplitude  as  discussed 
in  the  previous  section,  and  one  component  is  delayed  with 
respect  to  the  other  one  by  SO  ns  wh^  was  the  digitization 
time,  hence  corresponds  to  an  error  of  only  one  syn¬ 
chronization  pulse.  The  data  refer  to  a  start  time  of  33.2  fu 
into  the  reconneGtion  event  and  display  the  phasor  B,,,'  for 
delay  tunes  0  £  r  ^  3  ps.  The  vector  tip  sweeps  out  a  right- 
banM  tpktl  in  ihlie  irtiose  projection  into  the  S'  (dane  is  a 
drde.  At  five  discrete  time  steps  the  vector  By,,'  is  indicated. 
The  figure  also  ptesentt  values  for  s,  the  ratio  of  the  sonima- 
jor  to  seminiinor  axis,  and  the  phaee  difference  between  the 
vector  components  measured  between  successive  maaiina 

B _ and  hi  one  of  the  cycles.  All  the  test  wave  modes  in 

the  vichiity  of  the  spectral  peak  (Figure  9)  are  right-laaded 
circular  wtotler  waves. 


The  random  modes  have,  in  general,  three  wave  vector  com¬ 
ponents  that  must  be  known  in  order  to  find  the  polarization 
in  the  plane  orthogonal  to  k.  This  would  require  that  CSF  be 
measured  throughout  a  volume  so  that  one  can  perform  a 
three-dimensional  spatial  Fourier  transformation.  However,  it 
is  often  sufficient  to  consider  two  orthogonal  planes  if  certain 
symmetries  in  the  wave  propagation  exists.  For  example,  one 
may  consider  the  polarization  of  all  the  waves  along  the  lines 
k,  =  ±0.126  cm"'  (A,  =  ±50  cm)  in  the  k,-k,  plane.  Figure 
8a.  These  are  modes  that  travel  at  small  angles  with  respect  to 
the  uniform  magnetic  field  The  spectrum  in  the  k,-k, 
plane  is  roughly  symmetric  about  k^  and  as  an  approximation 
we  choose  for  the  average  wave  number  <k,>  «  0  for  purpose 
of  investigating  the  ptiarization  of  modes  in  the  k„  k,  plane. 
The  polarization  at  different  values  k,  should  be  right-hand 
circular  for  waves  that  fit  the  dispersion  relation  (equation  (7)). 
More  than  80%  of  the  modes  tested  are  confirmed  to  be 
right-hand  polarized.  For  example.  Figure  1  la  shows  a  phasor 
diagram  for  a  random  mode  with  k  2:  k,  ~  0.63  cm"'  with 
dbptidty  s  •  0.94  in  the  first  cyde,  i.e.,  with  nearly  right-hand 
circular  polarization.  Most  of  the  random  waves  exhibit 
dianges  in  s  and  orientation  of  the  nugor  axis  as  a  function  of 
delay  time  t.  The  CSF  (Figure  6)  and  resulting  k  spectra  are 
weakly  time  dependent,  and  it  is  possible  that  the  tiurd  com¬ 
ponent  k,  entoi  in  time  and  affects  the  calculated  pdariza- 
tion.  Whether  indeed  all  the  waves  ate  right-hand  polarized 
can  only  be  resolved  by  performing  a  ftiUy  three-dimensional 
Fourier  analysis  that  permits  a  proper  transformation  of  B 
into  the  (done  orthogonal  to  k.  We  interpret  the  neariy  dreu- 
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lar  poUrizations,  wen  mainly  for  paraUd  wavekngths  of  order 
10  cm  as  casM  of  whiMler  wavm  with  h, «  Oi 

EzaminatioB  of  aMdaa  that  am  more  ohlk|ue  reiaforoes  this 
picture.  Figure  llh  is  a  phasor  of  a  wave  propagatiiig  at 
$  =  26°  with  respect  to  die  magnetic  fidd  The  phasor 
diagram  is  very  complicated  and  is  best  ^tpracialed  in  a  ster- 
eo8co|M  view,  as  shown.  The  polarization  appears  highly  el¬ 
liptical  and  is  seen  to  reverse  as  a  function  oft.  The  projection 
onto  the  lower  grid  resonUes  the  fluctuating  ccmponent  of 
the  hodogram  m  Figure  2e  as  if  everything  save  the  loops 
visible  in  Figure  7b  have  been  averaged  away.  TUa  type  of 
phasor  is  re|»esentative  of  all  waves  propagating  at  large 
angles  with  respect  to  B^.  This  does  not  necessarily  mean  that 
these  waves  are  dhptical  but  that  there  may  be  inaufficient 
information  to  transform  B  into  the  proper  frame.  If  there 
were  more  data  points  available,  a  statistical  analysis  for  the 
polarization  as  a  function  of  propagation  angle  could  be 
undertaken  in  order  to  check  this  assumption.  The  present 
data,  however,  strongly  suggest  that  the  random  waves  are 
whistler  waves  propagating  at  all  angles  up  to  the  resonance 
cone  (0  *  cos"  ’  (<u/(a,)  =  88°). 

The  polarization  sampling  procedure  was  relegated  to 
waves  moving,  on  the  average,  at  small  angles  with  respect  to 
the  background  magnetic  field.  The  spectrum  of  in 

the  k^-k,  plane  has  one  intriguing  highly  oblique  mode  (Figure 
8c).  An  estimation  of  its  polarization  is  possible,  but  best  fol¬ 
lows  a  discussion  of  the  three-dimensional  spectra. 

Three-Dimensional  Wave  Propagation 

Wave  number  data  acquired  in  two  orthogonal  (k„-k,)  and 
{k,-k,)  planes  make  a  three-dimensional  reconstruction  of  the 
wave  vector  structure  possible.  Such  a  diagram  is  shown  in 
Plate  5  as  a  stereoscopic  pair  for  <B,,B,]>.  Wave  number 
spectra  are  shown  as  contour  plots  projected  to  the  bottom 
Had  to  the  side  of  a  rectangular  volume  in  k  space  (|kj  ^  0.69 
cm~S  |k^  |kj  ^  1.35  cm~‘).  Each  point  a  in  the  volume  is 
tested  to  insure  the  mode  associated  with  it  is  well  above  the 
noise;  CSF  [(fc„  k,)„  {k„  k,)J  ^  0.4S,  where  the  largest  value 
of  the  combined  spectra  is  normalized  to  be  1.00.  The  mode  is 
then  projected  into  three-dimensional  k  space  as  a  yellow 
marker  with  size  proportional  to  the  amplitude  to  which  it 
corresponds.  Also  drawn  in  magenta  are  the  dispersion  sur¬ 
faces  lc(d)  (equation  (7))  for  obliquely  propagating  whistler 
waves  at  a  fixed  density,  in  argon,  n,  =  4.9  x  10‘Vcm^>  iURg- 
netic  field,  =  10  G,  and  frequency  to/2n  =  1  MHz.  A  verti¬ 
cal  line  through  the  symmetry  axis  (k,  ^  k,  =  0)  is  drawn  to 
aid  the  eye.  A  blue  diamond  marks  the  origin  in  k  space. 

The  density  within  the  neutral  sheet  is  nonuniform  (part  2) 
and  may  be  2.5  times  larger  at  the  second-order  contact  points 
of  the  separatrix  than  in  the  center  of  the  neutral  sheet.  Waves 
in  lower  density  regions  in  turn  give  rise  to  points  in  Plate  5a 
that  lie  doser  to  k,  s  0  for  nearly  perpendicular  propagation. 
The  magnetic  fidd  is  constant  to  within  10%  along  the  y 
direction  but  increases  in  the  x-z  plane  as  one  moves  away 
from  the  neutral  sheet.  In  this  situation  the  wave  dispersion 
surfaces  are  mpected  to  sweep  through  a  volume.  The  curves 
that  are  drawn  in  Rate  5a  may  be  intmpieted  as  a  reasonable 
outer  envdope  for  these  experimental  parameters.  That  the 
majority  of  the  modes  fall  within  this  voiuine  of  k  space  is 
strong  evidence  diat  the  magnetic  turbulenoe  is  comprised  of 
whistler  waves. 

The  large  number  of  modes  in  wave  number  space  pre¬ 
cludes  any  fiirtber  investigation  of  the  CSF  component  shown 
m  Plate  5a.  This  is,  however,  not  always  the  case.  There  are  far 


fewer  modes  observed  for  tensor  element  <B,,B,2>  as  illus¬ 
trated  in  Plate  56.  The  dominant  mode  located  at  k^  »  —0.62 
cm"*,  k,  »  0.26  cm"‘,  and  k,  =  —0.26  cm~‘  is  an  oUique 
mode  with  its  counterpart  in  tensor  component  <B,iBx2>-  Its 
peaks  are  dearly  seen  in  the  contour  lines  projected  onto  the 
lower  and  side  [rianes. 

If  the  fluctuating  magnetic  field  C^F  had  been  measured 
everywhere  within  a  volume  of  real  space,  the  polarization 
analysis  would  be  tedious  but  straightforward.  The  three- 
dimensional  Fourier  transform  into  wave  number  space  and 
filtration  of  the  mode  of  interest  would  proceed  analogous  to 
the  technique  described  in  the  previous  section.  After  transfor¬ 
mation  into  a  plane  perpendicular  to  k,  the  polarization  may 
be  traced. 

Since  the  CSF  is  mapped  on  ordiogonal  planes  only,  the 
following  procedure  is  used  to  estimate  the  polarization.  Both 
planes  have  an  axis  (kj  in  common.  The  CSF  modes  are 
filtered  out  of  the  orthogonal  k  space  planes  and  the  compo- 
nente  vectorially  summed,  i.e,  <B„B,2>  =  <B„B,2>(k,p  k,) 
+  <B„B,2>(k,p  k,),  where  ot  =  z,  x.  The  data  from  each  of  the 
orthogonal  planes  were  of  comparable  magnitude,  but  their 
initial  phase  had  to  be  equalized  since  they  were  acquired  on 
separate  runs.  This  did  not  affect  the  polarization  calculation. 
Finally,  we  note  that  the  probe  in  common  for  this  CSF  pair 
was  the  movable  probe  that  is  equivalent  to  replacing  t  by 


The  data  were  transformed  by  two  successive  rotations  in  k 
space  such  that  k"  =  (0,  k,",  0)  and  8,2’  =  0  since  k  >  B  =  0. 
(For  simplicity  we  drop  the  prefix  B,,.)  The  transformation  is 
given  by 


where 


B/  =  B,/cos  9,  +  B,  ^  sin 


k,  sin  02 

(k,'  +  k/)‘'* 


-I-  cos  02 


) 


(10) 


Here  k  is  the  wave  number  in  the  original  laboratory  frame. 
Note  that  in  the  x-y  plane  k,  =  0  and  (10)  reduces  to  the 
transformation  given  by  (9). 


Fig.  12.  Phasor  generated  from  single  CSF  modes  digitally  fil- 
wred  from  the  three  dimensional  k  spectra  (as  shown  in  Plate  Sh)  and 
transformed  sudi  that  k  1  B.^.  The  hi^y  obUque  mode  a  «  oos"' 
kf/k  »  64°  is  a  right-hand  farmed  a^  nearly  circular  poiarieed 
whistler  (t  »  33.2 /a). 
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The  transformed  tfota  plotted  as  a  phasor  u  Figttre  12  idei- 
tify  the  dominuit  mode  to  be  a  ri^t-haad,  nearly  ckcular 
polariaBd  whistler  wave.  Althm^  the  agreement  here  could 
he  fortuitous,  we  view  this  as  being  uniflcefy  in  li^t  of  the 
nperunent  in  swveyug  modes  in  which  k  is  not  known  (sudi 
as  is  shown  in  Figure  1  lb). 

The  wtustlsr  waves  are  not  the  only  component  of  the  tur¬ 
bulence  within  the  neutral  sheet  We  have  performed  pralinii- 
nary  fluctuation  studies  at  highre  frequencies  and  observed  a 
wide  range  of  dectrostatic  modes. 

Ion  Sound  and  Lanomudi  TuanuLENCs 

Using  two  electrostatic  probes,  we  have  measured  with 
analog  corrdation  techniquet  the  frequency  and  wave  number 
spectrum  of  modes  along  within  the  neutral  sheet  The 
results  are  summarized  in  Figure  13.  The  cross  spectral  func¬ 
tion  at  2  MHz  eiubles  the  straightforward  evaluation  of  a 
single  wavelength  (unlike  the  whistler  case).  From  diflerent 
measurements  of  m  and  k  we  obtain  a  linear  dispersion 
characteristic  for  ion  acoustic  waves.  It  is  well  known  that 
short  wavelength  ion  sound  modes  as  can 

produce  significant  scattering  of  dectrons  which  manifests 
itself  as  an  enhanced  resistivity  (part  4). 

The  microwave  emission  GHz)  from  the  recon¬ 

nection  event  is  monitored  with  a  horn  antenna  located  out¬ 
side  the  plasma  column  and  connected  to  a  sensitive  micro- 
wave  receiver.  Figure  14  shows  a  large  burst  of  electro¬ 
magnetic  radiation  during  reconnection.  In  this  case  the  event 
is  triggered  late  in  the  discharge  to  make  the  signal  easily 
distinguishable  from  that  due  to  the  primary  (ionizing)  dec¬ 
trons.  In  a  separate  experiment  IWhelan  and  Stemel,  1981], 
electromagnetic  radiation  at  the  plasma  frequency  is  observed 
to  be  generated  when  electron  plasma  waves  scattered  off  ion 
acoustic  fluctuations.  The  Langmuir  waves  were  driven  unsta¬ 
ble  by  an  electron  beam. 


0  to  20  » 


FR£aU€IKr  f  (MHz) 


Fig.  13.  Characttrittici  of  dectrostatic  fluctuations  duthtg  recon- 
aection.  (a)  The  frequency  spectrum  of  density  fluctuations  extends  iq> 
to  the  ion  plasma  frequency  lategraMd  fluctuation  isvnb  aie  on 
the  order  all  at  5%.  (h)  Wave  aiimher  spretrum  of  low- 

bemieney  fluctMtiow  exhiiiits  the  diepersieu  of  the  ion  sound  waves. 
The  spatial  noise  properties  ate  obtafosd  from  two-probe  crom  power 
spectrd  fmetions  C,]  (dr),  an  exmnplc  of  which  a  inserted. 


MtCROW/WE  EMISSION 


12  GHz 


Timt  t 


Reconnection 


Fig.  14.  Microwave  emission  (arb  units)  at  the  plasma  frequency 
observed  with  an  external  horn  connected  to  a  radiometer.  The  back¬ 
ground  radiation  is  due  to  primary  electrons  within  the  dc  discharge. 
For  purpose  of  darity  the  reconnection  event  is  switched  on  12  ms 
after  the  discharge  is  initiated.  The  large  pulse  of  radiation  is  pre¬ 
sumed  to  be  associated  with  Langmuir  turbulence  driven  by  electrons 
accelerated  within  the  neutral  current  sheet. 


In  the  recomtection  event  a  runaway  tail  of  electrons,  con¬ 
fined  to  the  neutral  sheet,  is  produced  by  the  induced  electric 
field.  Langmuir  waves  grow  at  the  expense  of  the  free  energy 
in  the  tail,  scatter  off  ion  acoustic  waves  or  convert  on  density 
gradients  to  generate  the  observed  microwave  emission. 

Electron  DiSTRiBunoN  Functions 

Non-Maxwellian  particle  distribution  functions  may  pro¬ 
vide  conditions  for  plasma  wave  growth.  If  these  waves  reach, 
large  amplitudes  particles  in  turn  may  be  scattered  or  trapped 
by  them.  If  the  full  distribution  function  were  known, /(v,  r,  t), 
most  instabilities  may  be  predicted  theoretically  or  by  simula¬ 
tions.  Although  a  complete  measurement  of  /  is  a  massive 
undertaking  in  data  acquisition,  several  steps  toward  it  have 
been  made.  One  would  like  to  compare  the  observed  wave 
phenomena  and  particle  scattering  with  distribution  functions 
to  see  if  they  are  consistent.  A  strong  modification  of  the 
distribution  at  a  given  phase  velocity  may  be  reflected  by 
strong  wave  activity  there.  In  this  section  some  preliminary 
measurements  will  be  discussed;  details  will  be  presented  in  a 
subsequent  article  in  this  journal. 

The  electron  distribution  functions  are  measiued  with  a 
novel  velocity  analyzer  [Stenzef  et  al.,  1983<i,  c]  that  transmits 
electrons  within  a  small  acceptance  angle  (AO  =  0.036  sr).  The 
essentially  one-dimensional  velocity  of  the  transmitted  elec¬ 
trons  is  analyzed  with  a  gridded  retarding  potential  energy 
analyzer.  The  analyzer’s  mounting  permits  a  continuous  rota¬ 
tion  about  the  two  polar  angles  0,  so  as  to  obtain  f(v,  0, 
^)»/(t).  This  involves  repeated  measurements  at  typically 
200  different  an^.  The  device  may  be  moved  to  different 
spatial  positions  within  tire  reconnection  r^ion.  Temporal 
resohition  is  provided  by  a  fast  voltage  sweep  (At »  2  /a)  on 
the  Mtadyar  grid.  In  prindifie,  it  is  posrible  to  determine  the 
full  elcotm  distribution  function  fjr,  r,  t). 

The  analyzer  output  is  digitized  with  a  100-MHz  A/D  con¬ 
verts.  At  every  a^fo  10  current  voltage  traces  are  recorded. 
The  on-line  may  processor  (fifferentiates  the  I-V  curves  and 
smooths  out  the  d^tal  nmse.  It  then  calculates  the  ifistri- 
butiott  ftnetioa/a  wdl  as  the  rms  fluctuations  at  every  angfo. 
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Tke  dMtrUwtioii  Amctioos  aie  displayed  at  three- 
dknentioiial  steieo  puts,  for  ease  of  viewing,  m  Plale  6.  The 
upper  tuifooe  shows  the  ditfribution  fonctfon  f(v,i,  ej  taken 
3U  t»  into  Oe  reconnactkw  event,  at  the  e^  of  the  neutral 
sheet  (z  -22.5  cm,  PhUe  6a)  and  widun  it  (z  -  7.5  cm,  Plate 
Here  sn  denotes  vdocity  along  the  uniform  guide  field 
The  velodty  range  is  tP±,  Sut  ^  d  x  10*  cm/s.  The  distribiition 
fimction  is  plotted  logarithmically  over  three  orders  of  mag¬ 
nitude,  10"*’s:/(*,  z«7.5  cm)il0"“  cm~*  s*, /(e-0, 
z  •  -215  cm)  -  1.3  /(»  -  0,  z  «  7.5  cm).  In  drawing  these 
graphs  it  is  ttot  tbs  curve  is  symmetric  about  the 

asmuthal  ai^  d-  Full  three-dimensioniJ  data  sets  indicate 
that  the  symmetry  is  only  approximate.  The  dbtribution  Amc- 
tidn  is  examined  by  using  computer  generated  films  of 
successive  surfaces  of  /(v)  so  constant.  On  the  plane  projected 
immediately  below  the  distribution  surface  are  ten  (green)  con¬ 
tours  log  /(v) «  const  (3  dB/contour).  We  observe  a  high- 
energy  electron  tail  (Plate  66)  that  carries  the  plasma  current 
and  extends  to  energies  in  excess  of  75  eV.  Since  most  of  the 
energetic  electrons  in  the  discharge  have  energies  eV^  ^  45 
eV,  we  conclude  they  have  been  iteely  acoderated  by  the  dec- 
tric  field. 

The  shaded  violet  projections  on  the  rear  surfaces  indicate 
what  would  be  measured  by  probes  directly  along  {f^^)  across 
{fj  the  uniform  bacfcgrou^  field.  Anisotropic  distributioiu 
caimot  be  resolved  by  using  single-sided  planar  Langmuir 
probes  sinoe  they  integrate  f(v,  0, 4)  over  A0,  Ad  a  in  vdoc- 
ity  space. 

At  the  current  sheet  edge  near  the  contact  points  of  the  Y 
shaped  separatrix  (Plate  1,  Figure  5)  there  is  a  vertical  mag¬ 
netic  field  component  (B,  ~  2  G)  that  prevents  runaway  tail 
formation  (Plate  6a). 

The  velodty  space  rms  fluctuations  of  the  distribution  func¬ 
tions  are  shown  as  blue  contours  in  the  lowest  plane  of  Pbte 
6.  The  Doppler  shiAed  whistlers  have  phase  velodties  |C|||  » 
|<i)  -  <a,|/k,|  (3:2.5  X  10*  cm/s)  that  fall  into  the  velodty  range 
of  large  rms  fluctuations.  The  phase  vdodty  of  beam  driven 
Langmuir  waves  also  (all  into  this  regime  m/fcn »  ti  10* 
cm/s.  No  detailed  measurement  of  this  high-frequency  turbu¬ 
lence  has  been  made  thus  far. 

The  ion  acoustic  modes  have  phase  vdodties  <o/k  =  c,  =  3.5 
X  10*  cm/s  that  place  them  near  the  origin  of  the  plots  in 
Plate  6,  where  measurement  errors  limit  the  resolution  of/(v). 

Summary  and  Conclusions 

We  have  approached  sorting  out  the  complex  net  of  inter¬ 
woven  processes  that  constitute  magnetic  fidd  line  reconnec¬ 
tion  with  measurements  taken  from  three  perspectives.  The 
first  levd  is  macroscopic  in  which  ensemble  averaged  proper¬ 
ties  such  as  density,  temperature,  and  embedded  electric  and 
magnetic  fidds  were  mapped  out  as  a  function  of  space  and 
time.  This  is  somea^t  analogous  to  Ae  fluid  amwoa^  in 
thcQwsfawi  pl^dcs,  and  constitutes  the  bulk  of  meararements 
deserfoed  in  JGR  ptpen  1-5.  ABhough  tuthulHine  may  aumi- 
fest  itself  in  the  form  of  anomafous  resistivity  and  enhanced 
fluid  scattering  (ports  3  and  4),  a  seormd,  statistical  approaA  is 
necesrary  to  probe  the  underlying  structure  of  it  hi  this  work 
extensive  digital  iteta  aequbWon  and  datistical  oompulntioa 
have  been  brou^  to  bear  on  magnetic  turboleiioe  and  en¬ 
abled  the  mdblAttg  of  the  fitequnaey  and  wave  number  spec- 
tnnn  of  uaeubk  whistler  waves.  ^atidicBl  studies  of  the  wave 
pofarfasation  vis  i  vh  the  cross  special  tensor  indicate  the 
waves  are  right  handed  and  nearly  orcular.  Analyris  of  whis¬ 


tler  polarizatidtts  made  on  the  earth’s  railaoe  of  waves  whidi 
entered  dm  irmoaphere  [Twruda  and  /keda,  1979]  are  also 
surveys  hom  vAfoh  the  direction  of  the  incoming 

wave  must  also  be  edimaleA  The  propogadon  dueetkm  k  and 

pohkrization  for  fiuity  monochromatic  waves  have  been  asoer- 
tam^  from  hOdoiriuiis  and  flux  gate  magnetometers  on  the 
ISEE  satellites  is  they  traversed  the  bow  shock  [Hqppr  et  ol., 
19823-  Corrdatkm  of  niagaetic  fiehls  measured  by  the  sateliite 
pair  was  used  to  find  the  phase  vdodty,  again  for  the  uncom¬ 
mon  case  of  neariy  numoduomatic  signals.  Dispersive  ^ects 
make  such  two  point  oorrelarions  difficult  since  the  spacecraft 
are  always  more  than  several  wavelengths  apart.  Neverthdess, 
these  waves  were  found  to  be  oblique  whistlers  driven  unsta¬ 
ble  by  ion  beams  reflected  off  the  bow  shock.  Whistlers  seem 
to  be  ubiquitous  in  space  plasmas,  having  thus  far  been 
observed  in  the  bow  shock,  auroras,  plasma  sheet,  magneto- 
sheath,  and  solar  wind,  as  well  as  in  Ae  environs  of  Jupiter. 

In  our  laboratoiy  plasma  we  have  established  a  magnetic 
field  topology  similar  to  that  of  magnetotail  models.  A  possi- 
bte  difference  is  our  axial  field  component  which  at  times 
may  exist  in  Ae  magnetotail  but  is  usually  ignored  in  theoreti¬ 
cal  and  numerical  models  [Somierwp,  1979].  It  establishes 
magnetized  electrons  in  Ae  current  sheet  and  whistler  eigen- 
modes  throughout  the  dissipation  region.  Our  observations 
Aow  Aat  Ae  whistler  wave  activity  correlates  wiA  the  neu¬ 
tral  sheet  current;  for  exami^  the  wave  intensity  nearly  van¬ 
ishes  at  current  reversal  (/,  a:  0  at  t  »  80  fa). 

In  this  experiment,  we  used  two  probe  cross  correlation 
techniques  for  analyzing  magnetic  turbulence.  This  powerful 
technique  has  not  yet  been  employed  in  space  physics  since  it 
requires  moAer-daughter  satellite  systems.  From  the  corre¬ 
lation  measurements,  we  have  established  Ae  wave  dispersion 
ai(k)  and  magnetic  field  polarization  B(t).  This  was  feasible 
through  the  use  of  Agital  data  processing  involving  fast  Cou¬ 
rier  transforms  and  digital  filtering  on  Ae  spatial  and  tempo¬ 
ral  dependence  of  magnetic  fluctuations.  A  dominant  mode 
found  in  two  components  of  the  CSF  lent  itself  particularly  to 
examine  the  B  field  polarization  that  was  found  to  be  right- 
handed.  A  procedure  to  estimate  the  wave  numbers  in  three 
dimensions  for  this  right-hand  circularly  polarized  wave  as 
well  as  the  oAer  components  of  the  noise  enabled  comparison 
with  the  Aspersion  surfaces  of  Ae  whistlers.  Full  three- 
Amensional  analysis  will  be  possible  with  increased  computer 
memory  and  mass  storage.  Higher  order  correlations  have 
already  been  used  to  estimate  wave-wave  coujAng  and  trans¬ 
port  [Powers,  1974;  Roth  et  a/.,  1978]  in  a  simpler  geometry 
and  could  be  extended  to  this  case. 

Having  established  that  the  observed  random  waves  are 
oblique  whistlers  sets  Ae  stage  for  future  experiments  Aat  will 
explore  how  (or  ift  the  whistters  feed  back  on  Ae  reconnection 
process  via  particle  dynamics.  The  wave  frequendes  are  well 
above  the  ion  but  b^w  Ae  electron  cyclotron  frequency: 
fOd  *  e.>|j|  a:  (0  <  o>„  <  (o^  «.  (Opr  In  this  frequency  r^une 
km  motion  is  affected  by  the  wave  and  is  resonant  wiA  it  at 
0)  =  <oui-  'Fite  ions  may  also  respond  to  space  dtarge  fields 
that,  result  from  nonhoear  interaction  at  electrons  wiA  Ae 
whistlers.  It  is,  however,  not  yet  clear  if  Ae  whktlers  are  re¬ 
sponsible  for  ^  large  km  scattering  |xeviouriy  observed  (part 
3).  Also  unexpiainod  is  Ae  observation  that  the  number  of 
modes  in  k  Space  dilliers  for  dements  in  Ae  CSF  tensor.  This 
presumaUy  depends  upon  Ae  details  of  growA  and  satu¬ 
ration  Aat  have  not  as  yet  bem  expioted  experimentafiy.  To 
Ae  authors’  knowled^  no  theory  for  whisriert  under  these 
conditions  exists. 
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A  stronfer  sutemeot  can  be  matte  about  the  electrons 
whose  Dopfiier  shifted  phase  veloctty  e||  -  (<»  -  <e„ViE||  (2-S  x 
10*  aa/»)  fidls  into  the  veiocity  ^aoe  oecupied  by  the  fast 
electnms  (Vlate  6^  It  k  hkdy  these  ekctrons  provi^  the  free 
cnctfy  to  drive  the  waves  umtable  and  in  turn  are  scattered 
by  them.  A  possiUe  mechanism  could  be  filamentation  within 
tte  cmrent  sheet  (part  S)  in  adiich  individiial  current  channels 
act  as  localized  antennas.  Current  fihunentatioa  and  turbu¬ 
lence  have  been  observed  numerically  in  an  MHD  simulation 
[Afatthaeus  end  Montgomery,  1981]. 

In  an  interesting  partide  simulation  of  a  low  average  beta 
{fi  —  »  0.03)  neutral  sheet  ITaiuka  and  Sato,  1981] 

electron  heating,  anomalous  resistivity,  and  particle  transport 
developed  within  a  plasma  with  initial  electron  and  ion  drifts. 
Lower  hybrid  waves  were  observed;  the  nonlinear  devel- 
Ofnnent  includes  both  electromagnetic  modes  and  distortions 
in  the  magnetic  field. 

Spacecraft  observations  of  scattering  of  energetic  electrons 
by  whistlers  have  been  made  in  the  subauroral  magnetosphoe 
[Carpenter  and  LaBeUe,  1982].  Electrons  trapped  by  obliquely 
propagating  whistlers  may  be,  in  turn,  responsible  for  the  gen¬ 
eration  of  fast  electron  beams  (H.  Matsumoto,  private  com¬ 
munication,  1983). 

It  is  possible  to  unravel  the  interplay  of  particles  and  waves. 
It  must  be  done  at  the  most  fundamental  level,  the  detailed 
measurement  of  particle  distribution  functions  and  the  phase 
space  correlation  of  wave-particle  data.  As  reported  here,  tech¬ 
niques  have  been  developed  to  measure  /(v,  r,  t)  and  its  veloci¬ 
ty  space  fluctuations,  wUch  lay  the  groundwork  for  the  next 
phase  of  this  investigation.  As  data  become  available,  rapid 
theoretical  progress  in  turbulence  will  surely  follow  [Montgo¬ 
mery,  1982]. 
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Measurement  and  Instability  Analysis  of  Three-Dimensional 
Anisotropic  Electron  Distribution  Functions 

W.  Gekelman  and  R.  L.  Stenzel 

University  of  California,  Department  of  Physics.  Los  Angeles,  California  90024 
(Received  14  February  I98S) 

A  velocity  analyzer  with  high  angular  resolution  (</ft/4ir-10~'')  has  been  used  to  measure  the 
electron  distribution  function  f(v.9,^)  in  a  neutral  sheet  of  a  large  laboratory  plasma.  Contours 
of /(v)  -const  form  surfaces  in  velocity  space.  As  /(v)  decreases,  anisotropic  tails  form  detached 
surfaces  which  send  filaments  out  toward  the  main  body.  The  highly  anisotropic  structures  occur  at 
the  Doppler-shifted  phase  velocity  of  observed  whistler-wave  turbulence.  Instability  analysis  of 
/( v)  is  consistent  with  the  observed  whistler-  and  Langmuir-wave  turbulence. 

PACS  numbers:  52.25. -b,  52.35. Mw,  52.35. Ra 


The  interplay  between  waves  and  particles  is  one  of 
the  most  fundamental  processes  in  plasma  physics.  In 
a  well-diagnosed  laboratory  plasm.a  it  is  now  possible  to 
measure  the  dispersion  and  propagation  uidc)  of  un¬ 
stable  waves  as  well  as  the  three-dimensional  particle 
distribution  functions.  Such  basic  measurements  are 
of  general  interest  in  and  may  be  extrapolated  to 
fusion  and  space  plasmas  (loss-cone  distributions, 
beams,  and  double  layers). 

The  experiments  are  performed  in  a  large  linear- 
discharge  plasma  ((«,)  =  10'Vcm^  kT,  =  10  eV,  He 
pressure  =  10“*  Torr,  1-m  diam  x2-m  length).'  In  ad¬ 
dition  to  a  dc  axial  magnetic  field  Boy  — 10  G,  a 
transverse  magnetic  Held  with  a  vacuum  neutral  line 
is  applied  by  the  pulsing  of  parallel 
currents  (20  kA,  100  /xs)  through  two  external  elec¬ 
trodes.  In  vacuum  the  topology  of  this  field  is  an  X 
point;  the  field  lines  in  (he  center  of  any  transverse 
plane  intersect  at  90°.  When  plasma  is  present,  how¬ 
ever,  the  changing  magnetic  field  induces  a  current 
which  assumes  the  form  of  a  sheet  ( Ajr  -40  cm  wide, 
Ar-5  cm  half  width,  /,  =  1000  A).  The  measured 
ensemble-averaged  magnetic-field  topology  in  the 
plasma  [Fig.  (la)]  is  that  of  a  classical  neutral  sheet. 
The  electron-distribution  function  is  measured  with  a 
highly  directional  velocity  analyzer^  which  can  be 
placed  at  any  position  x  within  the  neutral  sheet  or  in 
the  strong  transverse-field  regions  [|x|  ^15  cm]  on 
its  border.  The  mounting  of  the  analyzer  on  a  double 
shaft  provides  for  variation  in  the  polar  angle  9  while 
the  probe  rotation  yields  the  variation  in  tft.  Every 
second  (the  experimental  repetition  time)  the  digitized 
current- voltage  trace  ( - 10  ns)  is  transported  by 
a  high-speed  network  from  the  data-acquisition  com¬ 
puter  (LSI  11/23)  to  a  VAX  11/750  computer  and 
on-line  array  processor  (MAP  300).  The  data  are 
smoothed  by  use  of  convolution  techniques  and  dif¬ 
ferentiated;  the  result  is  proportional  to  The 

plasma  potential,  which  can  fluctuate  from  shot  to  shot 
within  the  turbulent  current  sheet,  is  calculated  and 
/( u)  is  shifted  accordingly.  The  procedure  developed 


ensures  that  the  distributions,  averaged  over  an  en¬ 
semble  of  ten  shots,  measured  at  432  angles  match  at 
the  origin.  After  completion  of  a  run  the  data  can  be 
shipped  over  a  second  network  for  further  analysis  and 
graphics  using  a  larger  (Cray)  computer. 

T^e  display  of  /(v„v^,vj)  presents  a  fundamental 
problem  since  four  orthogonal  axes  are  required.  A 
reduced  distribution  /(vh.Vj.)  can  be  displayed  as  a 
surface^  but  this  defeats  the  purpose  of  the  measure¬ 
ment  by  imposing  symmetry.  If  one  notes  that  in 
three  dimensions  “cuts”  through  /(v,  0, 0)  produce 
surfaces  of  constant  /(v),  the  distribution  Inunction 
may  be  visualized  by  inspection  of  a  sequence  of  these 
surfaces.  In  this  case  /(v)  is  measured  on  the  center 
of  the  current  sheet  (x-0,  z-0)  30  /xsec  into  the 
reconnection  event.  In  Fig.  2(a)  the  surface 
C— /(v)//(0)  -2.5  xl0“^  is  a  sphere  indicating  that 
the  bulk  of  the  distribution  is  Maxwellian.  As  C  —  0 
the  configuration  surface  reflects  the  presence  of  faster 
particles.  The  detached  surfaces  in  Fig.  2(b)  denote 
beamlets  of  fast  particles  (an  electron  beam  with  velo¬ 
city  greater  than  a  Maxwellian  main  body  would  give 
rise  to  two  spheres  m  velocity  space).  A  continued 
scan  to  values  of  Cas  lov/  as  10~^  is  possible  since  the 
electron  currents  are  as  large  as  /,aii//iot  ~-0.3  and  the 
digitizers  have  0.4%  accuracy.  (Transformation  into 
velocity  space  entails  a  division  by  v’.) 

As  C  decreases  several  additional  beamlets  appear. 
The  beamlets  appear  to  “grow”  filaments  which  con¬ 
nect  them  to  each  other  and  the  main  bixly.  Inspec¬ 
tion  of  these  surfaces  indicates  that  some  fast  particles 
have  velocities  directed  oblique  to  the  local  magnetic 
field.  The  continuous  development  of  these  structures 
as  well  as  rotations  about  various  axis  in  velocity  space 
is  best  shown  in  a  computer-generated  movie.^  If  the 
probe  is  moved  into  a  region  of  appreciable  transverse 
field  (Ul  SI 5,  z-0)  no  fast  particles  are  observed; 
the  surfaces  are  all  spherical. 

Although  anisotropies  in  velocity  space  have  been 
previously  observed  by  our  group,’  new  techniques  in 
data  reduction  and  the  ability  to  store  and  process  huge 
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FIG.  I.  (a)  Linear  vector  plot  of  Bi  showing  the  neutral-sheet  magnetic  field  topology  for  a  ten-shot  ensemble  average. 
Sketch  shows  location  of  velocity  analyzer,  (b)  Three-dimensional  wave-number  spectrum  at  f  —  1  MHz  for  one  component 
of  the  cross-spectral-function  tensor  The  bottom  and  side  plots  are  wave  spectra  in  the  ky-k,,  and  k,-ky  planes.  The 

markers  with  size  proportional  to  amplitude  are  projections  of  the  peaks  in  k  space.  The  bounding  surfaces  are  the  cold-plasma 
dispersion  curves  for  oblique  whistlers  at  10  G,  n,“4.9x  lO'Vcm^.  The  vertical  line  is  drawn  through  kx~  k,  =  0.  The 
diamond  in  the  center  is  at  k~0.  The  prolific  mode  structure  falls  within  the  limiting  dispersion  surfaces. 


arrays  of  data  allow  us  to  report  these  detailed  struc¬ 
tures  for  the  first  time. 

The  number  of  particles  per  unit  solid  angle  is 
v^f(v)  which  strongly  weights  the  energetic  ones  in 
calculation  of  the  moments 

JJ** Jo  ji^'*>“f^'»,9,^)v^sm9d9d4>dv. 

The  ratio  of  tail  to  main  body  density  n„,j|/nioi  “0.015 
is  calculated  for  a-0.  The  higher  moments  are,  for 
a  —  1,  current  density  y,-0.72,  J, 

--0.18  A/cm^  for  a -2,  mean  energy  Tr(E) 
“  10.5  eV  with  off-axis  components  10~^  times  axial 
ones;  and  for  a -3,  heat  flux  Q,-1.8,  0^“91.1, 
0,  -  -  7.6  W/cm^.  The  local  current  is  predominantly 
in  the  direction  of  the  axial  magnetic  field  and  agrees 
with  other  measurements.  The  temperature  of  10.5 


eV  indicates  that  the  fast  electrons,  V||^4xl0* 
cm/sec,  have  run  away  in  response  to  the  induced 
parallel  electric  field.^  The  axial  heat  flux  is  carried 
mainly  by  the  runawa''s. 

We  believe  that  the  complex  structure  observed  in 
velocity  space  may  be  qualitatively  explained  by  two 
effects.  In  the  center  of  the  current  sheet  the  trajecto¬ 
ry  of  a  single  collisionless  electron  is  a  spiral  as 
B  —  ByQ.  However,  electrons  arriving  at  the  detector 
originate  at  a  variety  of  spatial  locations  and  have  com¬ 
plicated  orbits’  which  may  not  average  out  to  yield  a 
uniform  velocity  distribution  in  the  beam  direction. 
Perhaps  more  interesting  are  wave-particle  interactions 
which  redistribute  electrons  whose  phase  velocities 
match  those  of  unstable  waves.  There  is  a  great  deal 
of  evidence  for  this  in  the  neutral  sheet. 
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FIG.  2.  Surfaces  in  velocity  space  of  /,( Vi.  v,.  Ui) 
-C  —const,  /,(v-0)  -1.0.  Probe  is  located  in  the  center 
of  the  neutral  sheet,  fl^-10  G.  (a)  C -const -0.025. 
The  surface  is  spherical  indicating  that  the  bulk  of  the  parti¬ 
cles  is  Maxwellian.  (b)  C— 0.022.  (c)  C— 0.021.  The  dis¬ 
tortions  and  beamlets  mark  the  presence  of  electrons  in  the 
electric  field,  £^IIB^,  of  the  current  sheet,  (d)  C— 0.019 
viewpoint  changed  to  —v,  direction  to  reveal  filamentary 
structure. 


The  ensemble-averaged  magnetic  field  is  somewhat 
deceiving;  within  the  neutral  sheet  the  magnetic  fields 
exhibit  large  fluctuations  — 1.  The  broad 

frequency  spectrum  peaks  near  the  lower  hybrid  fre¬ 
quency  (/lh  — 200  kHz)  and  extends  to  the  electron- 
cyclotron  frequency.  A  detailed  study  of  the  instability 
growth  was  not  done  since  it  fully  saturates  within 
several  wave  periods.  Extensive  correlation  measure- 


FIG.  3.  Plot  of  parameter  y4(£|i)  (solid  line),  Eq.  (2), 
evaluated  from  the  distribution  function  shown  in  Fig.  2. 
A  >  0.04  is  an  indication  of  sufficient  anisotropy  to  destabi¬ 
lize  parallel  whistler  waves.  The  dotted  curve  indicates  a 
positive  slope  exists  in  / /,  </£i  at  tail  energies  sufficiently 
high  to  destabilize  Langmuir  waves. 


ments  have  been  performed  and  the  cross-spectral- 
function  tensor  evaluated.^  The  measurements  indi¬ 
cate  that  at  fixed  to  the  wave  number  k  depends  on  an¬ 
gle  as  summarized  in  Fig.  I  (b).  Data  from  two  orthog¬ 
onal  planes  (x-y,x-z)  were  used  to  construct  the 
three-dimensional  wave-vector  spectrum.  Shown  is 
the  multimode  structure  of  one  component  of  the 
cross-spectral  tensor  at  a  fixed  frequency  (/ -  I  MHz). 
The  dispersion  surfaces  for  oblique  whistlers,’ 

Acc/<u==o>^/n>(iU(COsP  — cu),  (1) 


fill  a  volume  of  k  space  since  w,,  vary  throughout 
the  reconnection  region.'®  "  No  random  waves  are 
observed  outside  the  surfaces  described  by  Eq.  (1)  for 
the  experimental  extrema  (/V“A9xl0'^  cm“\  B 
- 10  G).  Tests  on  a  number  of  these  modes  have  es¬ 
tablished  that  they  are  right-hand  circularly  polarized. 

The  measured  phase  velocity  of  the  Doppler-shifted 
whistlers  =4.3  xlO*  cm/sec,  ( (kii)  =0.6 

cm~')  falls  into  the  range  of  the  observed  tail.  An  an¬ 
isotropy  test  which  indicates  that  the  measured 
electron-distribution  function  is  prone  to  a  parallel 
whistler  instability' ^  is  that  A  (£„)  >  ( ojf /<u  - 1 )  “ ' 
where 


/<(£,,) 


d//9£i  -  9//9£  II )  dE^ 


(2) 


This  expression  has  been  evaluated  and  /4(£||)  is 
shown  in  Fig.  3.  In  a  range  of  energies  the  instability 
criterion  /4(£||)  >  y4<.  =0.04  is  satisfied.  The  max¬ 
imum  is  at  V||  =4.7x10*  cm/sec,  close  to  the  phase 
velocity  of  parallel  whistlers. 

The  reduced  parallel  distribution,  g(uii) 
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“/ f{E)dEi  [denominator  of  Eq.  (2)1,  is  also  shown 
in  Fig.  3.  It  exhibits  a  positive  slope  which  satisfies 
the  criterion  for  Cherenkov  instabilities 
(oblique  whistlers,  Langmuir  waves,  etc  ).  Measure¬ 
ments  of  the  microwave  emission  with  horn  antennas* 
have  shown  a  strong  peak  at  the  plasma  frequency 
f  —  fp—\2  GHz,  indicative  of  beam-driven  strong 
Langmuir  turbulence.’’ 

A  more  general  expression  to  test  for  growth  of 
oblique  whistlers”  has  been  used  with  measurements 
of  two-dimensional  distribution  functions  in  the 
Earth’s  bow  shock.’'*  Although  whistler  noise  is  ob¬ 
served  in  the  shock  transition  region  it  is  not  possible 
to  measure  wave-number  spectra  and  wave  dispersion 
in  space,  and  distribution  measurements  are  insensi¬ 
tive  at  low  energies  within  the  main  body.  A  number 
of  interesting  computer  simulations  along  these  lines 
have  been  carried  out”  ’''  and  have  shown  the  impor¬ 
tance  of  trapped  particles  in  the  deformation  of  /(v). 
Roux  and  Pellat”  have  shown  analytically  that  a  beam 
or  hole  is  formed  in  the  distribution  function  at  the 
whistler  resonant  velocity.  We  are  presently  extending 
our  analysis  to  examine  the  growth  of  the  observed 
oblique  whistler  modes. 

In  this  work  we  have  pointed  out  that  the  structure 
of  the  electron  distribution  function  in  a  turbulent, 
current-carrying  magnetoplasma  is  rich  in  detail.  It  is 
long  known  that  particles  resonant  with  waves  are  scat¬ 
tered  in  velocity  space  as  they  interact  with  them.  We 
have  observed  this  to  occur  in  three  dimensions  when 
oblique  whistler-wave  turbulence  is  present.  Since  fas¬ 
ter  particles  heavily  weight  the  moment  calculations, 
the  presence  of  low-density  tails  can  greatly  affect 
transport  processes  such  as  heat  flux.  We  believe  that 
the  complex  anisotropies  discovered  in  velocity  space 
are  of  great  consequence  in  diffusion  and  transport. 
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ABSTRACT 

The  role  of  laboratory  experiments  to  the  understanding  of  current 
systems  in  space  plasmas  is  reviewed.  It  is  shown  that  laboratory 
plasmas  are  uniquely  suited  to  make  detailed  investigations  of  basic 
physical  processes  in  current-carrying  plasmas.  Examples  are  given  for 
double  layers,  current-driven  instabilities,  and  the  plasma  dynamics 
at  magnetic  neutral  points  during  reconnection.  Observations  of 
current  sheet  disruptions  show  the  coupling  between  local  plasma 
phenomena  (double  layers)  and  global  circuit  properties  (magnetic 
energy  storage). 


INTRODUCTION 

Although  laboratory  plasmas  and  astrophysical  plasmas  consider 
a  vastly  different  parameter  regime, unstable  current  systems  are 
encountered  in  both  cases.  Although  the  observational  methods  are 
opposite,  local  versus  remote,  there  is  evidence  for  common  physical 
processes:  instabilities,  plasma  energization,  reconnection.  The 
investigation  of  these  processes  from  various  viewpoints  contributes 
to  the  general  understanding  even  though  a  one-to-one  correspondence 
is  often  not  established. 

Carefully  designed  laboratory  experiments  establish  stable  plasmas 
with  simple,  known  boundary  conditions,  and  a  variety  of  diagnostic 
tools  for  in-situ  time  and  space  resolved  measurements  of  fields,  waves, 
particle  distributions.  These  features  are  especially  suited  for  ob¬ 
serving  localized  processes  such  as  double  layers  which  have  been  well 
established  in  the  laboratory  (Torven,  1979;  Sato,  1982).  Similarly, 
electrostatic  and  short-scale  length  electromagnetic  turbulence  has  been 
successfully  studied  in  many  laboratory  plasmas  carrying  currents 
(Gekelman  and  Stenzel,  1978;  Hollenstein  and  Guyot,  1983;  Stenzel,  1977; 
Gekelman  and  Stenzel,  1983).  The  investigation  of  magnetic  reconnection 
requires  plasmas  of  larger  dimensions,  on  the  order  of  the  ion 

47 

M.  R.  Kundu  and  G.  D.  Holman  feds.),  Unstable  Current  Systems  and  Plasma  Instabilities  in  Astrophysics.  4  7-60. 

©  1985  by  the  I AU. 


69 


48 


R.  L  STENZEL  AND  W.  GEK  ELMAN 


collisionless  skin  depth  c/iopj.  In  earlier  investigations  (Bajm  and 
Bratenahl ,  1980;  Frank,  1976;  Ohyabu  et  al.,  1974)  this  has  been  esta¬ 
blished  by  high  densities  (ne  =  10^^  cm”^)  and  small  physical  sizes 
(1  ...  10  cm)  which  makes  the  diagnostics  difficult.  Recent  reconnec¬ 
tion  experiments  at  UCLA  (Stenzel  and  Gekelman,  1979;  Gekelman  et  al., 
1982;  Stenzel  et  al.,  1983a)  have  the  opposite  approach;  large  scales 
(1  ...  2m)  and  medium  densities  (n^  =  10^2  cm-^)  so  as  to  optimize  the 
diagnostic  aspect.  In  these  experiments  a  classical  magnetic  neutral 
sheet  has  been  established;  reconnection  and  particle  energization  are 
observed  facts.  The  current  sheet  associated  with  the  magnetic  topol¬ 
ogy  of  a  flattened  X-point  exhibits  microscopic  instabilities  (magnetic 
turbulence)  and  macroscopic  instabilities  (current  disruptions).  The 
latter  case  is  most  interesting  since  reconnection,  current  sheet 
disruptions  and  double  layers  occur  simultaneously.  It  is  an  example 
of  the  complicated  but  probably  realistic  nonlinear  plasma  dynamics  at 
magnetic  null  points  which  also  requires  the  knowledge  of  the  global 
circuit  properties. 


EXPERIMENTAL  CONFIGURATION 

The  UCLA  reconnection  experiments  are  performed  in  a  device  shown 
schematically  in  Fig.  1.  A  linear  discharge  plasma  column  is  generated 


Fig.  1.  Schematic  picture  of  the  experimental  arrangement,  (a)  Cross- 
sectional  view  showing  parallel  plate  electrodes  with  pulsed  currents 
Ig  and  magnetic  field  lines  without  plasma,  (b)  Side  view  of  the 
device  with  main  electrodes,  currents  (Ig,  Ig)  electric  fields 
(t  =  A  -  v<t)p),  and  magnetic  fields  (S  =  Bj_  +  Byo).  The  coordinate 
system  common  in  magnetospheric  physics  has  been  adopted  where  y  is 
along  the  neutral  line  (device  axis),  x  is  along  the  horizontal  neutral 
sheet,  and  z  is  normal  to  the  sheet. 


with  a  1  m  diameter  cathode.  Detailed  plasma  diagnostic  tools  are 
employed  in  conjunction  with  a  state-of-the-art  digital  data  processing 
system.  Time  and  space  resolved  probe  measurements  of  fields  (f,  t) , 
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4)1asiTia  properties  (ng,  T,  (ji,  v)  and  distributions  [f(v,  r,  t)]  are 
performed.  The  analog  data  are  digitized  with  100  MHz,  32K  A-D  con¬ 
verters,  evaluated  on-line  with  an  array  processor  linked 
to  a  Cray  computer. 


(a)  K.(x,z)  t  =  27/Jsec  He  ,  3x10'^  Torn 


X  (cm) 


Fig.  2.  Measured  magnetic  fields  and  flow  field  during  reconnection, 
(a)  Transverse  field  Sj_  (x,  z)  at  a  fixed  time  t,  position  y  =  137  cm, 
showing  the  classical  neutral  sheet  topology,  (b)  Transverse  ion 
flow  velocity,  vj_  (x,  z),  normalized  to  the  local  sound  or  Alfven 
speed,  Cj  =  (kTg/m,j  )i/2  =  c^  =  (B2/2yQnm,j  Vertical  compression 

and  horizontal  jetting  are”in  qualitative  agreement  witn  reconnection 
models. 


The  plasma  of  typical  density  ng  =  10^^  cm"^  and  temperature 
Tp  10  eV  is  uniform,  quiescent,  essentially  collisionless,  and  highly 
reproducible  in  pulses  of  duration  tp  =  5  msec,  repeated  every  t^.  =  2  sec 
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After  preparing  the  plasma  in  a  uniform  axial  magnetic  field  (By  =  lOG) 
a  pulsed  transverse  magnetic  field  =  lOG)  is  applied.  Its 

topology  vacuum  contains  an  X-type  neutral  point  along  the  axis  of 
the  device.  In  the  presence  of  plasma,  electron  currents  are  induced 
which  flow  preferentially  in  regions  of  Bj_  =  0.  They  are  so  large 
(I  =  lOOOA)  as  to  modify  the  field  topology  self-consistently,  forming 
a  Dungey  (1958)  type  neutral  sheet  during  the  interval  of  rising 
applied  fields  ( 0  <  t  s  lOOusec)  shown  in  Fig.  2a.  The  corresponding 
current  sheet  (J  =  v  x  t/ug)  has  a  thickness  (az  =  5cm)  between  the 


ARGON  .  10'^  TORR 


TIME  t  (10  ijsec/div) 


Fig.  3.  Disrupted  center  plate  current  Ia(t)  (top  trace)  and 
instantaneous  plate  voltage  Va(t)  (bottom  trke)  measured  in  situ.  At 
each  disruption  a  voltage  spike  is  generated  owing  to  the  distributed 
circuit  inductance  L.  Note  that  Ldl^/dt  »  =  lOV. 


collisionless  electron  skin  depth  (c/upg  =  0.6  cm)  and  the  ion  Larmor 
radius  (r^i  =  30  cm).  Thus,  we  are  investigating  the  fine  structure 
of  the  diffusion  region  which  is  beyond  the  resolution  of  present 
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•solar  observations  and  has  neither  been  observed  in  the  magnetosphere. 
Reconnection  occurs  in  our  experiment  at  a  rate  indicated  by  the 
axial  electric  field  along  the  separation  (Ey  =  0.5  V/cm)  (Sonnerup, 
1979)  or  the  normalized  inflow  velocity  of  tne  fluid  into  the  field 
reversal  region  (v/v;^  =  0.3)  shown  in  Fig.  2b  (Vasyliunas,  1975).  For 
uniform  plasmas  and  modest  current  densities  (normalized  electron 
drifts  vj/Vg  =  0.1)  the  current  sheet  is  macroscopically  stable  on  a 
time  scale  Tong  compared  with  the  Alfven  transit  time  across  the  sheet 


CURRENT  DISRUPTIONS  AND  DOUBLE  LAYERS 

The  stability  of  the  current  sheet  with  respect  to  increasing 
current  densities  at  the  separator  has  been  investigated.  This  is 
accomplished  by  raising  the  potential  of  the  central  portion  of  the 
end  anode  on  which  the  plasma  tenninates  axially.  When  monitoring 
collected  current  Ig  to  this  center  electrode,  we  find  that  at 
increasing  current  densities  (v^j/Vg  >  0.3)  spontaneous  current 
disruptions  develop  (see  Fig.  3a).  The  cause  for  this  current  switch- 
off  has  been  inferred  from  detailed  diagnostics  of  the  local  plasma 
properties.  During  a  disruption  the  plasma  potential  rises  in  the 
perturbed  current  channel  to  a  value  much  larger  than  the  dc  potential 
(Vdc  =  lOV)  applied  to  the  end  plate  (see  Fig.  3b).  Simultaneously, 
the  plasma  density  decreases.  These  processes  have  a  finite  axial 
extent.  In  particular,  the  plasma  potential  exhibits  an  abrupt  axial 
drop  (ac#)  =  30V  in  A  y  =  5mm  =  lOOAp),  i.e.,  a  double  layer  is  formed 
in  the  region  where  the  current  is  disrupted. 


Fig.  4.  Schematic  diagram  of  the  important  elements  in  the  physical 
model  for  the  disruptive  instability 


The  large  positive  plasma  potential  can  only  be  explained  in 
terms  of  the  circuit  properties.  Fig.  4  shows  schematically  that  the 
extended  current  path  has  a  distributed  inductance  L  which,  in  the 
presence  of  a  current  drop  results  in  an  inductive  voltage  rise  Ldl/dt 
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This  voltage  drives  the  plasma  potential  positive  and  causes  an 
expulsion  of  ions,  hence  a  density  and  current  drop.  The  current 
decrease  in  turn  reinforces  the  inductive  voltage  so  that  an  explosive 
disruption  develops.  The  current  lost  in  the  center  of  the  original 
sheet  is  redirected  to  the  sides. 

At  the  double  layer  particles  are  accelerated  on  expense  of  magne¬ 
tic  field  energy  stored  at  different  locations  in  the  current  system. 

The  kinetic  energy  first  resides  in  particle  beams.  Subsequent  beam- 
plasma  instabilities  transfer  the  directed  energy  into  waves  and  heat. 
For  example,  microwave  emissions  at  the  electron  plasma  frequency 
are  caused  by  electron  beam-plasma  interactions  (Whelan  and  Stenzel , 
1981). 

Many  of  the  observed  phenomena  of  current  sheet  disruptions  may 
apply  to  space  plasmas  as  well.  The  dynamic  current  modifications 
during  substorms  and  solar  flares  accelerate  particles  to  high  energies. 
The  formation  of  field  aligned  potential  structures  is  well  known 
from  auroral  physics  but  it  is  possible  that  nonstationary  double 
layers  arise  from  inductive  voltages  where  the  energy  storage  is 
remote  from  the  region  of  dissipation.  Type  III  solar  radio  bursts 
involve  beam-plasma  instabilities  possibly  analogous  to  the  present 
emission  process. 


TURBULENCE  IN  CURRENT  SHEETS 

Even  in  the  absence  of  current  disruptions  the  macroscopically 
stable  neutral  sheet  exhibits  a  significant  level  of  turbulence. 
Currents  in  collisionless  plasmas  create  various  fluid  and  kinetic 
instabilities  (Das,  1981).  The  first  step  in  the  observation  of 
turbulence  is  the  mode  identification  from  the  frequency  and  wave- 
number  spectrum.  Then,  the  instability  mechanism  has  to  be  isolated 
which  involves  particle  distribution  measurements.  Finally,  the  effect 
of  the  instability  on  the  macroscopic  transport  properties  is  of 
interest.  Some  of  these  questions  have  been  addressed  in  the  present 
experiment. 

Three  characteristic  modes  have  been  identified.  An  enhanced 
level  of  microwave  emissions  in  the  range  of  plasma  frequencies 
(6  s  f  i  12  GHz)  is  observed  during  reconnection  (see  Fig.  5  top). 

These  nonthermal  emissions  are  due  to  electron  plasma  waves  excited 
by  energetic  electrons.  A  second  spectrum  of  waves  below  the  ion 
plasma  frequency  has  been  investigated  (Fig.  5  bottom).  With  two-probe 
cross-correlation  measurements  the  dispersion  u (k)  is  obtained  which 
identifies  the  noise  to  consist  of  ion  acoustic  modes  (lo  =  kc,). 
Finally,  low  frequency  magnetic  fluctuations  above  the  lower  hybrid 
frequency  (fj^j,  =  200  kHz)  are  observed  (Fig.  6a)  and  studied  in  depth. 
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FREQUENCY  t  (MHz)  WAVENUMBER  k  (cm'’) 

Fig.  5.  Turbulence  spectra  in  a  current  sheet.  Top:  Microwave 
emission  near  the  electron  plasma  frequency  is  greatly  enhanced  during 
the  reconnection  pulse.  Bottom:  (a)  Electrostatic  turbulence  spectrum 
below  the  ion  plasma  frequency  (-25  MHz)  sampled  during  reconnection. 

(b)  Dispersion  relation  u(k)  obtained  from  measured  cross-correlation 
functions  (see  insert)  identifies  spectrum  to  consist  of  ion  acoustic 
modes . 


The  nonstationary  magnetic  turbulence  is  analyzed  by  ensemble 
averaging  two-probe  vector  cross-correlations.  Digital  Fourier 
analysis  in  time  and  space  reveals  a  multitude  of  modas  which  fall 
near  the  average  dispersion  surface  for  whistler  modes  (Fig.  6b).  By 
selecting  individual  modes  the  polarization  properties  B(Ar,  At) 
are  found  to  be  right  hand  circular,  confirming  that  the  modes  are 
indeed  whistlers. 


Fig.  6.  Magnetic  turbulence  in  a  current  sheet,  (a)  Frequency  spectrum 
(lower  hybrid  CO jih/2 11  ==  0.1  MHz;  electron  cyclotroii  toce/2ii  -  30  MHz). 

(b)  Three-dimensional  wave  vector  k  space  with  unstable  modes  of 
the  magnetic  turbulence.  The  cross  spectral  function  <Bxi  at  f  =  1 
MHz  has  been  measured  in  two  orthogonal  planes  and  spatially  Fourier 
transformed.  The  data  points  are  scattered  around  the  theoretical 
dispersion  surfaces  of  whistler  waves. 
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Since  whistler  instabilities  can  arise  from  velocity  space 
anisotropies  we  have  investigated  the  electron  velocity  distribution 
function.  These  measurements  are  performed  with  a  novel  velocity 
analyzer  with  good  angular  resolution  (An  =  10"2  sterad)  (Stetizel  et  al., 
1983b).  The  three-dimensional  distribution  f(vx,  Vy,  v,)  is 
displayed  in  Fig.  7a  as  surfaces  of  constant  value  T,  which  are  nested 
with  the  maximum  (f^ax  “  lO"^**  cm“^  sec^)  near  the  origin  (v  =  0). 
Anisotropies  cause  deviations  from  spherical  surfaces,  clearly 
visible  for  the  contour  f  =  6  x  10"*®  cm"®  sec^  displayed  separately 
in  Fig.  7b. 


Electron  Distribution  logF(v) 


Fig.  8.  Comparison  of  ensemble  average  distribution  <f(v)>  and  root- 
mean-square  fluctuations  f(v)(.^s.  Upper  panel  shows  the  mean  logiQ 
<f(V||,  V|)>  and  two  one-dimensional  cuts,  logio  0)>  and 

logiQ  <fto,  Vj_)>.  Middle  plane  shows  contours  of  constant  <f(V||,  v^)>. 
3  dB  per  contour.  Lower  plane  shows  fluctuations  f^^s  (V||,  Vj_), 

'4  X  10"*®cm"®  sec"®^  per  contour.  Note  anisotropic  distribution  of 
fluctuations  with  approx,  maxima  where  f(v)  exhibits  strong  gradients. 


The  elongation  along  the  magnetic  field  is  produced  by  energetic 
electrons  which  runaway  in  the  electric  field  along  the  separator. 
Runaway  electrons  are  observed  only  inside  the  neutral  sheet.  At  the 
sheet  edge  the  normal  field  component  prevents  the  free  acceleration 
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•of  electrons. 

Kinetic  instabilities  can  be  predicted  both  from  the  shape  of  the 
average  distribution  function  (e.g.  3f/3v  >  0)  and  by  fluctuations  in 
the  distribution  involving  wave-particle  resonances.  From  a  large 
ensemble  of  repeated  measurements  we  have  determined  both  the  mean 
value  and  the  rms  fluctuations  which  are  displayed  in  Fig.  8  in  a 
two-dimensional  velocity  space.  Comparing  the  contour  plots  of  f^ean 
{v|i,  Vj_)  and  f^ms  (V||.  vj_)  we  find  enhanced  fluctuations  near  the 
anisotropic  runaway  electrons.  The  particle  velocities  correspond 
to  the  Doppler  shifted  phase  velocities  of  unstable  whistlers,  and  their 
fluctuation  spectra  overlap.  Thus,  the  origin  of  the  whistler  wave 
turbulence  appears  to  lie  in  the  observed  nonequilibrium  distri¬ 
butions  in  the  current  sheet. 

The  complete  knowledge  of  the  distribution  function  f(v,  r,  t) 
provides  us  with  the  space-time  resolved  fluid  properties  of  which 
we  have  calculated  the  particle  density,  current  density,  mean 
energy  and  heat  flow.  Transport  parameters  have  been  studied  with  the 
help  of  an  array  probe  which  provides  the  bulk  electron  temperature 
with  high  time  resolution  (Wild,  et  al.,  1983).  Correlation  analysis 
of  temperature  fluctuations  in  the  current  sheet  indicate  that  heating 
occurs  in  bursts  of  short  durations  (at  =  3)jsec)  generated  in  narrow 
regions  (s  3  cm  diam. )  propagating  and  diffusing  across  the  field  at 
approximately  the  Alfven  speed.  Heat  bursts  correlate  with  current 
pulses.  The  cross-field  heat  conductivity  is  larger  than  the  classical 
value.  The  microscopic  heating  and  transport  processes  in  collisionless 
plasmas  undoubtedly  involve  plasma  turbulence  effects,  but  further 
investigations  are  required  to  understand  the  processes  in  detail. 
Laboratory  plasmas  are  highly  suited  for  this  task. 
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DISCUSSION 

Sonnerup:  What  is  the  Ion  gyroradius?  What  Is  the  reconnection 
rate? 

Stenzel:  The  calculated  ion  Larmor  radius  is  r^j  ^  30  cm  for  Argon, 
Tel  ^  10  cm  for  Helium  at  <B>  2^  20  G.  The  reconnection  rate,  as 

determined  by  the  observed  upstream  inflow  velocity,  is  given  by 

M  =  —  -v  0.2. 

Van  Hoven:  What  is  the  driving  electric  field  strength 
when  you  show  disruptions?  Is  there  any  E-field  threshold  for  this 
behavior? 

Stenzel:  We  are  in  the  regime  where  the  electron  drift  velocity 
approaches  the  thermal  velocity  or  E/Ep^gj^^gj.  2i  1.  There  are  no 
disruptions  observed  for  vj/vg  v  0.1;  they  build  up  in  the  range 

0-2  ^  vj/Vg  <  1. 

Benford:  Your  6  GHz  radiation  seems  to  be  evidence  for  strong 
electromagnetic  emission  near  uip.  What  is  the  efficiency  of  this 
radiation,  and  how  does  it  compare,  roughly,  to  microwave  emission  from 
flaring  loops  and  type  III  bursts? 

Stenzel:  Microwave  signals  at  f  v  fp  2l  ^  GHz  are  observed  both  as 
electrostatic  signals  in  the  beam— plasma  region  and  as  electromagnetic 
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^^aves  outside  the  plasma.  The  mode  conversion  on  density  gradients  and 
scattering  off  ion  acoustic  modes  has  a  rather  low  efficiency,  estimated 
to  be  Pem/^es  10”**.  I  have  pointed  out  the  qualitative  analogy  between 
type  III  solar  radio  bursts  and  our  (Op  emission  which  occur  during 
impulsive  reconnection  events. 

Sato:  Can  you  distinguish  between  double-layer  acceleration  and  the 
acceleration  by  reconnection  (slow  shocks,  for  example)? 

Stenzel:  Yes,  the  acceleration  by  the  double  layer  is  mainly  in  the 
direction  of  the  separator  (y-axls)  while  the  acceleration  due  to  the 
reconnection  is  in  the  direction  of  the  neutral  sheet  (x-axls) . 

Kundu:  With  regard  to  6  GHz  plasma  emission,  I  believe  that  the 
time  structure  of  these  microwave  bursts  is  less  than  10  v^s.  Does  it 
mean  that  in  an  astrophyslcal  object  such  as  the  solar  flare  plasma, 
one  should  be  able  to  produce  such  fine  time  structure  bursts  intrin¬ 
sically,  except  for  the  fact  that  other  effects  may  modify  such  time 
structures? 

Stenzel:  The  time  scale  for  the  current  distruption  and  associated 
microwave  bursts  is  roughly  given  by  the  ion  transit  time  through  the 
perturbed  current  channel,  here  t  5  cm/(5xl0^  cm/sec)  ^  10  nsec.  If 
a  solar  flare  of  similar  plasma  parameters  would  exhibit  a  multitude  of 
unstable  current  filaments  of  comparable  scale  size  one  may  Indeed 
expect  the  emission  to  be  randomly  modulated  at  the  above  mentioned 
time  scale. 

Birn:  In  magnetotail  simulations  the  asymmetry  of  the  configuration 
with  gradients  parallel  to  the  current  sheet  and  a  magnetic  field  com¬ 
ponent  perpendicular  to  it  seems  to  have  an  important  effect.  Could 
you  study  such  effects  with  your  device? 

Stenzel:  Yes,  such  studies  could  be  done.  We  could  establish  such 
asymmetries  by  inclining  the  current-carrying  plates  at  an  angle  with 
respect  to  one  another. 

Dum:  You  pointed  out  the  Importance  of  the  electrodynamics  of  the 
plasma  and  external  circuit  for  Double  Layer  (DL)  formation.  I  think 
this  is  a  very  significant  point  for  our  going  beyond  the  usual  electro¬ 
static  laminar  solution  of  the  one  dimensional  Vlasov  equation  that  is 
offered  as  an  explanation  even  for  such  highly  dynamic  processes  as 
auroras.  Nevertheless,  particles  are  also  important.  What  do  you 
think  are  the  particle  boundary  conditions  essential  to  DL  formation 
such  as  injection  or  reflection  from  electrodes,  etc.  Ir.  the  very 
recent  experiment  of  Guyot  and  Hollenstein  (Phys.  Fluids  1983)  DL 
formation  seems  to  depend  on  whether  one  electrode  is  heated  or  not. 

How  do  you  trigger  DL  in  a  controlled  manner,  as  you  mentioned? 

Stenzel:  The  anode,  which  is  the  only  relevant  boundary,  absorbs 
all  electrons;  thus  we  have  no  control  over  the  injected  particle  dis¬ 
tributions  as  in  Hollensteins  work. 

I  mentioned  that  we  now  trigger  current  disruptions  (not  necessarily 
double  layers)  in  a  controlled  manner.  This  is  done  by  pulsing  a  thin 
slab  of  magnetic  field  (82)  transverse  to  the  current  sheet  (Jy)  which 
causes  a  total  current  disruption.  Large  inductive  voltages  drop  off 
across  the  slab,  possibly  forming  again  double  layers. 

Bratenahl:  Did  I  understand  that  you  are  now  looking  for  spontaneous 
development  of  double  layers;  that  is,  double  layer  development  without 
artificial  stimulation? 
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Stenzel:  The  work  reported  did  describe  spontaneous  development  of 
double  layers  without  artificial  stimulations.  We  are  now  working  on 
controlled  disruptions  of  the  entire  current  sheet. 

D.  Smith:  You  mention  that  when  a  double  layer  forms,  a  3-D  potential 
results  with  radial  ion  acceleration.  Could  you  explain  in  more  detail 
how  this  occurs? 

Stenzel:  The  potential  structure  is  three-dimensional,  since  only  the 
perturbed  central  part  of  the  current  sheet  rises  to  a  large  positive 
potential  over  a  finite  axial  length.  Since  the  ions  are  unmagnetized, 
they  are  accelerated  in  both  radial  and  axial  directions.  This  leads  to 
a  large  ion,  hence  density  loss  from  the  perturbed  current  channel. 

The  pump-out  of  plasma  ultimately  causes  the  current  loss. 

Henriksen :  Does  the  electron  distribution  have  a  characteristic 
shape  in  velocity  space? 

Stenzel:  Yes,  we  have  made  extensive  measurements  of  the  electron 
distribution  function  and  found  asymmetric  distributions  even  without 
current  dusruptions.  The  current  sheet  contains  a  large  number  10%) 
of  runaway  electrons  which  are  accelerated  by  the  inductive  electric 
field  along  the  separator.  These  particles  strongly  influence  the 
transport  processes  and  provide  a  source  of  free  energy  for  driving 
various  microinstabilities.  ' 
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ABSTRACT 

Magnetic  reconnection  has  been  studied  in  a  laboratory  experiment  designed  to  model  the 
basic  two-dimensional  neutral  sheet  configuration.  However,  the  focus  has  been  put  on  the 
inner  region  of  the  neutral  sheet  where  the  ions  are  effectively  unmagnetized  and  MHD 
concepts  are  violated.  In  this  parameter  regime  driven  reconnection  is  governed  by  the 
fast  dynamics  of  electrons.  In  true  neutral  sheets  (B^  -  0)  the  current  is  carried  by 
electrons.  Thin  current  sheets  (Az  i  c/tupg)  rapidly  form  multiple  X  and  0  points  due  to 
the  onset  of  the  collisionless  electron  tearing  mode.  Magnetic  energy  is  transported  along 
the  separator  at  the  speed  of  whistler  waves  rather  than  Alfven  waves.  Due  to  space  charge 
separation  the  reconnection  electric  field  Ey  is,  in  general,  not  constant  along  the 
separator  but  localized  near  boundaries,  nonuniformities  in  density  and  magnetic  fields 
which  limit  the  current  ly.  This  leads  to  localized  particle  acceleration,  formation  of 
anisotropic  velocity  distributions  and  instabilities.  Reconnection  and  energization  can  be 
spatially  separated  which  shows  the  importance  of  investigating  both  the  global  current 
system  as  well  as  critical  local  plasma  properties.  Experiments  of  current  sheet  disrup¬ 
tions  are  performed  which  demonstrate  the  processes  of  magnetic  energy  storage,  transport, 
conversion  and  dissipation.  Double  layers  and  shock  waves  can  be  produced  by  current 
disruptions.  The  laboratory  experiments  show  new  dynamic  features  of  reconnection 
processes  not  considered  in  MHD  models  yet  relevant  to  narrow  current  sheets  or  the  center 
of  thick  sheets. 


INTRODUCTION 


Since  its  introduction  the  concept  of  magnetic  field  line  reconnection  has  been  formulated 
in  terms  of  MHD  theories  /1-A/.  This  was  based  on  the  conception  that  all  gradient  scale 
lengths  in  space  plasmas  are  large  compared  with  the  particle  gyroradll  and  that  the  time 
scales  for  magnetic  field  configurations  to  change  are  slow  compared  to  gyroperiods.  How¬ 
ever,  in  the  vicinity  of  magnetic  neutral  points  where  first  the  ions  become  unmagnetized 
the  basic  MHD  concepts  such  as  frozen-ln  magnetic  fields,  Alfven  wave  propagation,  space 
charge  neutrality  are  all  violated.  Thus,  the  central  region  of  neutral  sheets  is  governed 
by  different  physical  processes  than  the  outer  global  region.  While  some  processes  such  as 
single  particle  orbits  have  been  investigated  /5/  a  self-consistent  treatment  for  the  inner 
region  of  neutral  sheets  does  not  seem  to  exist.  This  region  is  also  difficult  to  observe 
from  satellites  in  space  due  to  the  short  transit  time  and  the  variability  of  the  plasma. 

On  the  other  hand  laboratory  plasmas  are  highly  suited  to  establish  neutral  sheets  less 
than  an  ion  Larmor  orbit  thick.  High  resolution  three-dimensional  measurements  of  magnetic 
fields,  electric  fields,  plasma  properties,  particle  distributions  and  waves  can  and  have 
been  performed  under  controlled  repeatable  conditions  /6/.  The  laboratory  experiments  can 
extend  the  parameter  regime  for  reconnection  to  electron-dominated  thin  current  sheets 
which  would  be  difficult  to  Investigate  in  space.  Such  thin  current  sheets  are  also  likely 
to  exist  in  space  plasmas,  possibly  embedded  in  larger  scale  structures  as  indicated  by 
rapid  intense  magnetic  fluctuations  near  null  regions.  Other  cases  for  thin  current  sheets 
are  auroral  arcs  /?/  and  magnetic  cavities  created  by  barium  cloud  releases  /8/. 


We  observe  that  thin  current  sheets  are  dominated  by  the  dynamics  of  electrons.  Magnetic 
field  perturbations  propagate  at  whistler  wave  speeds  rather  than  Alfv6nic  speeds.  Thin 
current  sheets  are  found  to  be  unstable  to  the  fast  electron  tearing  mode.  Space  charge 
electric  fields  produce  electron  Hall  currents  and  a  magnetic  field  component  By  along  the 
separator.  The  resultant  magnetic  field  configuration  approximates  that  of  flux  ropes 
embedded  in  a  neutral  sheet,  however,  with  scale  lengths  small  compared  to  the  ion 
gyroradius  r^,^  and  inertial  length  c/tOpj^.  An  external  magnetic  field  component  By 
magnetizes  the  electrons  throughout  the  current  sheet  and  permits  electron  currents  to  flow 
in  regions  of  normal  magnetic  fields  B^.  While  the  contribution  of  Che  ions  to  the  plasma 
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currents  is  negligible  the  ions  are  important  with  respect  to  setting  up  space  charge 
fields  to  which  they  respond  self-consistently.  For  quasi-stationary  conditions,  the  ions 
stream  across  the  separatrix  as  predicted  by  fluid  theory  but  the  acceleration  mechanism 
differs  from  fluid  models,  i.e.  the  magnetic  7  x  ^  force  is  exerted  on  the  current-carrying 
electrons  while  the  unmagnetized  slow  ions  are  accelerated  by  space-charge  electric  fields. 

The  properties  of  the  global  current  system  including  boundary  conditions  have  been 
investigated.  It  is  found  that  the  reconnection  electric  field  Ey  is,  in  general,  not 
uniformly  distributed  along  the  separator  but  localized  in  regions  of  low  conductivity, 
e.g.  density  depressions,  normal  magnetic  fields,  sheaths  and  double  layers.  Such 
localized  electric  fields  give  rise  to  localized  particle  accelerations.  The  Injection  of 
fast  particles  into  the  ambient  plasma  results  in  anisotropic  distributions,  velocity-space 
instabilities  and  plasma  turbulence.  The  role  of  these  Instabilities  is  seen  in  the 
dissipation  of  the  free  energy  in  the  distribution  function,  i.e  conversion  of  streaming 
energy  into  heat.  The  dynamics  of  a  current  system  has  been  studied  by  a  controlled 
disruption  of  a  current  sheet.  The  processes  of  energy  storage,  release,  transport, 
conversion  to  particle  kinetic  energy  and  dissipation  into  thermal  energy  have  been 
followed.  The  understanding  of  dynamic  reconnection  processes  requires  knowledge  of  both 
global  and  local  properties  of  the  current  system.  These  can  be  obtained  in  the  laboratory 
but  the  extrapolation  of  the  laboratory  results  to  current  systems  in  space  remains  a 
difficult  task. 

EXPERIMENTAL  SETUP 

As  in  the  magnetosphere  the  reconnection  process  in  the  laboratory  is  driven  by  an  external 
energy  source,  the  solar  wind  in  space,  an  external  current  system  in  the  laboratory. 

Figure  la  shows  schematically  the  experimental  arrangement.  A  pulsed  plasma  column  is 
generated  by  a  dc  discharge  between  aim  diameter  cathode  and  an  anode.  The  plasma  of 
parameters  listed  in  Figure  lb  is  uniform,  quiescent,  essentially  collisionless,  and  highly 
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Fig.  1.  (a)  Schenaclc  diagram  of  the  experimental 
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set-up  in  a  side  view  (left 
Basic  plasma  parameters. 


reproducible  in  pulses  of  duration  tp  ^  3  msec  repeated  every  t^  1.3  sec.  The  plasma 
column  is  immersed  in  a  uniform  axial  magnetic  field  0  <  ByQ  <  100  G.  For  the  driven 
reconnection  experiments  a  time-varying  ^  100  psec;  magnetic  field  (0  <  Bj_  <  20  G) 

is  applied  transverse  to  the  plasma  column.  This  field  is  produced  by  pulsing  primary 
currents  1^  through  two  parallel-plate  electrodes  adjacent  to  but  Insulated  from  the  plasma 
column.  The  two  primary  currents  close  via  the  cylindrical  chamber  wall.  Due  to  the 
symDetry  of  the  current  system  there  is  a  magnetic  neutral  line  along  the  axis  of  the 
device.  Typical  magnetic  field  lines  in  vacuum  are  shown  in  the  cross-sectional  view  of 
the  device  (Figure  la,  right  figure).  Since  the  space  between  the  plates  is  filled  with  a 
highly  conducting  plasma  a  secondary  current  I2  is  Induced  anti-parallel  to  the  primary 
current.  The  plasma  current  is  mainly  carried  by  electrons,  it  flows  preferentially  along 
the  neutral  line  (separator)  and  closes  via  the  outer  chamber  wall.  In  a  linear  device  the 
electrons  are  supplied  by  a  source  (cathode)  and  absorbed  by  a  sink  (anode).  The  plasma 
current  is  so  large  (I2  -  lOOOA)  that  the  vacuum  magnetic  field  topology  of  the  curl-free 
X-point  is  fundamentally  .:hanged  resulting  in  a  neutral  sheet,  magnetic  Islands,  or  tearing 
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mode  configurations.  We  emphasize  the  importance  of  the  current  closure  which  is  usually 
not  discussed  in  two-dimensional  reconnection  models.  For  example,  if  in  the  present 
experiment  the  electron  supply  is  absent  (e.g.  in  the  afterglow)  the  plasma  current  drops 
to  the  low  value  of  the  ion  saturation  current  (li/le  -  (“e  ^i/m^  «  10“^ J  and  the 

reconnection  proceeds  almost  as  in  vacuum,  irrespective  of  the  high  conductivity  of  the 
plasma. 

Detailed  plasma  diagnostic  tools  are  employed  in  conjunction  with  a  high-speed  digital  data 
processing  system.  Vector  magnetic  fields  ^(r,t)  obtained  in  situ  with  a  movable  probe 
at  up  to  3000  spatial  locations  (Ar  »  2  cm)  and  1000  temporal  points  (At  s  100  nsec).  From 
repeated  measurements  statistical  averages  are  (mean,  standard  deviation,  correla¬ 

tions).  Distribution  functions  are  measured  with  a  modified  retarding  potential  analyzer 
which  filters  particles  through  a  passive  microchannel  plate  and  thereby  obtains  high 
directional  sensitivity  (AD/4ir  =  lO"^^).  The  small  detector  (->3  mm  radius)  can  be  moved  in 
real  space,  rotated  at  each  position  through  the  two  orthogonal  spherical  angles  6,  ^  so  as 
to  obtain  from  the  differential  particle  flux  the  three-dimensional  distribution  function 
f(v,6,^)  or  f (v3(,Vy,V2) «  Electron  and  ion  phase  space  measurements  f(^,‘]^,t)  are  time 
resolved  to  within  a  few  microseconds.  Ensemble  averages  yield  fluctuations  in  velocity 
space.  Measurements  of  such  multi-dimensional  functions  produce  large  data  flows  (N  >  10^ 
numbers)  which  can  only  be  handled  by  digital  techniques.  The  analog  traces  are  therefore 
digitized  with  100  MHz,  32  K,  8  bit  A-D  converters,  evaluated  on-line  with  an  array 
processor  and  a  VAX  11/730  computer  linked  to  a  Cray  computer  for  further  analysis. 

Besides  particles  and  fields,  plasma  waves  are  investigated.  Frequency  spectra  (10^  <  f  < 
10^  Hz)  and  wavevectors  ic  are  measured  using  two-point  cross  correlation  techniques. 

In  order  to  facilitate  comparison  between  the  geometries  of  the  laboratory  and  the 
magnetotail  we  have  chosen  the  same  coordinate  system  where  x  is  along  the  neutral  sheet,  z 
is  normal  to  the  sheet,  and  y  is  along  the  neutral  line  or  separator. 

EXPERIMENTAL  RESULTS:  THE  NEUTRAL  SHEET 

During  the  time  of  rising  external  flux  (0  <  t  <  100  usee)  the  magnetic  field  topology  in 
the  x-z  plane  assumes  the  shape  of  a  classical  neutral  sheet.  Figure  2a  shows  a  vector 
field  (B^,  B2)  vs.  (x,z)  which  indicates  the  long  (Lx  -  40  cm),  thin  (Lx  «  3  cm)  region  of 
vanishing  transverse  fields  embedded  in  the  reverse-field  geometry.  Normal  to  the  plane 


Fig.  2.  (a)  Magnetic  field  vector  map  (Bx,  B2)  in  the  x-z  plane  (y  ■  137  cm 

from  cathode),  showing  the  classical  neutral  sheet  topology.  The  field 
configuration  is  quasi  stationary,  i.e.  unchanged  as  long  as  external  flux  is 
supplied.  (Byo  ^  20  G  *  const),  (b)  Magnetic  field  variation  Bx(z)  across  the 
neutral  sheet  yields  the  current  density  profile  Jy(z)  ”  (dB3(/dz)/po*  The 
half  width  A  «  3  cm  is  smaller  than  characteristic  ion  scale  lengths  (c^)p^  - 
60  cm,  12  cm),  but  larger  than  electron  lengths  (c/uipe  ==  0.7  cm,  r^^  - 

0.4  cm).  The  current  is  carried  by  electrons. 
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there  is  a  uniform  field  component  By^  20  G.  Figure  2b  displays  at  x  *  0  the  com¬ 
ponent  vs.  z  and  its  derivative  which  yield  the  current  density  Jy  =  The  half 

width  A  ^  5  cm  of  the  current  sheet  is  seen  to  be  small  compared  to  characteristic  ion 
scale  lengths  (c/oipj^  -  60  cm*  r^^  -  12  cm  in  helium)  but  larger  than  the  corresponding 
electron  scale  lengths  (c/ajpg  =  0.7  cm,  r^.^  =  0.4  cm)  and  the  length  of  the  classical 
diffusion  region  -  0.3  cm).  Such  narrow  current  sheets  are  typical  when  the  current 

is  dominated  by  electrons  [Je/Jl  ^  (mi/me)^^^j  which  is  always  the  case  near  a  true  neutral 
sheet  (B2  *  U).  A  finite  By  component  enables  electrons  to  flow  in  regions  of  finite  B2 
which  broadens  the  electron  current  sheet  (A  >  t^^)*  An  electron  current  sheet  of  width 
Ag  -  c/upg  carries  as  much  current  as  an  ion  current  sheet  of  width  A^  ==  c/cup^,  provided 
the  currents  are  driven  by  the  same  electric  fields  Ey.  Thus,  the  role  of  electrons  and 
the  inner  structure  of  the  neutral  sheet  are  very  important  to  the  reconnection  problem. 


Long  thin  neutral  sheets  are  subject  to  tearing  instabilities  /9/.  We  have  established  an 
experimental  configuration  which  is  particularly  suited  to  study  current  sheet  instabili¬ 
ties  since  it  establishes  known  boundary  and  Initial  conditions.  As  shown  in  Figure  3  a 
fine  wire  grid  (transparency  >  80%)  with  a  long  (Lj^  -  75  cm),  narrow  (1  <  <  4  cm)  slot 

is  inserted  in  the  plasma  in  the  midplane  (y  ^  100  cm).  At  a  desired  time  (f  *  0)  the 
grid  is  biased  to  a  large  negative  potential  (~Vg  -  200  V  »  kT^)  which  constricts  the 
current  sheet  t^  the  width  of  the  slot.  The  subsequent  space-time  evolution  of  the 
magnetic  field  B/„  ,,  ~  ,.\  is  measured  from  which  the  current  density  *  7  x  i/u 

is  calculated. 


Fig.  3.  Experimental  arrangement  for  studying  the  stability  of  thin  current 
sheets.  The  electron  current  is  forced  to  flow  through  a  narrow  (Az  -  1  cm) 
long  (Ax  =  75  cm)  slot  in  a  grid  at  y  =  100  cm  pulsed  to  a  large  negative 
potential  (e|Vg|  >>  kTg),  Tearing  mode  growth  in  the  x-z  plane  and  propagation 
effects  in  the  y-dlrectlon  are  observed. 


Figure  4  displays  the  magnetic  topology  B^^  which  has  developed  multiple  X  and  O-points 
within  At  8  ^isec  after  constricting  the  current  to  flow  through  a  1  cm  x  75  cm  slot.  The 
time  and  spatial  scales  Indicate  that  the  spontaneous  tearing  process  is  due  to  the 
collisionless  electron  tearing  mode  /lO/  whose  growth  rate  is  predicted  to  be 


y 


VtT  ®oy 


(1 


(1) 


where  ^he  wavenumber  ot  the  modes  (kjj  :  2x/20  cm“^),  Vg  the  electron  thermal  speed  (Vg 

1.5  X  10  cm/sec),  Bqx  Ihe  asymptotic  component  away  from  the  current  layer  (Bg^  == 

8G).  Boy  Is  tl'e  field  component  along  the  separator  (Bgy  =  20  G),  is  the  electron 
inertial  length  (Ig  -  ~  1  is  the  halfwidth  of  the  current  layer  (i  = 

2  cm).  The  predicted  growth  time  for  the  experimental  conditions  is  =  1*3  psec, 
consistent  with  the  observed  fast  island  growth. 

Figure  5  demonstrates  the  rapid  axial  propagation  of  the  current  system.  Contours  of 
constant  current  density  Jy  are  displayed  in  the  y-z  plane  (x  =  0)  at  different  times  At 
after  the  constriction  of  the  current  channel  by  a  4  cm  x  75  cm  slot  located  at  Ay  -  0. 

One  constant  amplitude  contour  enhanced  as  a  heavy  solid  contour  represents  approximately 
the  front  of  the  constricted  current  sheet.  It  propagates  along  By-  °  20  G  at  a  speed  v  s 
10^  cm/sec  which  is  much  faster  than  the  Alfven  speed  (vj^  »  7  x  lo^ cm/sec)  and  more  ^ 
closely  given  by  the  group  velocity  of  a  whistler  wave  packet  with  frequency  determined  by 
the  switching  time  (At  =  2  psec  =  f”*). 
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TEARING  INSTABILITY 


I  {  cm  i 


Fig.  4.  Magnetic  field  topology  for  a  thin  current  sheet  (A  =  2  c«  = 

3  c/Upe)  exhibiting  multiple  X  and  o’polnts  as  a  result  of  electron  tearing 
modes.  The  measurement  Is  taken  at  Ay  •  10  cm  from  the  center  grid  and  At  • 
8  psec  after  the  injection  of  the  narrow  current  sheet. 

PROPAGATION  Or  CURRENTS 


Oyly.z)  =  const. 


5040  30  20nsO*OX20KO 

—  ay  (an)  —ay  (cm)  GRID 

Fig.  S.  Propagation  of  a  current  sheet  along  the  separator.  At  At  ■  0  a  narrow 
current  sheet  is  Injected  at  the  slotted  (Az  -  4  cm)  grid  located  at  Ay  -  0. 

The  thick  curve  In  the  contour  plot  of  constant  current  density  Jyf™0,y,z) 
roughly  Indicates  the  front  of  the  current  sheet  which  propagetes^at  Initial 
speeds  Vy  »  Vy  S 
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V  -  =  2  Cl)  u  c  =  3  X  10^  cm/sec  (2) 

8 II  cape 

Thus,  the  common  assumption  that  currents  or  magnetic  field  perturbations  propagate  at  the 
Alfven  speed  does  not  hold  near  magnetic  neutral  regions  where  the  Ions  are  effectively 
unmagnetized.  There,  the  electrons  remain  magnetized,  carry  currents  and  transport 
magnetic  energy  along  the  separator  with  their  characteristic  modes,  i.e.  whistlers.  This 
result  is  very  important  when  we  consider  the  disruption  of  a  current  sheet  and  the  release 
of  stored  magnetic  energy. 

THE  RECONNECTION  ELECTRIC  FIELD  AND  PARTICLE  ENERGIZATION 

Two-dimensional  reconnection  models  assume  that  there  are  no  variations  along  the  neutral 
sheet  (3/3y  ~  0).  Such  models  cannot  address  the  problems  of  current  closure  and  the 
realistic  situation  that  global  current  systems  are.  In  general,  not  uniform  along  the 
entire  current  path.  This  has  important  consequences  on  the  reconnection  electric  field  E^ 
which  In  2-D  models  Is  assumed  constant  along  y  and  is  given  by  the  flux  transfer  rate 
3Ay/3t.  In  general  3-D  current  systems,  Faraday's  law  only  determines  tlM^ll^e  Integral  of 
the  reconnection  electric  field  to  be  given  by  the  magnetic  flux  change, a  E'dl  -  -3i()/3t, 
and  It  depends  entirely  on  the  circuit  properties  how  the  electric  field  Is  distributed 
along  the  current  path.  An  important  result  of  the  laboratory  experiments  Is  the  finding 
that  the  reconnection  electric  field  £„  can  be  highly  localized.  This  fact  Is  significant 
for  the  understanding  of  particle  acceleration  during  reconnection. 

Figure  6  shows  in  schematic  form  our  global  current  system  and  gives  three  examples  of 
localized  electric  fields.  Figure  6a  indicates  that  during  quasi-stationary  reconnection 
experiments  (1^  -  const.)  the  reconnection  voltage  (analogous  to  cross-tail  potential 

difference)  Is  observed  to  drop  off  mainly  (-80Z)  at  the  plasma  boundary  near  the  cathode. 
The  bulk,  of  the  highly  conducting  plasma  shields  Itself  from  the  applied  Inductive  electric 
field  by  setting  up  an  opposing  space  charge  electric  field,  hence  behaves  similar  to  the 
highly  conducting  external  current  closure.  It  is  the  cathode  sheath  which,  In  this 
current  system,  limits  the  total  current  and  has  the  highest  effective  "resistance"  V/I. 

The  consequence  of  this  localization  of  Ey  is  that  the  particle  energization  occurs  mainly 
at  the  plasma  boundary  rather  than  In  the  bulk  of  the  plasma.  Reconnection  and  energiza¬ 
tion  are  spatially  separated,  a  concept  never  considered  In  2D  theories.  It  Is  not 
inconceivable  that  similar  situations  also  arise  in  magnetospheric  current  systems;  for 
example,  flux  changes  due  to  reconnection  In  the  magnetotail  may  set  up  localized  electric 
fields  at  the  boundary  layer  and  In  auroral  regions  such  that  the  energy  dissipation  is  not 
In  the  tall  Itself  but  in  other  parts  of  the  extended  magnetospheric  current  system. 

A  second  situation  where  localized  electric  fields  have  been  observed  is  shown  In  Figure 
6b.  At  large  current  densities  spontaneous  fast  current  disruptions  have  been  observed 
/ll/.  When  the  electron  drift  velocity  approaches  the  thermal  velocity  a  spontaneous 
density  depression  leads  to  a  current  decrease,  a  localized  electric  field  which  In  turn 
enhances  the  density  depression,  and  the  process  runs  away  into  a  major  current  disruption. 
The  Inductive  voltage  associated  with  the  rapid  current/flux  change  (V  =  Ldl/dt  =  dji/dt) 
drops  off  at  the  location  of  the  current  disruption  forming  a  transient  potential  double 
layer  /12/.  Again,  electrons  and  Ions  are  locally  accelerated  gaining  kinetic  energy  on 
expense  of  released  magnetic  energy  stored  in  other  parts  of  the  current  system.  The 
Injection  of  particle  streams  Into  the  ambient  plasma  creates  a  variety  of  instabilities 
which  eventually  convert  the  directed  kinetic  energy  into  thermal  energy. 

A  third  example  of  localized  fields  and  particle  accelerations  Is  depicted  In  Figure  6c. 

In  this  experiment  which  will  be  described  In  more  detail  below,  the  electron  current  sheet 
Is  disrupted  by  a  pulsed  slab  of  normal  magnetic  field  (»  Byg  =  6G)  located  in  the 

center  of  the  device.  The  Inductive  voltage  generated  by  the  loss  of  magnetic  flux  due  to 
rapid  reconnection  drops  off  across  the  slab  of  normal  B2,  the  magnetic  "switch".  In  the 
same  region  the  particles  are  energized. 

Finally,  an  example  of  the  energized  particles  Is  shown  In  Figure  7.  The  electron 
distribution  Is  measured  with  a  directional  velocity  analyzer  /13/  in  the  neutral  sheet 
during  quasi-stationary  reconnection.  Figure  7a  shows  basic  current-voltage  traces  of  Che 
analyzer  which  represent  Che  Integrated  electron  flux  vs.  energy.  One  can  clearly 
distinguish  a  population  of  low  energy  bulk  electrons  and  a  group  of  energetic  electrons 
screaming  away  from  Che  cathode.  The  energy  of  the  tall  electrons  exceeds  the  dc  cathode 
potential  (-40  V)  by  approximately  the  reconnection  voltage  (di;)/dc  -  100  V).  Since  these 
fast  electrons  are  observed  already  near  the  cathode  end  of  Che  device  they  have  been 
accelerated  by  a  localized  electric  field.  With  increasing  distance  from  Che  cathode  the 
tails  are  scattered  and  the  bulk  temperature  increases  although  kT^  ~  10  eV  «  d6/dt. 

Figure  7b  Indicates  Che  electron  distribution  function  in  the  middle  of  the  device, 
displayed  as  a  contour  of  constant  value  f(v)  -  10“^^  cm“^  sec^  in  three-dimensional 
velocity  space  ^  •  (Vx,Vy,V2).  The  energize^  electrons  injected  into  the  neutral  sheet 
topology  are  scattered  in  velocity  space  by  E  +  v  x  B  forces  and  form  complicated  tall 
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LOCALIZED  ELECTRIC  FIELDS 
(a)  Quasi -stationary  Reconnection 


Sheath  Potential 
■AV=‘-d0/dt=jEdl 


(b)  Current  Disruptions  (  ,  n;  —  0) 


Double  Layer  AV  =  L  dl/dt 


(c)  Localized  Magnetic  Fields 

z 


y 


Fig.  6.  Examples  of  localized  reconnection  electric  fields  Ey  which  can  only  be 
explained  by  considering  the  global  current  system  rather  than  a  local  section 
containing  a  two-dimensional  current  sheet,  (a)  In  quasl-statlonary 
reconnection  (0  <  t  <  lOOpsec)  the  Induced  voltage  associated  with  the  flux 
transfer  across  the  separatrlx  mainly  drops  off  at  the  cathode  sheath.  Local 
particle  acceleration  by  Inertial  effects  permits  rapid  reconnection  In  Che 
otherwise  highly  conducting  current  system,  (b)  Spontaneous  fast  current 
disruptions  (AC  -  lOusec)  are  observed  at  large  current  densities  which  lead  to 
a  local  density  depletion  and  the  formation  of  a  double  layer.  At  the  double 
layer  Che  excess  magnetic  energy  is  converted  Into  particle  kinetic  energy. 

(c)  A  controlled  rapid  current  sheet  disruption  (At  -  lOpsec)  is  produced  by  a 
pulsed  normal  magnetic  field  B2(t)  which  prevents  the  free  electron  motion  along 
the  separator.  The  inductive  voltage  drops  off  across  the  region  of  where 
the  excess  magnetic  energy  Is  dissipated. 

distributions  /lA/.  The  free  energy  in  these  distributions  drives  various  Instabilities, 
e.g.  electron  plasma  waves  and  whistler  waves,  some  of  which  have  been  studied  in  detail 
/15/. 

The  laboratory  experiments  have  shown  the  Importance  of  investigating  both  the  global 
current  system  and  local  plasma  properties  in  order  to  understand  the  dynamics  and 
energetics  of  reconnection.  The  formation  of  localized  electric  fields  permits  rapid 
reconnection  to  occur  In  highly  conducting  plasmas  by  the  process  of  Inertial  reconnection. 
During  the  finite  transit  time  through  the  electric  field  structure  the  particles  are 
energized  without  undergoing  collisions.  The  field  structure  provides  dissipation  of 
magnetic  energy  in  a  current  system  without  significant  resistivity.  The  energization  of 
particles  can  be  localized,  sporadic,  and  usually  leads  to  anisotropic  distribution 
functions.  Such  features  have  also  been  observed  in  the  magnetotail  with  high  resolution 
electron  detectors  /16/. 

DISRUPTION  OF  A  CURRENT  SHEET 

Reconnection  attempts  to  explain  the  rapid  conversion  of  magnetic  energy  Into  particle 
kinetic  energy  which  is  observed  to  occur  In  solar  flares  and  magnetic  substorms.  Although 
in  the  laboratory  one  cannot  model  the  complicated  global  current  systems  encountered  in 
space  one  can  establish  simpler  but  related  configurations  and  Investigate  the  basic 
physics  of  magnetic  energy  storage,  release,  transport,  conversion  and  dissipation.  Such 
an  experiment  will  be  described  here. 

Figure  8  summarizes  the  events  observed  when  an  Initially  quasl-statlonary  current  sheet  Is 
suddenly  disrupted  in  order  to  release  the  excess  magnetic  energy  of  the  system.  The 
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(a)  INTEGRATED  ELECTRON  FUUX 


‘f 


Fig.  7.  Anisotropic  electron  distributions  produced  by  localized  electric 
fields,  (a)  Electron  flux  vs.  energy  observed  with  a  directional  particle 
analyzer  pointed  In  both  directions  along  the  separator.  From  the  localized 
field  Ey  near  the  cathode  a  flux  of  energetic  electrons  Is  ejected  which 
acquired  most  of  the  reconnection  energy  (ed4>/dt  100  eV).  (b)  Three- 
dimensional  electron  distribution  function  displayed  as  surface  of  constant 
value  f(Vjj,v  ,Vj)  •  10~^®  cm”®  sec^,  showing  the  distribution  of  tall  electrons 
after  propagation  along  the  neutral  sheet.  The  free  energy  in  the  distribution 
excites  Langmuir  and  whistler  turbulence. 


CURRENT  DISRUPTION 
DUE  TO  LOCALIZED  FIELD 


Cascading  Energy 

Deposition 

Fig.  8.  Schematic  diagram  describing  the  results  of  disrupting  a  current  sheet 
Jy  with  a  rapidly  pulsed  solenoldal  magnetic  field  82* 

disruption  mechanism  described  In  Figure  6c  establishes  reproducible  initial  and  boundary 
conditions  so  that  one  can  follow  the  complicated  events  in  space  and  time.  In  space 
plasmas  analogous  disruptions  nay  arise  from  large  amplitude  magnetic  waves  or  fast  tearing 
nodes.  We  find  that  after  applying  the  switch  magnetic  field  82  the  excess  magnetic  energy 
of  the  neutral  sheet  disappears  on  a  time  scale  of  whistler  wave  propagation  along  the 
separator  (~100  cm/lO^  cm  sec”^  *  10  psec).  The  Initially  laminar  current  sheet  cascades 
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Into  aaall  scale  filamentary  current  loops  which  generate  magnetic  helicity  /17/.  A 
potential  difference  corresponding  to  the  magnetic  flux  change  d4)/dt  drops  off  across  the 
magnetic  switch*  The  global  decrease  of  magnetic  energy  Id^/dt  causes  a  local  increase  in 
particle  energy,  i.c.  ions  are  accelerated  by  Ey  across  Bz  since  r^j^  >  B/VB,  electrons  are 
heated  by  cross-field  instabilities,  and  additional  plasma  is  created  by  ionization  of 
neutrals  at  high  electron  temperatures*  The  rapid  local  energization  causes  the  plasma  to 
expand  primarily  in  the  direction  of  the  ion  flow*  The  front  of  the  mass  flow  steepens 
into  a  large-amplitude  electrostatic  shock  wave*  It  slowly  dissipates  its  energy  by 
interaction  with  the  background  ions* 


The  events  summarized  in  Figure  8  will  now  be  supported  by  direct  observations*  Figure  9a 
displays  the  time  dependence  of  the  plasma  current  as  measured  at  the  anode.  After  it 
has  built  up  to  its  quasi-stationary  value  the  local  normal  magnetic  field  is  switched 
on  which  completely  disrupts  the  external  current  flow.  Figure  9b  shows  the  time 
dependence  of  the  dcsainant  oiagnetlc  field  component  B^  in  the  plasma  and,  for  comparison, 
in  vacuum.  The  formation  of  the  neutral  sheet  stores  excess  magnetic  energy  which  upon 
current  disruption  is  released  within  ac  <  10  ysec*  This  rapid  time  scale  cannot  be  ex¬ 
plained  by  Alfven  wave  propagation  (Ac  ~  100  psec)  or  resistive  dissipation  (At  >  1  msec), 
but  by  whistler  wave  propagation  of  rapid  magnetic  perturbations  such  as  shown  in  Figure  5 


Anode 

Current 


(120A/div.) 


Magnetic 

Field 

Bx 

(2G/div.) 


Time  t 

Fig.  9.  Time  scale  of  current  disruption  and  magnetic  energy  loss,  (a)  Plasma 
current  at  the  anode  which,  after  building  up  to  its  peak  value  is  disrupted  at 
t  =  25  psec  by  the  magnetic  switch,  (b)  Magnetic  field  Bx(t)  in  the  neutral 
sheet  region  (with  plasma)  and  in  vacuum.  Prior  to  the  disruption  excess 
magnetic  energy  is  stored  by  the  plasma  current.  After  the  disruption  it  is 
released  within  At  -  10  psec,  a  time  scale  too  short  for  convection  by  Alfven 
waves  or  resistive  diffusion  but  consistent  with  whistler  wave  propagation 
(Figure  5). 

Figure  10  shows  the  dynamics  of  the^current  system  during  the  disruption.  From  repeated 
measurements  of  the  magnetic  field  at  5000  spatial  locations  (Ar  ~  2  cm)  and  102A 

time  increments  (At  ■  50  nsec)  the  cut  rent  density  t)  “  *  ^^Po  calculated,  field 

lines  are  fitted  through  this  vector  field,  and  displayed  in  the  plasma  volume  between 
switch  and  anode  at  different  times  of  the  disruption  process.  Prior  to  the  disruption 
(t  <•  -0.6  psec)  the  current  flows  as  a  laminar  sheet  in  the  x-y  plane  (z  ^  0)  from  anode  to 
cathode.  As  the  electron  Inflow  at  the  left-hand  x-z  plane  is  inhibited  by  the  growing 
component  the  current  begins  to  circulate  within  the  plasma  volume  forming  temporally 
repeatable  but  spatially  random  patterns  of  small-scale  current  loops  and  filaments.  The 
process  of  cascading  from  large  to  small-scale  strucfrea  mlnrlHes  with  the  decrease  in 
the  strength  of  the  current  or  magnetic  energy,  however,  a  causality  between  the  two 
processes  is  not  easy  to  infer.  Similar  cascading  processes  have  also  been  observed  in 
other  laboratory  devices  and  compared  with  solar  flare  phenomena  /18/. 
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field  lines  7  mftS/u, 


Fig.  10.  Field  lines  for  the  current  density  vector  J/x.y.z  t)  different 
tines  of  the  current  sheet  disruption.  The  initially  lt’-’-6.6  usee)  laminar 
current  sheet  cascades  during  the  disruption  into  small  scale  current  loops  and 
filaments.  Simultaneously,  the  magnetic  energy  is  lost.  Due  to  the  controlled 
disruption  the  spatially  random  current  pattern  is  highly  repeatable  from  pulse 
to  pulse. 

From  the  knowledge  of  the  current  density  J  the  vector  potential  '  4^  j  ^(r,t)/ 

,r;;j^r  I  dV  has  been  calculated.  Since  is  also  known,  the  magnetic  helicity  K(c)  - 

A'?  dV  of  the  plasma  current  system  can’he  found.  Helicity  is  a  measure  for  the  linkage 
of  magnetic  flux  of  curtent  systems  and,  in  Ideal  HHD,  is  a  conserved  quantity  like 
magnetic  energy  /19/.  Figure  11  shows  the  observed  change  of  helicity  which  indicates  flux 
linkage  to  occur  during  the  current  sheet  decay.  Helicity  and  energy  decay  on  the  fast 
time  scale  of  electron  whistlers. 


MAGNETIC  HELICITY  K=/A-SdV 


-10  0  10  20  30 
TIME  it  (fjsrc) 


Fig.  11.  Temporal  variation  of  the  magnetic  helicity  of  the  plasma  current 
system.  The  current  sheet  disruption  changes  helicity  and  energy  on  a  time 
scale  of  whistler  wave  propagation  through  the  device. 

In  order  to  account  for  the  loss  of  magnetic  energy  the  plasma  parameters  have  been 
measured  throughout  the  plasma  volume.  The  only  region  where  significant  energization 
occurs  is  around  the  magnetic  switch.  Figure  12  shows  an  axial  plot  of  density, 
temperature  and  plasms  potential  during  the  early  disruption  (At  -  2  psec).  A  potential 
d^p  s  30  V  drops  off  across  the  region  of  normal  B.  which  locally  raises  the  electron 
temperature  to  kTg  “  10  eV,  and  increases  the  density  by  ionization.  With  increasing  time 
the  density  peak  grows  (Sn/n  >  lOOZ)  and  propagates  into  the  low  potential  region  due  to 
ion  acceleration. 
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Plasma  Parameters 
During  Current  Disrufjtions . 


Fig.  12.  Plasma  properties  n^,  Tg,  '(p  vs.  y  during  the  early  phase  of  the 
current  disruption.  A  potential  difference  A(jip  -  30V  drops  off  across  the 
region  of  normal  magnetic  field  which  leads  to  Ion  acceleration,  electron 
heating  and  a  density  Increase  by  Ionization,  all  on  expense  of  the  magnetic 
energy  of  the  disrupted  current  system. 


The  ejection  of  energized  plasma  Is  shown  In  Figure  13.  It  displays  the  Ion  density  vs. 
time  observed  at  Increasing  distances  Ay  from  the  magnetic  switch  toward  the  cathode.  The 
front  of  the  ejected  plasma  Is  observed  to  steepen  Into  an  electrostatic  shock  wave  of 
width  S  =  0.6  mm  =  15  Xp  propagating  with  sonic  Mach  number  M  -  Vy/Cg  =  1.5  along  By^. 

With  Increasing  distance  the  shock  amplitude  decays,  l.e.  wave  energy  Is  transf erred^to  the 
background  plasma.  Such  dissipation  processes  Involve  the  acceleration  and  reflection  of 
background  Ions.  Measurements  of  the  Ion  distribution  function  have  shown  the  tall 
formation  /20/.  At  the  foot  of  the  shock  enhanced  density  fluctuations  are  observed 
associated  with  reflected  Ions  exciting  beam-plasma  Instabilities. 

SUMMARY 


Magnetic  reconnection  has  been  studied  In  a  laboratory  plasma  with  electron  dominated 
current  sheets.  Reconnection  Is  viewed  as  a  property  of  an  entire  current  system  rather 
than  a  local  process.  Investigation  of  the  third  dimension,  l.e.  the  direction  along  the 
current  flow,  has  shown  that  reconnection  (flux  transfer)  and  dissipation  (acceleration  by 
localized  electric  fields)  can  be  spatially  separated.  The  energy  is  transported  along  the 
current  system  by  characteristic  electromagnetic  modes,  which  are  whistlers  for  electron 
dominated  current  sheets.  The  central  region  of  neutral  sheets  (As  l  r^i)  i»  controlled  by 
the  dynamics  of  electrons.  Examples  of  dynamic  current  systems  have  been  presented 
(current  constrictions,  partial  and  total  disruptions)  which  differ  In  the  details  of  the 
events  but  have  in  common  the  rapid  conversion  of  magnetic  energy  Into  particle  kinetic 
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Fig.  13.  The  local  energization  of  Che  plasma  at  the  magnetic  switch  leads  to  a 
mass  outflow  which  steepens  Into  an  electrostatic  shock  wave  propagating  along 
Byg  =>  6  G  toward  the  cathode,  (a)  Ion  flux  vs.  time  at  different  positions  Ay 
snowing  shock  formation  by  nonlinear  steepening,  (b)  Shock  wave  propagation  at 
sonic  Mach  number  M  •  Vg/Cg  =  1.5  (Vg  =  5.5  x  10^  cm/sec,  kTg  -  6  eV).  (c)  Shock 
wave  amplitude  in/n  »  100%  near  region  of  formation.  For  comparison  the 
waveform  of  the  switch  Is  shown,  (d)  Shock  thickness  Ay  •  Vg  At  =  0.6  nn  = 

15Xjj  of  a  fully  steepened  shock  wave  observed  on  two  time  scales  (1  ;jsec/dlv 
top,  0.2  usec/dlv  bottom).  Note  oscillations  at  the  shock  foot  due  to  reflected 
Ions  forming  beam-plasma  instabilities. 

energy  In  highly  conducting  plasmas.  Since  current  systems  In  space  also  differ  greatly 
there  may  be  many  unexpected  forms  and  signatures  of  reconnection  in  addition  to  the 
standard  fluid  motion  across  the  separatrix.  The  laboratory  experiments  suggest  to  look  at 
processes  along  the  separator,  l.e.  electric  fields  at  nonuniformities,  the  role  of 
electrons  In  transporting  energy  and  momentum,  the  role  of  a  By  component,  the  current 
closure,  and  correlations  between  magnetic  flux  variations  and  particle  energization  at 
different  parts  of  the  current  system.  Simultaneous  observations  from  several  spacecraft 
would  seem  a  prerequisite  for  Investigating  three-dimensional  dynamic  current  system. 
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Anomalous  Currents  to  an  Electrode  in  a  Magnetoplasma 


J.  M.  Llrrutia  and  R.  L.  Stenzel 
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(Received  31  March  1986) 

An  electrode  pulsed  to  a  large  positive  potential  ( V  »  kTje)  in  a  collisionless  magnetoplasma 
is  observed  to  draw  currents  far  in  excess  (lOx)  of  the  field-aligned  electron  saturation  current 
prescribed  by  probe  theory.  The  large  currents  are  inferred  to  be  due  to  anomalous  cross-Held 
transport  caused  by  current-driven  instabilities.  Electric  Helds  due  to  space-charge  separation  and 
anomalous  conductivity  lead  in  time  to  ion  expulsion  from  the  current  channel  and  hence  to  densi¬ 
ty  depletion  and  current  collapse. 

PACS  numbers;  52.40.Hf,  52.25.Fi,  52.35.Qz,  52.70.Ds 


The  current  to  an  electrode  biased  above  the  plasma 
potential  in  a  collisionless  plasma  is  a  long-standing 
problem  of  fundamental  interest  in  plasma  diagnos¬ 
tics'"^  and  current  systems  in  space.^  Earlier  reports^ 
have  derived  the  plasma  parameters  from  probe 
characteristics  but  have  not  considered  the  perturbing 
effects  of  saturation  currents  on  the  plasma.  In  this 
Letter  we  report  measurements  of  the  spatial  and  tem¬ 
poral  evolution  of  the  electron  saturation  current  sys¬ 
tem,  J(r,(),  within  the  plasma  as  well  as  detailed  diag¬ 
nostics  of  electric  field,  density,  and  temperature. 
From  the  observations  we  conclude  that  electrons 
drawn  by  a  highly  biased  (  V  »  kT^/e),  pulsed  probe 
raise  the  plasma  potential  in  the  flux  tube  roughly  sub¬ 
tended  by  the  electrode  giving  rise  to  electron  Hall 
currents,  expulsion  of  the  unmagnetized  ions  from  the 
flux  tube,  and  anomalous  radial  electron  currents. 
The  enhanced  parallel  current  near  the  probe  creates 
an  anomalously  low  parallel  conductivity  which 

results  in  resistive  electric  fields,  £||  “y||/cr(y)*|,  far 
outside  of  the  sheath  region.  Ions  are  accelerated  away 
from  the  probe  and  electrons  toward  it  resulting  in  a 
net  density  depletion,  —  0,  and  hence  the 

probe  current  collapses.  Recovery  of  the  current- 
driven  instability  causes  plasma  to  flow  back  to  the 
perturbed  region  and  the  process  repeats  periodically 
creating  pulsating  currents  in  steady  state.  The  large 
transient  current  pulses  penetrate  into  the  plasma 
along  Bo  at  near  electron  thermal  speeds,  but  since 
they  exist  only  over  durations  of  ion  transit  times 
across  the  current  channel  the  penetration  along  Bo  is 
limited  to  distances 


L 


kT,mi 


1/2 

'■/«- 


Where  is  the  radial  potential  drop  in  the  channel 
and  Tf,  the  flux  tube  radius.  Our  plasma  is  much  larger 
than  these  scales;  hence  the  results  are  not  boundary 
dependent  and  can  be  scaled  to  relevant  applications.^ 
The  experiments  are  conducted  in  the  quiescent  af¬ 
terglow  (r-I20  /tsec)  of  the  pulsed  (trep“2  sec, 
ton”4  msec)  plasma  (n,  =  2x10"  cm~’,  T,  —  ST, 


=  1.5  eV,  £o=30  G,  ;;„=:4xl0-''  Torr,  Ar) 
schematically  shown  in  Fig.  Ua).  The  electrode  is  a 
tantalum  disc  of  0.8-cm  radius  (r^)  insulated  on  one 
side.  It  is  biased  to  F  — 80  V  {»kTfle)  via  a 
transistor  switch  with  respect  to  the  grid  anode  and 
chamber  walls.  The  magnetic  field,  B(r,/),  induced  by 
the  current  is  measured  with  a  probe  consisting  of 
three  electrostatically  shielded,  small  (1-cm  diam), 
orthogonal  magnetic  loops.  The  probe  tip  is  capable  o.‘' 
scanning  a  suitably  large  volume  (0  <  /•  <  8  cm, 
0  <  z  <  55  cm)  about  the  electrode.  The  current- 
density  vector  field,  J(r./),  is  calculated  via 
J-VxB/fio-  The  local  plasma  parameters  («,,  kT,, 
<iipi)  are  obtained  with  a  small  ( 1  x  l  .5-mm^)  Langmuir 
probe.  The  probe  tip  can  also  scan  a  volume  about  the 
electrode.  The  probe  bias  is  changed  in  successive 
shots  and  the  probe  current-voltage  characteristic  as  a 
function  of  time  is  reconstructed  at  each  spatial  point. 
The  acquisition  and  handling  of  the  necessary  large 
data  set  is  done  by  computer.^ 
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FIG.  1.  (a)  Experimental  device,  (b)  Electrode  current  / 
vs  time  f  at  K  -  80  V. 
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The  normally  Hoating  electrode  is  pulsed  (rise 
time  —  200  nsec,  /on  ^  ^0  fisec)  to  the  supply  voltage 
The  probe  current  exhibits  a  temporary  increase  above 
the  expected  electron  saturation  current  ( ,,,  >  wrp 
^  ne[kTf/2vmfV^^  —  I  A)  whenever  the  supply  vol¬ 
tage  exceeds  the  plasma  potential.  Current  overshoots 
have  earlier  been  observed^  and  explained  theoretical¬ 
ly^-^  by  excess  ions  in  the  probe  sheath  but  here  we 
will  show  that  the  underlying  physics  is  different. 
After  the  first  current  overshoot  smaller  current  peaks 
appear  at  somewhat  regular  time  intervals.  This  im¬ 
plies  that  no  steady-state  (d/d/  =  0)  currents  exist, 
presumably  because  the  plasma  has  been  radically  per¬ 
turbed  by  the  first  overshoot.  This  Letter  concentrates 
on  the  first  overshoot  (/m»x=  10  A)  caused  by  an  80-V 
bias  [Fig.  Kb)]. 

The  penetration  of  currents  or  magnetic  fields  is  ad¬ 
dressed  by  observation  of  the  temporal  evolution  of 
magnetic  energy  density,  as  shown  in  Fig.  2. 

The  front  of  the  magnetic  energy-density  perturbation, 
arbitrarily  chosen  as  the  first  contour,  assumes  a  speed 
of  V  —  2  X  10^  cm/sec,  nearly  the  thermal  electron  drift 
speed  to  the  probe,  -  (A:7i/2ir/«,  — 2xi0’ 

cm/sec,  but  lower  than  the  speed  of  a  whistler  wave 
packet  of  2-/isec  duration, 

Vf  —  2c/<o^(«u<o„)'^^  —  9x  lO’  cm/ sec. 


FIG.  2.  Contours  of  constant  magnetic  energy  density  in 
the  r-2  plane  at  different  times.  Contour  separation  is  0.2 
nJ/cm^.  The  front  asymptotically  assumes  a  speed  of  2  x  10’ 
cm/sec.  The  electrode  is  shown  to  scale  in  the  first  panel; 
however,  it  is  located  at  z  —  0  cm. 
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and  much  faster  than  the  Alfven  spzed, 
v A  “  Bq/  ( fionm, )  =  3  X  1 0‘  cm/ sec, 

and  the  ion  sound  speed, 

Tf! m, )  =  2  X  1 0*  cm/ sec. 

The  region  of  maximum  energy,  which  is  also  the  re¬ 
gion  of  maximum  current  density,  does  not  extend  far 
(z  <  20  cm)  from  the  electrode.  The  axial  divergence 
implies  that  currents  are  not  confined  to  flow  along  the 
flux  tube  but  must  funnel  radially  in  from  the  entire 
plasma  into  the  electrode. 

Examination  of  the  current  density  reveals  the  ex¬ 
istence  of  both  parallel  and  perpendicular  currents. 
Figures  3(a)  and  3(b)  display,  for  example,  the  current 
density  (y||<0.7  A/rm^,  yj.<0.4  A/cm^)  in  the 
plane  z  —  5  cm  at  /  -  1.5  #isec.  The  direction  of  7^^  is 
mainly  azimuthal  to  the  flux  tube,  i.e.,  a  Hall  current, 
and  partly  radial  outside  the  flux  tube.  Consequently, 


FIG.  3.  (a)  Topological  map  of  Ji  for  a  limited  region 
about  the  origin  for  the  plane  z  -  5  cm  at  /  - 1.5  sec.  Con¬ 
tour  spacing  is  0.08  A/cm^  Maximum  value  is  0.8  A/cm’. 

(b)  Same  as  (a)  but  for  Jf  Maximum  value  is  0.4  A/cm’. 

(c)  Tnqectory  of  a  current  density  line  at  /  -  I.S  fisec  exhi¬ 
biting  the  right-handed,  conical  helix  of  varying  pitch  path. 

(d)  Current  at  the  boundaries  of  the  explored  volume  as  a 
function  of  time. 
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the  current  flows  along  nested  right-handed  conical 
helices  whose  apex  is  at  the  electrode.  Because  J  JJn 
is  a  function  of  space,  the  pitch  and  conical  angle  of 
the  helices  vary  continuously.  A  typical  example  out 
of  the  infinite  set  of  current-density  lines  is  shown  in 
Fig.  3(c).  While  the  current  density  describes  the  local 
properties  of  the  current  system  one  would  also  like  to 
know  the  integrated  currents  /  -  /  J  •  dK.  Figure  3(d) 
shows  that  the  current  across  the  surface  closest  to  the 
electrode  (/|i.i-scm)  evolves  in  time  as  the  externally 
measured  current  [Fig.  Kb) I.  Its  magnitude,  howev¬ 
er,  is  not  the  same  since  it  is  measured  at  r  -  S  cm. 
Current  flowing  through  the  opposite  surface 
(/„  ,_55j„)  is  initially  much  smaller  than  the  current 
at  z  —  5  cm.  However,  the  perpendicular  current  flow¬ 
ing  out  of  the  four  side  surfaces  of  the  rectangular 
volume  explored  [Fig.  3(c)l  account  for 
spite  of  the  low  magnitude  of  7^  boundaries. 

Hence,  current  is  funneled  into  the  electrode.  Only 
after  the  current  overshoot  begins  to  decay  does 
Im.z-55ciii  overtake  as  the  source  for  /||.,_scm- 
The  late  appearance  of  flux  across  the  plane  at  z  —  SS 
cm  is  clearly  due  to  the  propagation  of  the  current 
front.  The  closure  path  for  the  current  therefore 
evolves  in  time  and  probably  encompasses  the  entire 
plasma  volume.  It  is  appropriate  to  note  that 


FIG.  4.  Vector  field  map  of  the  net  electric  field, 
E(r,/)-  -  V,/*,  for  selected  times.  The  electrode  is 

shown  to  scale  in  the  bottom  panel;  however,  it  is  at  z»0 
cm. 


2  /  -  ^  J  •  dA  is  not  strictly  zero  at  all  times  ( =  5%  of 
/x  at  /-l.S  /isec)  indicating  that  current  measure¬ 
ments  may  have  a  similar  error. 

Upon  switchon  of  the  electrode  bias,  the  plasma  po¬ 
tential  is  observed  to  rise  with  negligible  delay 
throughout  the  entire  plasma  by  as  much  as  30 
V  and  near  the  probe  by  ^  40  V.  The  global  po¬ 
tential  shift  is  caused  by  the  return  electrode  (grid 
anode  and  chamber  walls)  which  must  draw  an  ion 
current  equal  to  the  electrode  electron  current  so  as  to 
maintain  charge  neutrality.  The  local  potential  gra¬ 
dients  near  the  electrode  are  observed  to  extend  well 
beyond  the  sheath  region  [thickness  6  =  (Ai^/ 
kr,)'^^Xfle==5X/v  — 015  mm].  Figure  4  shows  vec¬ 
tor  maps  of  the  net  electric  field,  E(r,r) 
-  -  V^- V(n^r,)//tf.  Inductive  electric  fields  are 
n^ligible  compared  to  the  electrostatic  fields  [8^4/ 
8/  a:  1  mV/cm].  One  can  observe  the  following 
features:  (a)  At  early  times  (r  —  1  /isec)  the  electric 
field  is  predominantly  radial  {E,^2  V/cm)  as  a  result 
of  an  excess  positive  charge  [(»/-«,)/«,  10“*]  in 

the  flux  tube  caused  by  electron  drain  to  the  probe; 
and  (b)  as  the  probe  current  peaks,  an  unexpectedly 
strong  parallel  field  [£,  ^4  V/cm  »  V(«kr,)/ne] 
builds  up  near  the  probe  (z  <  3  cm).  The  conse¬ 
quence  of  the  radial  field  is  an  azimuthal  drift  of 
the  magnetized  electrons  ( Vf  E,  x  Bo/£^  ^6x10^ 
cm/sec)  producing  the  Hall  currents  (y«= /icv«<  0.4 
A/cm^)  and  a  radial  outflow  of  the  unmagnetized  ions 
(/•„  a=  20  cm  »  fyi  =»  2  cm  »  Zj,  =«  2  mm)  causing  a 
small  radial  current  density  (/,  •»««,«=  5  mA/cm^), 


FIG.  5.  Space-time  dependence  of  electron  temperature, 
*r«,  and  density,  n,. 
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which  when  integrated  over  the  large  cylindrical  area 
involved  yields  significant  currents  —  2irr/,I„y, 
—  3.1  A).  The  Hall  current  is  expected  to  generate 
ion-sound  turbulence*  (v*»  c,,r,  »  Tj)  producing 
azimuthal-wave  electric  fields,  £«,  and  radially  inward 
electron  Bohm  transport^  (7^,  ^  /^^).  Thus  the  radial 
transport  of  ions  and  electrons  accounts  for  the  unex¬ 
pected  cross-field  currents.  The  large  parallel  electric 
field  produces  electron  heating  [or  increase  in  mean 
energy  if  /«(v,  r,/)  is  non-Maxwellian]  and  accelera¬ 
tion  of  ions  away  from  the  electrode  with  a  resultant 
density  depression.  Both  of  these  effects  are  shown  in 
Fig.  S  which  displays  kT,  and  n  vs  (r.r)  at  different 
times.  Local  heating  is  observed  during  times  and  lo¬ 
cations  of  strong  E-J  (/  «=  1.5  fisec,  z  <  3  cm).  The 
increase  in  IcT,  and  the  electron-ion  drift  velocity  pro¬ 
duces  ion  acoustic*  and  Buneman^  instabilities 
(wrf  =  1.2  X 10*  cm/sec  >  =  7  x  10^  cm/sec  »  c,, 

Tf  »  Ti )  which  have  been  tentatively  identified  in 
observed  density-fluctuation  spectra.  Such  instabilities 
produce  the  anomalous  parallel  conductivity"’  which, 
in  turn,  explains  the  large  parallel  electric  fields. 

In  summary,  currents  far  above  the  electron  satura¬ 
tion  current  can  momentarily  be  collected  in  a  magne¬ 
toplasma.  Such  large  currents  are  due  to  the  appear¬ 
ance  of  anomalous  conductivity  and  are  possible  only 
when  the  electrode  bias  is  rapidly  switched  above  the 
plasma  potential.  The  large  currents  are  terminated  by 
the  onset  of  anomalous  electric  fields  which  lead  to  ion 
expulsion,  and  hence  density  depletion,  from  the  vi¬ 
cinity  of  the  electrode.  The  observations  have  pro¬ 
found  bearing  in  the  use  of  rapidly  swept  Langmuir 
probes,  charging  of  spacecraft,  and  the  proposed  shut¬ 
tle  electrodynamic  tether  system.^ 
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